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Abstract 

This  invsstigation  sxasinss  ths  relationship  bstwssn  an  aircraft's 
lateral -directional  handling  qualities  and  its  eigenstructure. 

Historical  data  were  used  to  identify  eigenvalues  expected  to  yield 
satisfactory  handling  qualities.  Intuition  and  siaple  aatheaatical 
sodels  were  used  to  develop  desirable  eigenvectors.  For  conventional 
aircraft,  the  dutch  roll  eigenvector  was  shown  to  be  a  function  of  the 
roll  to  sideslip  (phi/beta)  ratio.  Candidate  eigenstructures,  developed 
using  the  eigenvalues  and  eigenvectors  predicted  to  yield  good  flying 
qualities  were  developed  which  varied  the  phi/beta  ratio  of  the  dutch 
roll  eigenvector. 

Control  laws  to  produce  the  candidate  eigenstructures  were  derived 
using  eigenstructure  assignsent  with  output  feedback.  The  control  laws 
were  based  on  a  12th  order  linear  nodal  of  the  lateral-directional 
dynaaics  of  the  YA-7D  DIGITAC  aircraft  at  two  flight  conditions;  0.4  and 
0.6  Mach  at  15,000  feet  pressure  altitude. 

The  handling  qualities  of  the  candidate  eigenstructures  were 
exaained  analytically  and  in  flight  testa.  Analytical  axanination  used 
sinplified  pilot  models  and  frequency  response  characteristics.  An  A-7D 
siaulator  verified  the  control  laws  produced  the  designed  open  loop 
response  and  were  safe  for  flight. 

Flight  testing  was  conducted  to  evaluate  the  lateral -directional 
handling  qualities  of  the  YA-7D  DIGITAC  at  15,000  fast  and  0.4  and  0.6 
Mach.  Specific  objectives  included:  confirn  the  pole  and  zero 
locations  of  the  lateral-directional  axes,  correlate  pilot  handling 
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quality  ratings  with  tha  phi/bsta  ratio  of  tha  dutch  roll  aoda,  and 
avaluata  aiganatructura  asaignaant  aa  a  control  law  aynthaaia  tachniqua. 
Two  Haad-up  Display  tracking  tasks  and  an  air-to-air  tracking  task  wars 
flown  to  avaluata  tha  cloaad  loop  handling  qualitias  of  tha  various 
aiganstructuras . 

Tachnical  problaas  pravantad  flight  tasting  at  0.4  Mach,  but 
aufficiant  data  wars  gatharad  at  0.6  Mach  to  aaat  all  othar  tast 
objactivas.  Tha  control  laws  darlvad  from  aiganatructura  asaignaant 
succasafully  achiavad  tha  daaignad  opan  loop  rasponaa. 

Coopar-Harpar  ratings  for  tha  cloaad  loop  taaks  ahowad  no 
corralation  with  tha  aagnituda  of  tha  phi. bata  ratio.  Howavar,  pilots 
conaiatantly  downgradad  control  laws  with  a  phasa  angla  of  120  dagraas. 
Tha  lataral  atick  dynaaica  of  tha  YA-7D  playad  an  important  rola  in  tha 
low  pilot  ratings  givan  to  tha  control  laws. 

Eiganatructura  asaignaant  was  found  to  ba  a  valid  tachniqua  for 
darivation  of  control  laws.  Tha  aasa  with  which  naw  control  laws  could 
ba  darivad  was  tha  aost  poaitiva  faatura.  Tha  major  diaadvantaga  was 
tha  tha  lack  of  intuitiva  intarpratation  of  tha  control  laws  which  aada 
problaa  idantification  and  siapla  changaa  difficult. 


AN  ANALYSIS  OF  LATERAL-DIRECTIONAL  HANDLING  QUALITIES  AND 
EIGENSTRUCTURE  OF  HIGH  PERFORMANCE  AIRCRAFT 

I.  Introduction 

Background 

It  has  long  baan  recognized  that  classical  aathods  of  control 
ayataa  analysis,  such  as  root  locus,  Boda,  and  Nyquist  tachniquaa,  hava 
several  drawbacks  whan  appliad  to  aircraft.  First,  tha  classical 
aathods  raly  on  analysis  in  tha  fraquancy  doaain,  and  tharafora  aust  ba 
conaidarad  an  indiract  approach  bacauaa  tha  tiaa  raaponsa  aust  ba 
intarpratad  froa  tha  fraquancy  raaponsa.  Also,  tha  classical  aathods 
bacoaa  cuabaraoaa  whan  appliad  to  Hultipla  Input-Multipla  Output  (MIMO) 
ayataaa.  Still,  tha  tachniquaa  hava  baan  in  uaa  a  long  tiaa  and  aany 
•pacifications,  including  aircraft  handling  qualities  design  objectives 
are  written  in  tarns  of  tha  fraquancy  raaponsa. 

Modern  control  thaory,  including  linaar-quadratic-Gausaian  (LOG) 
and  eigenatructure  aaaignaant  techniques,  use  atata  apace  nodels  which 
overcone  soaa  of  tha  liaitationa  of  tha  classical  approaches.  MIMO 
ayataaa  are  easily  handled  and  analysis  is  perforaed  directly  in  tha 
tiaa  doaain.  However,  LOG  still  suffara  froa  its  fundaaental 
aathaatical  raquiraaant  -  tha  scalar  parforaanca  criteria.  Whan  using 
LOG,  tha  designer  aust  reduce  all  tha  parforaanca  raquiraaents  to  a 
single  quadratic  criterion,  than  solve  a  vary  coaplax  nonlinear  system 
of  aquations.  While  great  strides  hava  baan  aada  in  recant  years  to 
characterize  aircraft  parforaanca  criteria  in  quadratic  fora,  LOG  still 


often  turns  out  to  ba  a  trial  and  error  approach  with  tha  designer 


adjusting  the  quadratic  parforaanca  critaria  until  a  raasonabla  solution 
is  obtainad.  Thus,  LOG  aay  also  ba  considarad  an  indiract  approach. 
Quadratic  parforaanca  critaria  rathar  than  aircraft  raaponsa  aust  ba 
spacifiad. 

Eiganatructura  assignaant  overcomes  this  last  liaitation.  MHO 
systaas  ara  aasily  handlad  and  daairad  aircraft  raaponsa  is  spacifiad 
diractly  in  taras  of  tha  closad  loop  aiganvaluas  and  aiganvactors.  Tha 
coaplex  aathaatical  aathods  nacassary  for  LOG  ara  also  aliainatad. 
Furtheraore,  tha  ordar  of  tha  closad  loop  systaa  is  not  incraasad  sinca 
aval labia  statas  ara  fad  back,  an  iaportant  considaration  in  ainiaizing 
tiaa  dalays  introducad  by  hlghar  ordar  affacta  which  dagrada  flying 
qualltiaa.  Howavar ,  tha  prica  paid  for  this  siaplicity  is  that  tha 
daairad  closad  loop  aiganvaluas  and  aiganvactors  aust  ba  spacifiad.  Tha 
daairad  aiganvaluas  ara  aquivalant  to  tha  closad  loop  polas  of  classical 
theory,  and  tharafora  a  graat  daal  of  aapirical  data  axists  which 
ralatas  aircraft  flying  qualltiaa  to  aiganvalua  placeaant  (or  pola 
location).  Unfortunately,  only  a  faw  auch  relationships  exist  for  tha 
aiganvactors. 

This  thesis  develops  approachs  to  aid  in  eigenvector  selection. 

While  tha  aaphaais  is  on  high  parforaanca  aircraft,  tha  techniques  apply 
to  other  aircraft  as  wall. 

AqaytppUPPf  gad.  Ua.4t«lUpnf 

A  nuaber  of  assuaptiona  have  bean  aada  to  siaplify  tha  coaplexity 
of  evaluating  several  possible  designs.  As  described  in  detail  later. 


tha  aathaatical  aodals  have  bean  aada  as  exact  as  possible.  Servo  and 


actuator  states,  pracoapansation  filters,  and  aanaor  dynaaica  hava  all 
baan  includad  to  aaka  tha  daalgns  aa  raaliatic  aa  poaaibla.  However,  no 
othar  attaapt  warn  aada  to  daaign  a  practical  controller.  Tha  feedback 
gains  are  baaed  on  point  designs,  using  linaarizad  aquations  of  notion 
and  linaarizad  aerodynaaics.  No  evaluation  is  aada  of  tha  sensitivity 
of  tha  controller  gains  to  changes  in  flight  condition.  As  a  result, 
gain  schedules  are  not  developed,  nor  is  tha  result  of  delating 
feedbacks  with  saall  gains  analyzad.  In  addition,  the  result  of 
failures  in  the  feedback  channels  on  vehicle  stability  is  not  discussed. 
Vhile  none  of  these  aaauaptiona  restrict  tha  results,  they  are 
considerations  that  would  be  necessary  in  the  complete  design  of  a 
practical  aircraft  controller. 

The  linearized  aerodynaaics  used  were  based  on  wind  tunnel  data  for 
a  clean  A-7D  aircraft.  The  YA-7D  D2GITAC  aircraft  used  in  the  flight 
tests  was  aodified  externally  with  a  Yaw,  Angle  of  Attack,  Pitot-Static 
Probe  (YAPS)  and  carried  alx  pylons  with  two  external  pods.  Tha  inertia 
effects  of  these  changes  were  modeled  as  closely  as  possible,  but 
corrections  to  the  aerodynamics  were  not  available.  Even  though  the 
aerodynaaics  were  not  precisely  modeled,  the  effects  were  saall  compared 
to  other  real  world  effects  and  did  not  significantly  affect  flight  test 
results.  A  more  complete  description  of  the  YA-7D  DIGITAC  aircraft  is 
given  in  Section  IX. 

The  last  assumption  applies  to  the  eigenvalues.  As  explained 
later,  they  were  picked  to  provide  the  best  flying  qualities.  However, 
there  is  generally  a  range  of  values  ovsr  which  the  flying  qualities 
will  be  the  same,  all  other  variables  being  held  constant.  Therefore, 


thera  is  sobs  flexibility  in  picking  ths  eigenvalues  which  can  be  used 
to  iaprovs  ths  "fit"  of  ths  associated  eigenvector.  This  procedure  is 
discussed  in  detail  in  Reference  2,  but  here  would  introduce  one 


f- 


additional  variable  in  the  atteapt  to  correlate  handling  qualities  and 
eigenspace  selection  and  therefore  has  been  avoided.  But  again,  in  an 
actual  aircraft  control  systea  design,  this  flexibility  could  be  used  to 
iaprovs  the  systea  perforaanca. 


The  principles  of  eigenstructure  aaaignaant  and  tha  difficulty  of 


aalacting  an  aiganatructura  to  yiald  good  flying  qualitiaa  ia  baat 
explained  rav laving  linaar  ayataa  thaory.  Conaidar  tha  ayataa  defined 
by 

a 

x  *  Ax  ♦  Bu  (1) 

where  x  ia  a  n  X  1  atata  vector 

A  ia  a  n  X  n  conatant  aatrix 
B  ia  a  n  X  a  conatant  aatrix 
u  ia  a  a  X  1  control  vector 

Tha  aolution  to  thia  linaar  autonoaoua  ayataa  ia 

x<t)  »  e*tx<0)  ♦  \  eA*Bu<t-T)dT  <2> 

u 

where  x<0>  ia  tha  initial  condition. 

eAt  is  the  atata  tranaition  aatrix  and  can  be  solved  in  teraa  of 
tha  ayataa  aiganvaluaa  and  eigenvectora.  Aaauaing  n  unrapaatad 
eigen valuea: 


.At  «  TeAtT-1  (3> 

where  A  ia  a  n  X  n  diagonal  aatrix  of  ayataa  aiganvaluaa,  \ i 

aAt  ia  an  n  X  n  diagonal  aatrix 

T  ia  tha  n  X  n  aodal  aatrix  of  ayataa  eigenvectors 
Making  thia  aubatitution,  tha  aolution  to  Eq  (1)  becomes: 


x<t> 


TaAtT-lx(O) 


T  ^  •AtT-lBu(t-r)dT 


<4> 


Equation  (4)  illuatrataa  tha  fundaaantal  idaas  of  aigenstructura 
aaaignaant.  Cloaa  axaaination  ravaala  tha  following: 

(1)  Tha  aiganvaluaa  dateraina  tha  dacav/orowth  rata  of  both 
tha  fraa  raaponaa  and  tha  for cad  rasponsa. 

(2)  Tha  aiganvactora  dataraina  tha  diatribution  or  ahapa  of 
tha  raaponaa. 

(3)  Tha  fraa  raaponaa  dapanda  on  tha  initial  condition,  x<0), 
which  in  turn  datarainaa  tha  dagraa  to  which  aach  aoda 
will  participata.  Tha  axact  aodal  diatribution  ia 
coaplicatad  by  tha  tara  T_1x<0). 

(4)  Tha  f oread  raaponaa  ia  of  couraa  dapandant  on  tha  control 
input  u(t).  Howavar,  tha  control  input  diatribution  aaong 
tha  aodaa  ia  difficult  to  pradict  bacauaa  of  tha  tara 

t-1b«.. 

Thua,  it  ia  avidant  that  if  faadback  ia  to  ba  uaad  to  iaprova 
aircraft  handling  qualitiaa,  aiganvactor  aalaction  aa  wall  aa  pola 
(aiganvalua)  placaaant  auat  ba  conaidarad.  It  ia  alao  claar  that 
aalaction  of  "optimum"  aiganvactora  to  raault  in  good  flying  qualitiaa 
will  ba  difficult  bacauaa  of  tha  T-1  tarma  in  tha  raaponaa.  Pilot 
control  inputa  will  ba  diatributad  among  tha  aodaa  by  tha  tara  T_1B, 
furthar  clouding  afforta  to  ralata  flying  qualitiaa  and  aiganvactor 


aalaction 


As  an  example  consider  the  lateral -directional  dynamics  of  the  A-7D 
at  0.4  Mach  described  by: 

x  *  Ax  ♦  Bu  (1) 

where  x  *  t  r  0  p  0)? 
u  *  [£as  $r  ^ 


A  * 


"  -.4652  3.732  .1415  0 

-.9978  -.1216  .0014  .0761 

1.082  -14.530  -1.944  0 

0  0  1.000  0 


(5a) 


B 


_  1.258  -2.719 

.006  .034 

-15.85  2.394 

0  0 


(5b) 


One  simple  measure  of  lateral-directional  flying  qualities  is  the 
amount  of  pilot  rudder  input  necessary  to  coordinate  a  turn.  That  is, 
large  sideslip  angles  resulting  from  an  aileron  input  will  usually 
degrade  pilot  opinion.  If  no  rudder  is  necessary  the  turn  is  perfectly 
coordinated.  Clearly  then,  it  is  desirable  to  eliminate  the  sideslip 
response  as  much  as  possible  to  improve  turn  coordination.  In  addition, 
without  justification  at  this  time,  the  dutch  roll  response  will  be 
decoupled  so  it  consists  of  a  pure  yawing  motion.  The  methods  described 
in  the  next  section  can  be  used  to  calculate  a  feedback  matrix,  F,  such 
that  u  -  Fx  will  produce  the  desired  eigenstructure.  The  gains  are 


F  =  f  .1762  -1.287  .0488  -.0055  1  (6) 

.7149  -2.452  .0804  -.0628  J 


After  augmentation  the  system  modal  characteristics  are: 


jr  ’jT’jrrjr^rr, 


1.000  ♦  01 
.112  + .3131 
0 
0 


1.000  -  01  1.0  1.0 

.112  -.3131  .0  -2.49 

0  -.329  -264.93 

0  13.17  105.96 


A  =  dlagl  -1.2  ♦  2.751,  -1.2  -  2.751,  -.025,  -2.50  ) 


The  eigenvalue*  show  a  well  damped  dutch  roll  and  atable  roll 
aubaldence  and  spiral  aodes.  The  corresponding  eigenvectors  show  the 
desired  decoupling  of  roll  response  In  the  dutch  roll  and  the 
elimination  of  any  sideslip  response  in  the  spiral  aodes.  For  reasons 
explained  later,  there  is  still  a  saall  coaponent  of  sideslip  and  yaw 
rata  in  the  roll  subsidence  eigenvector. 

The  sideslip  response  to  a  unit  step  aileron  input  is: 


0<t)  *  .186  *-l*2tmin<2.75t  ♦  50.633)  ♦  .061  e~2.5t  -  .202  (8) 


Thus,  although  nearly  all  sideslip  response  was  eliminated  from  the 
roll  subsidence  and  spiral  aodes,  the  turn  was  not  perfectly 
coordinated.  Even  though  the  free  response  would  not  show  any  sideslip 
response  to  an  initial  roll  rate,  the  forced  response  produces  a  steady 
state  sideslip.  As  stated  earlier,  this  is  due  to  the  control  nixing 
occurring  through  the  T"lBu  tera.  Physically,  the  adverse  yaw 
generated  by  the  step  aileron  input  excites  the  dutch  roll  mode.  It 
illustrates  that  augaenting  the  three  basic  aircraft  aodes  alone  nay  not 
be  sufficient  to  achieve  the  desired  flying  qualities.  Additional 
aodification  of  the  control  inputs  may  also  be  necessary. 


Elgenatructure  Assignment  using  Output  Feedback 


The  mathematical  baaia  for  eigenstructure  assignment  using  output 
feedback  ia  explained  in  detail  in  Reference  3.  The  purpoae  here  ia  to 
explain  the  fundamentals,  without  proof,  in  a  clear  underatandable 


Consider  the  linear  time  invariant  ay stem  described  by  the 
equations: 


x ( t )  *  Ax(t)  ♦  Bu(t) 


where 


y<t)  *  Cx(t) 


i)  |  t  Rn  ,  u  t  8*  >  It  (:  Dr  are  the  state,  control,  and 
output,  respectively; 

ii)  A,  B,  C  are  real  constant  matrices; 

iii)  rank  A  *  n  ,  rank  B  *  a  ,  rank  C  *  r; 

iv)  the  system  is  fully  controllable  and  fully  observable. 


The  linear  control  law  ia  of  the  form: 


u  ( t)  ■  Fy(t)  (10) 

It  can  be  shown  the  number  of  closed  loop  eigenvalues  that  can  be 
assigned,  and  number  of  asaociated  eigenvectors  that  can  be  partially 
assigned  is  equal  to  the  larger  of  a  or  r.  The  number  of  entries  in 
each  vector  exactly  chosen  is  the  smaller  of  a  or  r.  Since  for  typical 
aircraft  control  aysteas  r  is  greater  than  a,  to  increase  the  number  of 
assigned  eigenvalues  the  number  of  outputs  must  be  increased. 

Similarly,  to  increase  the  number  of  entries  in  each  eigenvector  that 
are  specified  (or  improve  the  "fit"),  the  number  of  independent  control 


input*  aust  be  increased.  It  follows  that  with  a  aingla  control  input 
the  designer  has  no  control  over  the  eigenvectors.  Substituting  Eq  (10) 
into  Eq  (1)  yields: 

• 

x(t)  •  (A  ♦  BFC)x(t)  (11) 

Assuae  (Ai)  are  the  r  desired  closed  loop  eigenvalues  and  vj.  are  the 
correaponding  closed  loop  eigenvectors.  Then  it  follows  they  aust 
satisfy: 

(A  ♦  BFC)yi  «  AiVi  (12) 

or 

yi  «  (Ail  -  A)-lBFCvi  (13) 

Let  wi  3  FCvi.  Then  after  substitution 

(Ail  -  A)‘*Bwi  •  yi  (14) 

The  iaplication  of  Eq  (14)  is  very  iaportant.  It  shows  the  desired 
closed  loop  eigenvector  Vi  aust  lie  in  the  aubspace  spanned  by  the 
coluana  of  (Ail  -  A)-1B.  »»i  represents  the  linear  coabination  of  the 

coluans  that  results  in  yi.  The  subspace  is  diaension  a  which  is  the 
rank  of  B,  the  nuaber  of  independent  control  variables.  Therefore,  the 
nuaber  of  control  variables  deteraines  how  large  the  subspace  is  in 
which  achievable  closed  loop  eigenvectors  aust  reside.  Thus,  as  stated 
earlier,  a  coaponents  of  the  desired  eigenvector  can  be  specified 
exactly.  And  it  can  be  concluded  that  if  the  eigenvectors  yi  are  chosen 
to  lie  precisely  in  the  subspace  spanned  by  the  coluans  of  (Ail  -  A)~lB, 
they  will  be  achieved  exactly. 

In  general,  however,  the  desired  eigenvectors  will  not  lie 
precisely  in  the  subspace,  and  therefore  cannot  be  exactly  achieved.  In 
other  words,  aore  than  a  coaponents  of  the  eigenvector  are  specified. 


An  achievable  eigenvector  is  chosen  as  the  projection  of  vj.  onto  the 
subspace  spanned  by  the  coluans  of  (K^I  -  A)-1B  to  ainiaize  the  sua  of 
the  squares  of  the  error.  This  is  aost  easily  accoaplished  as  follows: 

Let  Li  *  (Xj.1  -  A)-lB  (15) 

Then  the  achievable  eigenvector  is 

via  *  LiL*iVi  (16) 

where  L*i  is  the  generalized  inverse  of  Li.  It  is  best  to  calculate  the 
generalized  inverse  using  a  singular  value  decoaposition  to  avoid 
ill-conditioning  difficulties  (4) . 

This  procedure  also  allows  partial  specification  of  the 
eigenvector.  It  would  be  difficult  to  specify  each  coaponent  of  every 
eigenvector.  Instead,  only  those  coaponents  which  the  designer  wishes 
to  specify  and  their  corresponding  equetions  in  Li  are  used  in  the  least 
squares  solution  of  yj,a> 

With  via  calculated  it  is  easy  to  calculate  WJ,  using  Eq  (14).  wi 
is  then  used  to  find  the  unspecified  coaponents  in  Yia.  It  follows 
directly  froa  the  definition  of  wi  that 

F  =  W(CV)“1  (17) 

where  V  is  aade  up  the  eigenvectors  that  lie  in  the  achievable  subspace 
of  (Xj.1  -  A)~lB.  The  inverse  of  CV  will  exist  provided  the  desired 
eigenvectors  are  linearly  independent  and  the  aultiplicity  of  each 
desired  eigenvalue  is  less  than  or  equal  to  the  saaller  of  Rank  C  or 
Rank  B.  For  those  instances  where  the  inverse  does  not  exist,  there  is 
usually  a  design  flaw  such  as  noncontrollability. 
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While  output  feedback  doee  not  offer  the  design  flexibility  of  full 
state  feedback,  it  is  sore  practical  to  implement,  eliminating  the  need 
for  observers  to  estimate  those  states  that  cannot  be  measured. 
Furthermore,  notch  filters,  integrsl  compensation,  and  high  and  low  pass 
filters  can  be  included  in  the  closed  loop  design.  Feedforward  gain 
compensation  to  modify  the  control  inputs  can  be  left  for  post  design. 

Eigenstructure  assignment  using  output  feedback  does  suffer  fron  a 
couple  of  drawbacks.  First,  stability  is  not  guaranteed  for  those 
eigenvalues  not  assigned.  That  is,  the  free  poles  can  move  to 
unpredictable  and  sometimes  undesirable  places.  If  they  do  move  to  an 
undesirable  location,  there  is  no  set  method  to  vary  the  desired 
eigenstructure  to  move  them  to  more  satisfactory  positions.  Also,  thera 
is  no  penalty  on  the  control  effort  required.  The  desired 
eigenstructure  may  be  ochieved  at  the  cost  of  saturating  control 
deflections  or  rates.  While  these  are  legitimate  design  concerns,  they 
did  not  cause  difficulties  in  this  research.  Methods  exist  (3)  to 
overcome  these  drawbacks,  but  full  state  feedback  (or  estimation)  is 
required. 

The  procedure  described  above  was  programmed  in  FORTRAN  and  used  in 
the  design  of  the  control  laws  used  throughout  this  report.  A  copy  of 
the  program  and  a  description  of  the  implementation  (including  the 
Incorporation  of  complex  eigenvalues  and  eigenvectors)  is  included  in 
Appendix  A. 

Appendix  B  contains  a  sample  output  of  the  eigenstructure 


assignment  program 


III.  YA-7D  Linear  Hpdii 


Twelfth  order  aodels  were  used  to  describe  the  eyetea  dynaaice  as 
accurately  ae  possible.  The  linear  aodels  used  in  the  analysis  were 
carefully  developed  to  aatch  the  aircraft  dynaaics  closely  so  variations 
in  handling  qualities  due  to  changes  in  the  eigenstructure  would  not  be 
confused  with  aodeling  errors.  Higher  order  effects  with  frequencies 
less  than  100  radians/ second  (rad/sec)  were  included  in  the  aodel.  This 
cutoff  frequency  was  arbitrary,  but  insured  the  aircraft  response  in  the 
frequency  range  of  interest  was  adequately  aodeled.  Thus,  instead  of  a 
siaple  fourth  order  aodel  to  siaulate  the  lateral -directional 
characteristics  of  the  YA-7D,  a  12th  order  aodel  was  necessary.  In 
addition  to  the  four  basic  aircraft  rigid  body  states,  actuator  and 
servo  states  for  the  aileron  and  rudder  are  included,  adding  a  total  of 
six  additional  states.  High  pass  and  noise  filters  were  also  added  in 
the  precoapensation  design  of  the  control  systea  for  a  total  of  12 
states.  The  following  discussion  explains  the  developaent  of  the  state 
space  aodel  in  detail. 

Rigid  Body  Equations  of  Motion 

The  linear  perturbation  lateral -directional  aodel  used  in  this  report 
was  derived  froa  the  nonlinear  equations  of  notion  (6:259).  These 
linear  equations  assuae: 

i)  Stability  axes  are  used  with  the  X  axis  initially  aligned  with 
the  horizon. 

ii)  The  perturbations  froa  equilibriua  are  saall. 

iii)  The  aircraft  is  a  rigid  body  of  constant  nass. 


iv)  The  flow  la  quasi steady  and  the  aerodynaalca  can  be  described 
by  linear  terns. 

The  equations  of  notion  of  the  rigid  body  in  state  space  forn  are: 


x  =  Ax  ♦  Bu 

where  x  =  [  r  fi  p  $  ] T 
u  *  C  6as  fr 


A 


B 


where: 


Mr' 

MB' 

Mp' 

0  “ 

Yr/Uo-l 

Yv 

Yp/Uo 

g/Uo 

W' 

L0' 

Lp' 

0 

0 

0 

1 

0  . 

Mias' 

M6r' 

- 

YA««/Uo 

YAr/Uo 

LA  as' 

Ur' 

0 

0 

r  -  yaw  rate,  rad/sec 

0  -  sideslip  angle,  rad 

p  -  roll  rate,  rad/sec 

0  -  bank  angle,  rad 

*as  “  aileron/ spoiler  deflection,  rad 

6 r  ~  rudder  deflection,  rad 


(1) 


(18a) 


<18b) 


The  prine  on  the  aerodynaaic  derivatives  refers  to  a  conbined 
stability  derivative  that  includes  product  of  inertia  effects  (6:257). 
For  a  general  notion  or  input  quantity  designated  by  the  subscript  i. 


Li  ♦  (IXz/Ix>  Mi 

Li'  «  -  (19) 

1  -  (IX22/ixiz) 

Mi  ♦  <Ixz^z> 

Mi'  *  -  (20) 

1  -  (Ixz^/Ixlz* 


In  the  reaainder  of  the  report  the  prine  will  be  dropped  for 
convenience.  Except  where  noted  otherwise,  any  reference  to  a  stability 
derivative  refers  to  the  priaed  derivative. 


The  aarodynaaic  data  wara  taken  fro*  Raferancaa  7  and  8.  Two 
flight  conditiona  wara  chosen  to  investigate:  one  g  level  flight  at 
15,000  feat  pressure  altitude  (He)  and  0.40  and  0.60  Mach.  The  YA-7D 
exhibits  relatively  poor  lateral -directional  characteristics  at  0.4  Mach 
with  eoderate  adverse  yaw  and  poor  dutch  roil  characteristics.  The 
second  flight  condition,  0.6  Mach,  represents  a  sore  benign  condition 
where  the  YA-7D  has  little  adverse  yaw  and  slightly  improved 
lateral-directional  characteristics. 

The  aerodynaaic  data  was  for  a  clean  A-7D.  The  aircraft  used  in 
the  flight  tests  was  the  YA-7D  DIGITAC  described  in  Section  IX.  The 
primary  external  differences  were  the  addition  of  a  YAPS  head,  six 
pylons,  a  DIGITAC  data  pod  under  the  right  wing,  and  an  inert  MK-82  boeb 
under  the  left  wing.  Unfortunately,  exact  inertia  and  aerodynaaic  data 
for  this  configuration  were  not  available.  Mass  properties  froa 
Reference  7  were  corrected  for  the  pylons,  data  pod,  and  MK-82  and  were 
used  with  the  nondiaensional  aerodynaaic  coefficients  to  calculate  the 
diaensional  stability  derivatives.  No  corrections  were  aade  to  the 
aerodynaaic  coefficients  for  the  external  vehicle  differences  and 
reference  angle  of  attack  increase  (because  of  the  increased  weight  of 
the  YA-7D  DIGITAC). 

The  diaensional  stability  derivatives,  vehicle  aaas  properties,  and 
flight  conditions  are  Hated  in  Table  I  and  Table  II  on  the  following 


pages 


Table  I 


YA-7D  Aerodynaaic  and 

Vehicle  Data 

(0.4  Mach) 

Mach  No. 

.40 

Weight 

27,180  lbs 

Altitude 

15,000  ft 

lx 

20,708  slug  ft? 

Load  Factor 

1.0 

*z 

83,020  slug  ft? 

Configuration 

Clean 

Ixz 

-7,412  slug  ft? 

C.G.  x  MAC 

28.71 

a  Tria 

9.55  Deg 

Nr' 

-0.4260 

1/sec 

NB' 

3.349 

l/sec2 

*p' 

0.1592 

1/sec 

Was' 

0.9953 

1 /rad-sec? 

Mr' 

-2.541 

l/rad-sec? 

Lr' 

0.8392 

1/sec 

L0' 

-11.25 

l/sec2 

. 

Lp' 

-1.516 

1/sec 

Ldas' 

-12.23 

l/rad-sec2 

L4r* 

1.882 

l/rad-sec? 

T  r 

0.0021 

1/rad 

Yv 

-0.1133 

1/sec 

Yp 

0.0013 

1/rad 

Yda. 

0.0058 

1 /rad-sec 

Ydr 

0.0316 

1 /rad -sec 

g/Uo 

.0761 

1/sec 

Actuator  and  Servo  States 

Tha  aileron  and  rudder  actuators  have  bandvidtha  of  20  rad/aec  and 
33.33  rad/aec,  respectively.  They  are  driven  by  identical  94  rad/sec, 

(  *  0.7  servos.  In  addition,  the  YA-7D  uses  an  aileron-rudder 
interconnect  to  reduce  the  adverse  yaw  due  to  an  aileron/spoiler 
deflection.  It  ia  aechanized  through  the  autoaatic  flight  control 
aystea  by  feeding  back  aileron  deflection  (aeaaured  at  the  roll  servo) 
to  the  rudder  channel.  The  croaafeed  gain  ia  noraally  a  function  of 
horizontal  tail  deflection  and  at  0.4  Mach  and  15,000  feet  has  a  value 
of  0.20.  Because  of  the  saall  aaount  of  adverse  yaw  at  0.6  Mach,  the 
croaafeed  gain  was  set  to  zero.  No  atteapt  was  aade  to  further  optiaize 


Table  II 


YA-7D  Aerodynaaic  and 

Vehicle  Data 

(0.6  Mach) 

Mach  No. 

.60 

Weight 

27,180  lbs 

Altitude 

15,000  ft 

lx 

19,900  slug  ft2 

Load  Factor 

1.0 

Iz 

83,828  slug  ft? 

Configuration 

Clean 

Ixz 

-1,978  slug  ft? 

C.G.  *  MAC 

28.71 

a  Tria 

4.63  Deg 

Mr' 

-0.4979 

1/aec 

N*' 

4.243 

l/sec? 

Mp' 

0.0485 

1/sec 

Mas' 

-0.1145 

l/rad-sec? 

*Ar' 

-4.889 

1 /rad -sec? 

tr' 

0.7462 

1/sec 

Ul' 

-20.30 

l/sec2 

Lp' 

-2.323 

1/sec 

Lias' 

-26.54 

l/rad-sec? 

Lir' 

4.662 

l/rad-sec2 

Tr 

0.0020 

1/rad 

Yy 

-0.1507 

1/sec 

Yp 

0.0008 

1/rad 

YAas 

0.0105 

1/rad-sec 

Yir 

0.0422 

1/rad-sec 

g/Uo 

.05072 

1/sec 

the  gains.  The  croaafeed  can  be  added  directly  In  the  aervo  equations 
and  the  state  space  aodel  of  the  servos  and  actuators  becomes: 


where 


x  »  Ax  ♦  Bn 


(1) 


*  *  t  <5  as 

Ar  V1 

v2  ^3 

V4  ]Y 

u  »  I  A  sac 

Arc 

(the  coaeanded 

deflections) 

-  -20 

0 

20 

0 

0 

0  " 

0 

-33.33 

0 

0 

33.33 

0 

* 

0 

0 

0 

1 

0 

0 

0 

0 

-8883 

-131 

.9  0 

0 

0 

0 

0 

0 

0 

1 

_  0 

0 

1777 

0 

-8883  -131 

.9  _ 

(21a) 
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B 


0  0 

0  0 

8883  0 

0  0 


0  8883  J 


(21b) 


where  V1...V4  -  serv o  states,  rad 

Aasc  -  cossanded  aileron/spoiler  deflection,  rad 
i rc  "  cossanded  rudder  deflection,  rad 


Sensor  Dvnaslcs 

The  YA-7D  DIGITAC  Digital  Multimode  Flight  Control  Systes  (DFCS) 
had  four  quantities  available  for  feedback;  body  roll  and  yaw  rates, 
roll  angle  fros  the  Inertial  Measuresent  Systes  (IMS),  and  lateral 
acceleration  seasured  by  dual  acceleroseters  mounted  4.78  feet  forward 
of  the  center  of  gravity  (CG) .  Therefore,  four  sensors  have  dynamics 
which  say  have  to  be  considered  in  the  state  space  model;  roll  rate 
gyro,  yaw  rate  gyro,  roll  angle  inertial  measuresent  unit,  and  lateral 
acceleroaeter .  The  roll  rate  gyro,  yaw  rate  gyro,  and  lateral 
acceleroaeter  have  second  order  dynaaics  with  frequencies  of  314 
rad/sec,  126  rad/sec,  and  100  rad/sec,  respectively.  These  are  all  high 
enough  that  any  errors  induced  by  ignoring  their  dynanica  are 
negligible.  The  IMS  sensing  roll  angle  also  has  second  order  dynaaics, 
but  at  a  such  lower  natural  frequency,  31.4  rad/sec.  However, 
experience  showed  these  gains  to  be  relatively  aaall,  so  the  dynamics 
will  also  be  ignored. 

While  sensor  dynaaics  can  be  ignored,  sensor  position  corrections 
cannot.  The  roll  rate  and  yaw  rate  gyros  aeasure  rates  in  body  fixed 
axes.  But  the  state  space  aodel  is  in  stability  axes,  which  differ  fros 
body  fixed  axes  by  the  equilibriua  angle  of  attack.  Therefore,  the 
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aaaaured  roll  rate  will  consist  of  coaponanta  of  both  roll  and  yaw  rate 
as  aaaaured  In  stability  axes.  Siailarly,  tha  aaaaured  yaw  rata  will 
consist  of  coaponsnts  froa  both  ratas.  In  addition,  the  yaw  rate  gyro 
is  tilted  down  2.5  degrees.  The  transforaation  is: 


[pg  "I  I"  cos  a  -sin  a  1  f  P  1 

rg  J  L  sin(a  -  2.5)  cos(a  -  2.5)  J  L  r  J 


where:  pg  -  roll  rate  gyro  output,  rad/ sec 

p  -  stability  axis  roll  rate,  rad/sec 
rg  -  yaw  rata  gyro  output,  rad/sec 
r  -  stability  axis  yaw  rate,  rad/sac 

The  correction  to  the  aaaaured  roll  rata  is  applied  directly  in  the 

output  aatrlx  C  using  Eq  (22).  However,  since  the  yaw  rate  is  filtered 

prior  to  feedback,  the  correction  is  applied  differently.  This  will  be 

explained  in  aore  detail  in  the  next  section. 

The  lateral  acceleroaeter  incurs  significant  error  because  it  is 
not  located  at  the  canter  of  gravity,  but  4.78  feet  forward.  The 
acceleroaeter  is  coaaonly  placed  forward  of  tha  center  of  gravity  so  the 
neasured  acceleration  is  aore  of  a  true  representation  of  the  sideslip. 
It  picks  up  a  significant  contribution  froa  the  yaw  acceleration  which 
will  tend  to  cancel  the  contribution  to  lateral  acceleration  froa  the 
rudder,  leaving  sideslip  the  doainant  contributor  (6:483-487).  A 
siailar  correction  because  of  the  use  of  stability  axes  aust  also  be 
aade  here  as  in  the  section  above.  The  equation  for  the  lateral 
acceleration  at  a  distance  X0  forward  of  the  center  of  gravity  and  Z0 
below  it  is 


ay  *  Uq0  -  UqT  ♦  X0r  -  ZqP  -  g# 


In  body  axes  the  YA-7D  lateral  acceleroaeter  is  4.78  feet  forward  of  the 
CG.  Therefore: 


Xo  *  4.78  cos  a 

Zo  a  -4.78  sin  a 


(24) 


The  development  of  this  portion  of  the  state  apace  sodel  is  not  yet 
complete,  however,  since  a  noise  filter  must  be  added  to  the 
accelerometer  output.  This  will  be  explained  in  more  detail  in  the  next 
section. 


Predeslqn  Comoensetion 

A  “high  pass“  filter  was  added  in  the  yaw  rate  feedback  channel. 

The  transient  yaw  rata  must  be  isolated  from  the  steady  state  response 
so  the  control  system  will  not  fight  the  steady  state  yaw  rate 
associated  with  a  turn.  This  is  classically  done  by  adding  a  “high  pass" 
or  washout  filter.  The  transfer  function  is: 


rwo  • 

-  - -  (25) 


rg  s  ♦  1 

where  r*o  is  the  washed  out  yaw  rate  and  rg  is  the  yaw  gyro  output. 

This  filter  will  have  the  effect  of  passing  signals  with  frequencies 
higher  than  one  rad/sec  while  suppressing  lower  frequency  signals,  thus 
achieving  the  seired  effect.  Its  implementation  in  the  state  space 
model  is  complicated  slightly  because  the  rate  filtered  is  the  yaw  rate 
sensed  by  the  gyro,  not  the  stability  axis  yaw  rate  used  in  the  model. 
As  discussed  in  the  previous  section,  the  gyro  yaw  rate  will  consist  of 


a  linear  combination  of  the  stability  axis  yaw  rate  and  roll  rate.  When 


this  is  tsksn  into  consideration,  ths  resulting  differential  equation  of 
the  washout  filter  is 


rwo  =  cos<a  -  2.5)  r  *  sin(a  -  2.5)  p  -  rwo  (26) 

When  the  necessary  cosputations  are  aade,  the  washout  filter  in  state 
space  fora  at  0.4  Mach  becomes: 

« 

x  *  Ax 


x  »  t  r  0  p  0  <Sas  i  r  ^wo  3T 


A  = 


-.3198  1.943  -.0281  0  -.5133  -2.291  -1 


(27) 


At  0.6  Mach,  the  state  space  equation  becoaes: 


A  » 


-.4698  3.486  -.0379  0  -1.101  -4.712  -1 


(27b) 


For  siaplicity,  only  the  entries  for  the  washout  filter  have  been 
included  in  the  A  aatrix.  The  other  entries  have  been  discussed 
elsewhere  and  are  not  changed  by  the  addition  of  the  high  pass  filter. 

To  eliainete  high  frequency  noise  froa  the  lateral  acceleration 
signal,  a  low  pass  or  "noise  filter”  was  necessary.  The  break  frequency 
of  the  filter,  10  rad/sec,  was  chosen  to  eliminate  aost  unwanted  noise 
without  having  a  great  affect  on  the  low  frequency  dynamics.  The  final 


fora  of  the  filter  becoaes: 
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(28) 


s  ♦  10 


where  af  la  the  filtered  acceleroeeter  signal.  The  atate  apace 
representation  at  0.4  Mach  of  the  filtered  lateral  acceleration  aignal 


corrected  for  position  becoaea: 


x  *  Ax 


where  x  ■  t  r  0  p  0  6a a  <*r  af  ]T 


■17.53  -245.11  22.50  0.  147.34  -4.62  -10 


At  0.6  Mach,  the  atate  space  representation  becoaea: 


A  «  -19.59  -539.88  31.39  0.  332.48  -20.91  -10 


(29a) 


All  the  necessary  portions  of  the  YA-7D  state  apace  aodel  have  now 
been  developed.  The  resulting  12th  order  aodel  is  given  in  Table  III 
and  Table  IV  on  the  following  pages  for  0.4  and  0.6  Mach,  respectively. 
A  spectral  decoaposition  of  the  aodels  ia  included  in  Appendix  B. 


WSVI 


Tabla  III 


YA-7D  Stata  Spaca  Modal  <0.4  Mach) 

x  «  A*  ♦  Bu  y.  *  Cx 
x  *  t  r  0  p  0  tfr  vi 

v2  v3  v4  rwo  af  ]T 

u  *  [  iaec  <rc 
X  *  t  r^o  «f  M  •  1T 


.4260 

3.349 

0.1592 

0 

0.9953 

-2.541 

.9979 

-.1133 

0.0013 

0.076 

0.0058 

.03159 

0.8392 

-11.25 

-1.516 

0 

-12.23 

1.882 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

-20 

0 

0 

0 

0 

0 

0 

-33.33 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

.3198 

1.943 

-.0281 

0 

-.5133 

-2.291 

17.53 

-245.1 

22.50 

0 

147.34 

-4.623 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

33.33 

0 

0 

1 

0 

0 

0 

-131.9 

0 

0 

0 

0 

0 

1 

0 

0 

-8883 

-131.9 

0 

0 

0 

0 

-1 

0 

0 

0 

0  -1 

0 

,  0 
0 
0 

20 

0 

0 

-8883 

0 

1777 

0 

0 


-r.  -v- -r  . 
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YA-7D  Stata  Spaca  Modal  (0.6  Mach) 


&  *  t  r  0  p 


u  *  C  6aac  ire 
H  *  I  r„o  «f 


"-.4979 

-.9980 

0.7462 

0 

0 

A«  0 

0 
0 
0 
0 

-.4698 
_ -19.59 


x  ■  Ax  ♦  Bu  y,  =  Cx 


0  Aaa  *1 


v2  v3  v4  Two  •£  1T 
)T 


pg  M  )T 


4.243 

0.0485 

0 

-0.1145 

-4.889 

0 

-.1507 

0.0008 

0.0507 

0.0105 

.04219 

0 

20.30 

-2.323 

0 

-26.54 

4.662 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

-20 

0 

20 

0 

0 

0 

0 

-33.33 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

-8883 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1777 

3.486 

-.0379 

0 

-1.101 

-4.712 

0 

539.9 

31.39 

0 

332.54 

-20.91 

0 

0 

0 

0 

0 

0  " 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

33.33 

0 

0 

0 

1 

0 

0 

0 

0 

-131.9 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

-8883 

-131.9 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

0 
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Fortunately,  there  is  a  great  deal  of  empirical  data  to  aeeiat  in 
the  selection  of  the  desired  closed  loop  eigenvalues.  Military 
Specification.  Flying  Qualities  of  Piloted  Vehicles  -  MIL-F-8763C  <9) 
and  the  background  information  in  Reference  10  have  extensive  charts 
relating  the  eigenvalues  of  the  three  basic  aircraft  lateral-directional 
modes  to  pilot  opinions.  In  the  following  discussion,  the  eigenvalues 
are  picked  to  yield  Level  1  flying  qualities  (equivalent  to  less  than 
3.5  on  the  Cooper-Harper  scale).  They  were  chosen  for  Class  IV  aircraft 
(fighter-  interceptor,  attack,  etc)  and  Category  A  Flight  Phase 
(air-to-air  combat,  ground  attack,  etc).  These  are  the  most  stringent 
criteria  and  apply  to  most  of  the  A-7D  missions. 

Roll  Subsidence  Mode 

Because  of  its  importance  to  flying  qualities,  there  is  a  large 
empirical  data  base  relating  this  eigenvalue  to  pilot  opinion.  The 
selection  is  a  little  clouded  because  other  factors,  such  as  the  stick 
force  shaping,  aileron  power  (  Lias  ) *  and  adverse  yaw  (  Ntfas  )  ell  play 
a  role  in  the  optimum  value.  Prefilters  in  the  YA-7D  roll  Command 
Augmentation  System  (CAS)  slow  the  response  and  must  also  be  considered. 
However,  it  is  clear  there  are  upper  and  lower  bounds  on  this 
eigenvalue.  If  it  is  too  large,  the  roll  response  will  likely  be  rated 
too  sensitive  and  abrupt  and  can  lead  to  "roll  ratcheting,"  a  type  of 
pilot  induced  oscillation.  If  the  value  is  too  small,  the  response  will 


be  "sluggish". 


Paragraph  3. 3. 1.2  of  HIL-F-8785C  raquiraa  a  minimum  value  of  1.0 
for  the  roll  mode  eigenvalue  (9:23).  Although  a  maximum  value  is  not 
included  in  current  requirements,  the  supporting  data  suggests  that  when 
the  roll  mode  eigenvalue  is  greater  than  3.0  there  is  a  constant 
degradation  in  pilot  opinion.  Several  examples  of  eigenstructure 
assignment  in  the  literature  (2  and  5)  have  used  4.0  for  the  roll 
subsidence  node  eigenvalue,  but  this  is  likely  too  high.  A  value  of  2.5 
is  near  optimum. 

Because  of  the  prefilters  in  the  CAS,  a  faster  time  constant  must 
be  chosen.  With  this  consideration,  the  optimum  value  for  this 
eigenvalue  is  chosen  to  be: 

Xri  ■  “4.0  (30) 


Spiral  Mode 

The  spiral  eigenvalue  has  little  direct  effect  on  flying  qualities. 
MIL-F-8785C  does  not  even  require  stability,  specifying  a  minimum  time 
to  double  of  12  seconds  (9:23).  While  there  are  some  data  that 
indicates  there  should  be  a  limit  on  the  degree  of  positive  spiral 
stability,  other  data  show  strong  positive  spiral  stability  can  be 
beneficial.  However,  the  degree  of  positive  stability  is  indirectly 
limited  by  the  limit  on  aileron  forces  in  turns  and  the  required  roll 
maneuverability.  Thus,  from  a  handling  qualities  point  of  view,  the 
selection  of  this  eigenvalue  is  arbitrary  and  it  could  be  argued  that 
there  is  no  point  in  expending  the  control  effort  to  place  it.  However, 
the  elimination  of  sideslip  from  the  spiral  nodal  response  is  crucial  to 


a  coordinated  turn.  In  this  case,  the  eigenvalue  is  not  critical,  but 


the  eigenvector  ia.  Therefore,  the  spiral  eigenvalue  ia  arbitrarily 
choaen  to  be: 

*  -.025  (31) 

This  corresponds  to  a  tiae  to  half  of  28  seconds. 

Dutch  Roll  Mode 

Allowable  dutch  roll  oscillations  are  specified  in  teraa  of  ainiaua 
values  of  the  daaping  ratio  (Id**  natural  frequency  (wd>,  and  td*d:  the 
last  of  which  is  a  function  of  the  roll  angle  to  sideslip  ratio,  1 0/p  I 
when  that  ratio  is  large.  The  ainiaua  racoaaended  values  froa  paragraph 
3. 3. 1.1  of  HIL-F-8785C  (9:22)  are: 

Kin  id  Min  tdwd  Min  wd 

0.4  0.4  1.0 

When  wd^ I 0/0 >d  is  greater  than  20  (rad/sec) 2,  the  ainiaua  ldwd  should  be 
increased  by  0.014  (wd^ 1 0/0 Id  “20)  .  These  are  ainiaua  values  only, 

however  the  data  base  does  not  support  increasing  the  ainiaua  daaping 
ratio  above  0.4  (10:494-502).  Also  the  dutch  roll  eigenvector,  and 
therefore  1 0/0 ld»  1*  as  yet  unspecified.  Therefore,  the  dutch  roll 
eigenvalue  aust  be  picked  so  any  anticipated  values  of  1 0/0  < d  will  still 
yield  Level  1  flying  quelities.  With  this  in  aind,  the  value  selected 
for  the  dutch  roll  eigenvalue  is: 

Kd  =  -1.2  i  2.751  (Jd  s  0.4;  wd  «  3.0)  (32) 

This  value  of  \d  seats  Level  1  criteria  for  10/01  less  than  8.8. 


V.  Eigenvector  Selection 


As  discussed  In  the  Section  11,  the  selection  of  eigenvectors  is 
difficult  because  of  the  T~1  distribution  of  control  inputs  and  initial 
conditions  asong  the  nodes.  However,  if  the  aircraft  ia  to  have  good 
handling  qualities  it  nust  respond  with  aircraft-like  dynanics  (10:545), 
regardless  of  the  conplexity  of  control  laws.  In  other  words,  the 
augnented  vehicle  should  have  doainant  nodes  corresponding  to  the  dutch 
roll,  roll  subsidence,  and  apiral.  The  selection  of  eigenvalues  in  the 
last  section  assured  the  augnented  nodes  will  be  near  the  frequencies  of 
the  usual  rigid  body  nodes.  But,  in  addition,  the  augnented  node  shapes 
nust  be  sinilar  to  those  of  the  classical  nodes.  For  exanple,  an 
aircraft  whose  roll  subsidence  response  was  prinarily  in  yaw  would 
likely  have  unacceptable  flying  qualities.  This  idea  leads  to 
restrictions  on  the  acceptable  eigenvectors.  Furthernore,  because  of 
the  fundanental  kinenatic  relationship  between  roll  rate  and  bank  angle, 
it  would  nake  no  sense  to  specify  a  response  in  one  of  these  eigenvector 
conponents  and  none  in  the  other.  While  these  intuitional  insights 
provide  useful  bounds  on  the  eigenvectors,  and  are  particularly  useful 
in  choosing  the  roll  subsidence  and  apiral  eigenvectors,  they  are  little 
help  in  selecting  the  dutch  roll  eigenvector. 

In  the  following  discussion,  eigenvectors  will  be  selected  to  yield 
Level  1  handling  qualities.  Selection  will  be  based  on  the  intuitional 
ideas  developed  earlier  as  well  as  the  aodal  characteristics  of 
ainplified  fourth  order  lateral-directional  models. 


The  eigenvector  coaponenta  were  developed  in  Section  111.  and  are 
repeated  here  for  convenience: 

v  =  C  r  0  p  0  6os  ir  V1  V2 
v3  v4  rwo  «f  1T 

Only  the  first  four  coaponenta  will  be  used  in  the  eigenvector  design. 
Therefore,  in  the  following  diacuaaion  they  are  the  only  ones  shown. 

The  other  coaponenta  are  left  unspecified. 

Roll  Subsidence  Eigenvector 

This  eigenvector  has  the  aost  effect  on  the  sideslip  which  develops 
while  rolling  and  is  vary  iaportant  in  attaining  good  turn  coordination. 
Sideslip  aay  develop  in  a  turn  either  froa  yaw-roll  coupling  or  through 
adverse  yaw  froa  the  aileron  input.  The  roll  subsidence  eigenvector 
will  deteraine  the  aagnitude  and  sign  of  the  yaw-roll  coupling. 

Ideally,  the  response  should  be  in  roll  rate  and  bank  angle  with  no 
sideslip  or  yaw  response.  That  is,  a  roll  rate  should  not  produce 
yawing  aoaenta  or  aideforce.  The  roll  subsidence  eigenvector  becoaes: 

vr8  *  t  0  0  1  x  IT  (33) 

where  the  x  indicates  the  value  is  not  specified  and  is  arbitrary.  It 
is  not  necessary  to  specify  the  bank  angle  aode  content  because  of  the 
kineaatic  relationship  between  p  and  phi.  Since  there  are  three 
coaponenta  specified,  they  will  not  be  achieved  precisely.  As  explained 
in  Section  II,  only  a  coaponenta  of  each  eigenvector  can  be  exactly 
achieved.  Since  there  are  only  two  independent  lateral-directional 


control  inputs  and  three  desired  coaponenta  of  the  eigenvector  to  be 


placed,  they  will  be  projected  into  the  achievable  elgenepace  In  a  least 
squares  Banner.  Therefore,  this  sods  will  still  sake  a  very  small 
contribution  to  the  sideslip  and  yaw  rate. 


An  example  Illustrates  the  effects  of  this  eigenvector.  It  is 
possible  to  use  this  eigenvector  to  control  the  sideslip  and  improve 
turn  coordination,  even  with  adverse  yaw  present.  For  the  classical 
fourth  order  lateral -directional  system,  yaw-roll  coupling  will  be 
minimized  when  Np  3  g/UQ.  However,  for  aircraft  with  large  N6os  the 
aileron  will  still  introduce  significant  yaw,  producing  an  uncoordinated 
roll.  The  roll  subsidence  eigenvector  can  be  selected  to  belance  the 
moment  produced  by  the  adverse  yaw  with  moments  from  the  roll  itself 
through  the  augmented  Np.  Consider  again  the  simple  example  for  the 
A-7D  used  in  section  II.  In  that  example,  the  roll  subsidence 
eigenvector  given  by  Eq  (33)  was  used  in  the  eigenstructure  assignment. 
If  the  roll  subsidence  eigenvector  is  instead  specified  by 

vr*  »  I  1  0  -12.6  x  ]T  (34) 

the  sideslip  response  to  a  step  aileron  input  becomes 

(S(t)  *  .14  e‘1.2t  ain(2.75t  ♦  108.6)  -.116  e'2.5t  -  .039  (35) 

Figures  1  and  2  show  the  beta  time  response  for  a  step  aileron  input  for 
the  two  cases.  Clearly,  the  second  eigenstructure  produces  much  less 
sideslip  end  would  likely  receive  a  higher  pilot  rating.  The  roll  rate 
response  is  unaffected  by  the  change  in  the  roll  subsidence  eigenvector 
and  is  shown  in  Figure  3  for  both  cases. 
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Figura  3.  Roll  Rate  vs  Tina  for  Step  Aileron  Input 


While  it  la  theoretically  possible  to  iaprove  turn  coordination 
with  this  method,  it  ia  not  practical.  An  aileron-rudder  croaafeed  as 
used  in  the  A-7D  design  and  many  other  modern  aircraft  is  a  better 
method  to  reduce  the  adverse  yaw.  Although  the  sideslip  buildup  for  a 
forced  response  is  reduced  using  the  eigenvector  in  Eq  (34),  the  free 
response  to  a  roll  rate  disturbance  would  be  unconventional,  with 
sideslip  content.  This  ia  undesirable  and  would  violate  the  notion  that 
an  aircraft  should  behave  like  a  conventional  aircraft  to  have  good 
flying  qualities.  Therefore,  Eq  (33)  will  be  used  to  specify  the 
ei genstructur e . 
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As  in  the  case  of  the  roll  subsidence  eigenvector*  an  understanding 
of  the  dynamics  aakes  specification  of  this  eigenvector  relatively  easy. 
Since  a  stable  spiral  aode  was  specified*  there  will  be  ssall  node 
contents  in  yaw  rate*  but  the  priaary  response  should  be  in  roll  rate 
and  bank  angle.  Sideslip  content  auat  be  zero  to  prevent  a  continuous 
buildup  over  this  long  acting  aode  and  is  necessary  for  a  coordinated 
turn.  The  desired  spiral  eigenvector  is: 

v*p  *  [  x  0  x  1  ]T  (36) 

Since  there  are  only  two  coaponents  specified*  they  will  be  achieved 
exactly  in  the  eigenstructure  assignaent. 

Dutch  Roll  Eigenvector 

The  dutch  roll  eigenvector  is  the  aost  difficult  to  specify. 

Because  the  eigenvalue  is  a  coaplex  quantity*  the  corresponding 
eigenvector  aust  also  be  coaplex.  Therefore,  the  designer  aust  specify 
at  least  soae  of  the  real  and  iaaginary  quantities.  Or  in  other  words, 
the  relative  aagnitude  and  phase  relationship  <or  portions  of  each)  of 
the  eigenvector  coaponents  aust  be  designated. 

Physical  intuition,  which  was  helpful  in  the  design  of  the  real 
eigenvectors,  gives  less  guidance  in  this  case.  So  the  fourth  order 
rigid  body  lateral -directional  equations  of  notions  discussed  in  Section 
III  will  be  used  to  develop  aodal  response  ratios.  By  exanining  the 
nodal  response  ratios  for  this  siaplifled  aodel,  it  will  be  possible  to 


gain  insight  into  the  aakeup  of  the  dutch  roll  eigenvector.  Since  one 
of  the  eigenvector  coaponents  is  arbitrary*  there  are  three  aodal 


raA, 


response  rstios  to  be  snalyzed.  The  ones  chosen  are  (p/0)d»  (r/0)d»  and 
<0/0>d. 

The  examination  of  the  modal  response  ratios  of  the  simplified 
equations  of  notion  will  begin  with  (p/0)d>  This  is  a  siaple  kinematic 
relationship  which  can  be  determined  exactly  (6:82).  The  magnitude  of 
the  response  ratio  is  determined  by  the  natural  frequency: 


I  p/0  Id  =  Wd 


Similarly,  the  phase  angle  is  given  by: 


arg(p/0)  *  */2  ♦  sin-1  td  rad 


Thus,  with  positive  damping,  the  roll  rate  will  always  lead  the  bank 
angle  by  at  least  90  degrees.  It  is  important  to  note  that  the  ratio  of 
these  eigenvector  components  is  determined  entirely  by  the  dutch  roll 
eigenvalue.  Once  a  desired  damping  ratio  and  natural  frequency  are 
chosen,  <p/0)d  is  uniquely  specified.  Also,  since  this  is  a  kinematic 
relationship,  it  will  hold  for  higher  order  linear  systems  as  well. 

The  ratio  of  the  yaw  rate  and  beta  eigenvector  components,  (r/0)d» 
is  also  determined  by  a  relatively  simple  relationship  (6:418): 


(r/0)d  3  -  Ad  *  *v  ♦  g/U0  (0/fl)d 


The  designer  has  little  sore  flexibility  here  than  for  (p/0)d* 
This  is  because  without  direct  sideforce  control,  the  effective  Yv  can 
only  be  slightly  modified  by  the  closed  loop  controls.  Thus,  once 
(0/0)d  snd  the  dutch  roll  eigenvalue  are  chosen,  (r/0>d  is  determined. 


^.■wwwv. Tr.-y.'y.v.vjv.  v.-'v. *■ .  jw  y.v.v.v.- irrwjw^.  <r .t#w-’V'.:v'. 


For  conventional  aircraft  with  aodarata  <0/0>< j,  tha  following 
approxiaationa  give  good  raaulta  avan  with  highar  ordar  dynaaica: 


lr/0ld  *  w<i 


arg(r/0>  ■  coa"l  <id>  "  *  rad 


Intuitively,  this  aakaa  aanaa.  If  aidaforcaa  ara  ignorad  and  tha  dutch 
roll  ia  “anakay"  with  aaall  <0/0>d»  aidaalip  rata  ia  equal  to  yaw  rata. 

It  ia  claar  that  for  convantional  aircraft  (p/0)d  and  (r/0>d  ara 
largaly  datarainad  by  tha  dutch  roll  aiganvalua.  Thua,  tha  laat  aodal 
reaponaa  ratio,  <0/0>dr  will  ba  aoat  critical  in  dataraining  tha  daairad 
atructura  of  tha  dutch  roll  aiganvactor.  For  a  fourth  ordar 
lateral -directional  aystea  (6:416-417): 


<0/0>d 


-Lrla  -  tYy  ♦  (L0/Lr) ] ) 
a2  -  Lpa  -  (g/U0)Lr 


where  a  ia  evaluated  at  tha  coaplex  dutch  roll  root.  Tha  approxiaata 
aagnituda  reaponaa  ia 


1 0/0 Id  • 


(L0*YyLp)2  ♦  2\  dwpLr (L0+YvLr)  *  ( ) 2l  1/2 


w])2(Lp 2  ♦  2ld*Dl*P  *  vi)2) 


and  tha  correaponding  phaaa  angle  ia  approxiaataly 


arg<0/e>d 


(l-?d2)^2<l*pl*0  ♦  2{dWDL0  ♦  w[>2Lr) 

tan“l  - - -  (44) 

»D<2Td2“l>L*  *  tdWD2l*r  ♦  ldwDM?L0  ♦  wpLrLp 

whera  wd  la  tha  daapad  fraquancy  of  tha  dutch  roll. 

Thus,  tha  problea  of  salacting  tha  dutch  roll  aiganvactor  has  baan 
raducad  to  salacting  an  Moptinua“  value  of  (0/0)d»  which  Is  typically  a 
coaplax  function  of  Lfi,  Lr,  and  tha  systaa  eigenvalues,  Xd  and  Xra* 
(0/0>d  1®  addressed  In  MIL-F-8785C  only  indirectly.  First,  as  1 0/0 Id 
increases,  tha  ainiaua  acceptable  dutch  roll  daaping  and  fraquancy  are 
increased  as  discussed  in  Section  IV.  As  long  as  tha  dutch  roll  has 
reasonable  fraquancy  and  daaping  tha  aagnituda  of  <0/0>d  can  be  quite 
large.  For  tha  eigenvalue  selected  in  this  report,  1 0/0 Id  a®  large  as 
8.6  still  neats  Laval  1  criteria.  Presently  no  lower  Halts  on  10/0 Id 
exist  because  of  conflicting  data.  Sons  data  (10:417)  indicate  pilots 
like  dihedral  to  daap  tha  dutch  roll  by  ailerons  alone.  This  is 
typically  done  by  KC-13S  pilots  in  the  power  approach  configuration.  On 
the  other  hand,  pilots  evaluating  the  YC-15  with  zero  effective  dihedral 
liked  it,  consenting  on  the  uncoupled  yaw  response  to  rudder  (10:417). 
With  the  dutch  roll  well  controlled,  pilot  daaping  should  not  be 
necessary,  so  zero  dihedral  should  not  be  a  problea.  Thus,  1 0/0 Id  froa 
zero  to  approxiaately  eight  neat  Level  1  criteria. 

Paragraph  3. 3. 6. 3  of  HIL-F-8785C  (9:33)  restricts  the  phase  angle, 
arg(0/0).  Pilots  have  been  trained  to  expect  the  aircraft  to  roll  left 
when  left  rudder  is  applied.  This  requireaent  specifies  that  for  an 
increase  in  right  sideslip,  the  pilot  aust  add  zero  or  right  bank  to 
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aaintaln  equillbriua.  Thia  la  aqui valent  to  specifying  the  effective 
dihedral  be  poaitive.  Dihedral  in  flying  qualitiea  work  noraally  refera 
to  the  phaaing  of  the  roll  and  aidealip  notion  in  the  dutch  roll, 
arg(p/0)<j.  With  poaitive  dihedral,  left  rolling  acconpaniea  right 
aidealipping  and  vice  veraa  in  the  dutch  roll  oacillation.  Thia 
definition  ia  still  aoaewhat  anbiguoua.  A  nore  preciae  definition  which 
aakea  intuitive  aenae  fron  a  pilota  point  of  view  ia  to  define  poaitive 
dihedral  aa  p  leading  0  by  90  deg  to  270  degreea.  At  90  and  270  degrees 
a  pilota  rudder  input  will  cause  no  rolling  notion  fron  the  dutch  roll, 
while  at  180  degrees  the  rolling  nonent  will  be  greatest.  Thia 
definition  also  appears  to  conply  with  paragraph  3. 3. 6. 3  of  MIL-F-8785C 
since  a  left  rudder  will  produce  a  left  rolling  nonent  fron  the  dutch 
roll.  However,  it  nay  not  always  produce  an  overall  aircraft  left 
rolling  nonent  due  to  rolling  nonent  contribution  fron  the  rudder 
deflection  itself  which  opposes  the  dutch  roll. 

The  phase  angle  between  p  and  phi,  113.6  degrees,  can  be  calculated 
directly  fron  Eq  (38)  using  the  design  daaping  ratio.  Uaing  thia 
relationship,  it  is  easy  to  show  the  range  of  arg(0/0)d  for  positive 
dihedral  aa  defined  above  is  fron  336.4  to  156.4  degreea.  Because  of 
the  rudder  contribution  to  the  rolling  nonent,  conpliance  with  paragraph 
3. 3. 6. 3  requires  thia  range  of  phase  angles  to  be  decreased.  For  thia 
study,  arg(0/0>d  fron  0  to  120  degreea  was  arbitrarily  chosen. 

The  last  requirenent  of  H1L-F-8785C  puts  iaplicit  Units  on  (0/0>d 
by  restricting  the  oscillations  in  the  roll  response  to  a  step  aileron 
input  (9:24).  When  the  dutch  roll  content  in  the  roll  response  becoaea 
excessive,  preciae  roll  response  becoaea  difficult  and  pilot  opinion 
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will  be  downgraded.  Paragraph  3. 3. 2. 2  restricts  the  asount  of 


oscillation  in  the  roll  response.  This  indirectly  limits  the  position 
of  the  zeros  of  the  quadratic  numerator  of  the  phi  to  4aa  transfer 
function,  which  in  turn  determines  the  degree  of  dutch  roll  excitation 
due  to  aileron  deflection.  The  indirect  effects  of  (0/0)d  can  be  seen 
in  the  following  equations  (6:463): 

(w0/wd)2  =  l  -  (L0/N0)  <N6aa/Uas>  (45) 

2(f0w0-?dwd)  *  L0/N0(Np  -  g/UQ)  (46) 


For  nonzero  N£as»  the  position  of  the  zero  relative  to  the  pole 
changes  as  L0  and  Lr  vary.  Changes  in  L0  and  Lr  will,  of  course  affect 
<0/0)d«  For  N6as  equal  to  zero,  the  dutch  roll  pole  and  zero  will 
exactly  cancel  and  (0/0)d  will  have  no  effect  on  the  phi  to  £8  transfer 
function.  Intuitively,  if  there  is  no  adverse  or  proverse  yaw,  the 
dutch  roll  node  will  not  be  excited  by  aileron  inputs  and  there  will  be 
no  oscillations  in  the  roll  response  regardless  of  (0/0)d.  When  the 
numerators  do  not  cancel,  the  dutch  roll  content  occurs  primarily  in  yaw 
and  sideslip  if  1 0/0 1 d  is  low  (less  than  1.5)  or  primarily  in  the  roll 
axis  when  1 0/0 1 d  is  large. 

This  requirement  of  MIL-F-8785C  generally  applies  for  aircraft  with 
dutch  roll  frequency  and  damping  close  to  the  Level  1  boundary.  In  this 
case,  the  dutch  roll  frequency  and  damping  were  chosen  to  be  well  within 
Level  1  criteria,  even  for  large  1 0/0 1 d .  Thus,  the  dutch  roll  effects 
on  the  roll  response  will  be  small  because  of  the  high  damping.  Also, 
adverse  yaw  in  the  YA-7D  is  effectively  controlled  using  aileron-rudder 
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feedback.  Thus,  this  spscif ication  adds  no  direct  restrictions  on 
<#/0>d- 

Another  effect  for  nonzero  Nias  is  the  reduction  of  steady  state 
roll  performance  with  increasing  L0/N0.  This  follows  from  Eq  (45)  and 
(6:376) 

• 

Pas  *  Liaa^Lp  (w0/wd)2  *aa  (47) 

For  adverse  yaw  (N£as  <  and  negative  L0,  significant  degradation 
in  roll  performance  can  occur.  The  effect  of  (0/0)d  is  not  obvious 
because  of  the  phase  angle,  but  the  reduction  in  roll  performance  should 
be  greatest  when  p  leads  0  by  180  degrees. 

HIL-F-8785C  puts  no  other  restrictions  on  (0/0)d»  So,  compliance 
requires  10/0 Id  lsss  than  8.6  and  arg(0/0)d  from  0  to  120  degrees. 
Another  consideration  is  response  in  gusty  wind  conditions.  For 
increasing  (0/0)d»  more  of  the  dutch  roll  response  will  occur  in  roll. 
Therefore,  in  gusty  wind  conditions  the  pilot  will  have  to  work  harder 
to  maintain  a  desired  bank  angle.  Exact  limits  on  pilot  acceptance  of 
this  condition  are  difficult  to  Identify,  more  so  in  this  case  since  the 
dutch  roll  is  very  well  damped  which  reduces  the  gust  perturbations. 
Thus,  there  is  no  easy  way  to  include  these  considerations. 

One  other  criteria  for  selecting  the  dutch  roll  eigenvector  has  no 
direct  bearing  on  the  handling  qualities.  Zt  has  been  shown  that  the 
eigenstructure  will  be  most  tolerant  to  errors  in  modeling  when  the 
eigenvectors  are  orthogonal  in  the  atate  space  (11).  Selecting  (0/0)d  of 


zero  will  meet  this  criteria. 


VI.  Flight  Control  Svatea  Design 

YA-7D  Digital  Hultlaode  Flight  Control  Svatea 

The  Digital  Multiaode  Flight  Control  Systaa  (DFCS)  was  designed  to 
duplicate  the  operation  and  functions  of  the  standard  Automatic  Flight 
Control  Systee  (AFCS)  and  to  provide  pilot  selectable  "aultiaodes"  which 
ispleaented  specified  control  laws.  For  this  test,  the  standard  DIGITAC 
eultieodes  were  replaced  with  the  controls  derived  in  the  this  report. 
The  new  control  laws  were  individually  selectable  froa  the  cockpit  using 
the  Hultieode  Control  Panel  shown  in  Figure  4.  The  control  laws  were 
selected  using  switches  P2,  P3,  and  P4.  The  P7  switch  was  used  to  sake 
the  roll  CAS  prefilter  tine  constant  cockpit  selectable  froa  0.300, 
0.125,  or  0.025  seconds.  The  reaaining  switches  were  not  used.  The 
standard  control  nodes  provided  by  the  DFCS  are  described  in  the  Flight 
Manual.  USAF  Series  A-7D  Aircraft  <12>  and  included  the  standard  CAS 
nodes  of: 

1.  Yaw  Stability  Augaentation 

2.  Pitch  and  Roll  Control  Augaentation 

3.  Attitude  Hold 

4.  Heading  Hold 

5.  Heading  Select 

6.  Altitude  Hold 

During  the  teat,  the  pitch  CAS,  and  roll  CAS  described  in  the  following 
paragraph  were  on.  The  other  CAS  nodes  were  used  between  teet  points. 


but  were  off  during  tasting. 


Figure  4.  Multiaode  Control  Panel 


YA-7D  Roll  Coaaand  Augaentation  Syatea  Design 


"ha  roll  and  yaw  control  ayateas  are  shown  in  Figures  5  and  6, 
respectively.  The  YA-7D  used  a  coabination  of  electrical  and  aechanical 
systeas  to  feed  pilot  lateral  inputs  to  the  ailerons.  The  aechanical 
systea  was  a  conventional  direct  link  to  the  aileron  actuator,  and  was 
Halted  to  15  degrees  of  coaaanded  deflection.  The  electrical  CAS, 
liaited  to  10  degreea,  provided  a  aore  unifora  roll  response.  The  CAS 
used  lateral  stick  force  as  the  input.  Because  the  eaphasis  of  this 
investigation  was  on  eigenstructure  assignaent  and  not  on  control 
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calculated  to  yield  eiganatructurea  for  a  range  of  (#/0)<i  which  aet 
requireaents  of  MIL-F-8785C  (9) .  Flight  taating  waa  than  conducted  on 
all  candidate  eiganatructurea.  The  teat  matrix  ia  ahown  in  Table  V. 

Gain  aatrlcea  were  calculated  uaing  the  computer  program  in  Appendix  A 
and  the  linear  aodela  for  the  two  flight  conditiona,  0.4  and  0.8  Mach  at 
15,000  feet  He.  The  control  law  gaina  are  ahown  in  Tablea  VI  and  VII 
for  0.4  and  0.6  Mach,  respectively. 

Table  V 


Phi  To  Beta  Teat  Matrix 

Phi /Beta 
Magnitude  Ratio 

Phi/Beta 

0 

Phaae  Angle 

60 

120 

0 

X 

1.5 

X 

X 

X 

3.0 

X 

X 

X 

.V.\\V  _\r 


Control  Law 

Galna  (0. 

4  MACH) 

PHI/BETA 

RATIO 

FEEDBACK 

GAINS 

MAGNITUDE 

PHASE  ANGLE  XI 

K2 

K3 

K4 

K5 

K6 

K7 

K8 

K9 

0 

0  .4596 

.0401 

.1088 

.1049 

1.172 

.0798 

-.0835 

-.0063 

.2000 

1.5 

0  .4115 

.0122 

.1094 

.0124 

1.178 

.0708 

-.0868 

-.0072 

.2000 

1.5 

60  -.0567 

.0002 

.1841 

.0133 

.9958 

.0650 

-.0581 

-.0070 

.2000 

1.5 

120  -.0734 

.0245 

.1932 

.0157 

.9670 

.0725 

-.0514 

-.0061 

.2000 

3.0 

0  .3381 

-.0144 

.1146 

.0102 

1.172 

.0621 

-.0882 

-.0080 

.2000 

3.0 

60  -.5030 

-.0353 

.2489 

.0119 

.8448 

.0520 

-.0365 

-.0076 

.2000 

3.0 

120  -.6117 

.0054 

.2786 

.0163 

.7609 

.0642 

-.0193 

-.0060 

.2000 

Not#:  Ka  aqual  0.02  for  all  control  lawa. 

n 

n 


L*j 


Tabla  VII 


Control  Law 

Gains  (0. 

.6  MACH) 

PHI/BETA 

RATIO 

FEEDBACK 

GAINS 

IGNITUDE 

PHASE  ANGLE  K1 

K2 

K3 

K4 

K5 

K6 

K7 

K8 

K9 

0 

0 

.3091 

.0155 

.0485 

.0057 

.5725 

.0160 

-.0033 

-.0032 

.0000 

1.5 

0 

.2897 

.0086 

.0428 

.0049 

.6098 

.0184 

-.0048 

-.0031 

.0000 

1.5 

60 

.0301 

.0045 

.0516 

.0046 

.5593 

.0167 

-.0040 

-.0033 

.0000 

1.5 

120 

.0336 

.0109 

.0575 

.0052 

.5238 

.0164 

-.0026 

-.0033 

.0000 

3.0 

0 

.2511 

.0016 

.0381 

.0042 

.6472 

.0187 

-.0063 

-.0031 

.0000 

3.0 

60 

-.2293 

-.0051 

.0549 

.0036 

.5492 

.0155 

-.0046 

-.0034 

.0000 

3.0 

120 

-.2325 

.0066 

.0664 

.0048 

.4814 

.0148 

-.0019 

-.0035 

.0000 

Not*:  K*  *qual  0.02  for  all  control  laws. 
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Handling  qualities  are  those  qualities  or  characteristics  of  an 
aircraft  that  govern  the  ease  and  precision  with  which  a  pilot  is  able 
to  perform  tasks  required  in  support  of  an  aircraft  role  (13).  This 
definition  has  several  iaplicationa  which  are  important  to  understand. 

It  requires  the  pilot  in  the  control  loop,  therefore  evaluation  of  the 
open  loop  response  alone  is  insufficient  to  insure  acceptable  handling 
qualities.  <A  good  example  of  the  importance  of  analysis  with  pilot  in 
the  loop  is  the  dramatic  effect  tine  delays  in  the  forward  loop  can  have 
on  pilot  opinion.)  Therefore,  the  eigenatructure  of  the  lateral- 
directional  axes,  which  determine  stability,  is  only  one  portion  of  the 
total  control  system  which  affects  the  handling  qualities.  It  is  also 
important  to  realize  the  handling  qualities  must  be  evaluated  in  a  task 
required  in  the  aircraft  role.  For  fighter  aircraft  such  as  the  YA-7D, 
precision  in  air-to-air  or  air-to-ground  tracking  is  a  valid  criteria. 

Two  criteria  will  be  used  to  evaluate  the  handling  qualities 
analytically,  bank  angle  tracking  and  turn  coordination.  These  are 
typically  the  parameters  of  most  interest  in  lateral -directional  notion. 
Precision  in  bank  angle  control  is  necesaary  for  the  pilot  to  be  able  to 
precisely  orient  the  lift  vector  and  turn  the  aircraft.  In  addition  to 
good  bank  angle  control,  satisfactory  turn  coordination  is  necessary  for 
precision  in  nose  pointing.  Bank  angle  tracking  can  be  thought  of  as  an 
inner  loop  closure  with  heading  control  as  the  outer  loop.  If  the  inner 
loop  bank  angle  characteristics  and  turn  coordination  are  good,  the 
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pilot  will  ba  able  to  precisaly  control  both  lift  vector  and  nose 
poaition . 

In  the  following  analysis,  the  augaented  eiganatructurea  at  0.4  and 
0.6  Mach  will  be  used.  The  reaulting  12th  order  eigenatructures  are 
still  insufficient  to  completely  analyze  the  handling  qualities.  The 
two  linear  prefilters  in  the  CAS  are  of  low  enough  frequency  that  they 
aust  be  considered.  The  aechanical  system,  which  coaaanda  the  aileron 
power  actuator  directly,  has  a  gain  of  .01453  rad/lb.  In  state  space 
fora,  the  dynaaics  of  the  two  prefilters  and  aechanical  systaa  aay  be 
described  as: 


x  *  A  x  ♦  B  u 


x  *  l  FI  4 esc  <aa  *2  1T 


=  I  Fs  4rc  1T 


3  0 

0  0 
.291  0 

.  0  0 


where  FI  is  the  first  stick  state  and  Fs  is  the  stick  force.  The  other 
states  and  inputs  are  as  previously  defined,  and  only  the  changed 
portion  of  the  Aas  «nd  V2  states  has  been  shown.  The  14th  order 
eigenatructures  will  be  the  basis  for  analysis. 


Estimation  Techniques 

Handling  qualities  estimation  is  still  a  very  inexact  procedure  as 
evidenced  by  the  amount  of  data  in  MIL-F-8785C  (9>  which  does  not 
correlate  well  with  the  specif icatlona  and  the  amount  of  research  being 
conducted  on  handling  qualities  of  highly  augaented  aircraft. 


Nevertheless,  most  aircraft  which  aaat  the  specif i cations  will  have  good 
flying  qualities.  This  points  to  the  difficulty  in  analyzing  the  flying 
qualities  of  the  eigenstructures  developed  in  the  previous  sections. 

Even  after  eliminating  the  eigenstructures  which  do  not  meet  the 
criteria  of  MIL-F-8785C.  there  are  nany  left  that  do.  Differences  in 
the  handling  qualities  between  then  can  be  expected  to  be  snail. 

Nany  techniques  have  been  developed  for  analyzing  pitch  dynanics  of 
higher  order  systess,  including  pilot  nodels  such  as  the  Neal-Smith 
criteria  (10:228),  and  several  criteria  on  the  tlse  response  (10:233). 
The  tine  response  techniques  do  not  have  direct  application  to 
lateral -directional  dynamics,  and  pilot  models  are  somewhat  dependent  on 
the  criteria  assumed  for  the  pilot.  Hoh's  bandwidth  criteria  (14  and 
10:213)  is  a  technique  that  is  intuitively  attractive  and  can  be  adapted 
to  a  wide  variety  of  problems. 

One  measure  of  the  handling  qualities  of  an  aircraft  is  its 
stability  margin  when  operated  in  a  closed  loop  tracking  task.  In 
classical  control  theory,  the  maximum  frequency  at  which  such  closed 
loop  tracking  can  occur  without  threatening  stability  is  the  bandwidth, 
wbw*  The  pilot  can  be  thought  of  as  a  pure  gain  controller  with  unity 
feedback  as  in  Figure  8.  Although  a  simple  assumption,  an 


aircraft  with  good  handling  qualitiea  will  require  minimal  pilot 
compenaation .  It  ia  further  asauaed  the  pilot  will  close  the  loop  and 
add  gain  to  achieve  45  degrees  of  phase  margin  or  6  decibels  (dB)  of 
gain  margin,  whichever  produces  the  lower  bandwidth.  The  higher  the 
bandwidth,  the  more  precise  the  pilot  will  be  able  to  control  the  ayatea 
output  without  going  unstable.  Bandwidth  will  be  determined  from  the 
Bode  plot  using  the  criteria  in  Figure  9.  This  method  makes  the  system 
bandwidth  independent  of  the  actual  pilot  gain. 

Another  aspect  of  handling  qualities  is  the  ability  of  the  pilot  to 
increase  the  bandwidth  by  increasing  gain  while  adding  compensation. 


Figure  9.  Bandwidth  Defined  by  Gain  and  Phase  Margins 
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The  capability  to  do  thia  ia  tied  to  the  phaae  characteristics  above  the 
bandwidth  frequency.  Syateaa  which  have  a  rapid  falloff  in  phaae  can 
tolerate  only  aaall  gain  increases  and  require  large  pilot  coapensation. 
The  high  frequency  phase  characteristics,  which  are  a  function  of  the 
higher  order  dynamics,  can  be  represented  by  an  equivalent  amount  of 
pure  time  delay,  r«.  The  equivalent  time  delay  can  be  accurately 
estimated  from  the  frequency  response  using  the  equation: 

.  -  #1  -  180 

«  -  (49) 

57.3  wi 

where  01  and  wi  are  the  phase  angle  and  frequency  at  twice  the  frequency 
where  the  phase  angle  is  180  degrees,  widO  <10:218). 

The  third  criteria  for  satisfactory  handling  qualities  is  for  the 
magnitude  response  to  have  a  K/s  slope  in  the  region  of  the  bendwidth 
where  the  pilot  is  operating  at  maximum  gain.  A  "shelf"  in  the 
magnitude  response  in  this  region  ia  highly  undesirable,  indicating  a 
rapid  decrease  in  phase  margin  as  gain  is  increased. 

Hence,  the  f lying  qualities  in  a  particular  closed  loop  task  are 
dependent  on  three  parameters: 

1)  The  bandwidth,  wpw 

2)  The  shape  of  the  phase  curve  at  frequencies  above  wbw,  as 
measured  by  the  effective  time  delay,  r«. 

3)  The  shape  of  the  magnitude  curve  near  wbw 

The  bandwidth  required  to  achieve  satisfactory  (Level  1)  handling 
qualities  is  dependent  on  the  task.  In  a  bank  angle  tracking  task,  the 
pilot  is  closing  the  phi  to  stick  force  loop,  therefore  the  transfer 


function  of  interest  is  #/Fs.  Exact  categorization  of  the 
aigenstructuras  into  Laval  1,  Laval  2,  etc.  ia  difficult.  By  analogy  to 
a  pura  first  roll  order  response,  in  which  the  bandwidth  is  equal  to  the 
roll  subsidence  eigenvector,  it  would  be  reasonable  to  asauae  bandwidths 
between  one  and  three  radians/aecond  will  produce  Level  1  handling 
qualities.  Bandwidths  less  than  one  radian/second  will  be  too  sluggish, 
and  greater  than  three  radians/ second  too  abrupt.  Also,  Reference  10 
(10:407)  indicates  a  degradation  in  handling  qualities  when  re  exceeds 
100  nsec.  This  will  be  explained  in  aore  detail  later. 

Techniques  to  analyse  turn  coordination  are  not  as  wall  developed. 
The  pilot  loop  closures  are  not  as  clear  and  the  notion  involves  two 
axes.  One  Measure  of  turn  coordination  is  the  pilot  rudder  input 
necessary  to  Maintain  sideslip  at  zero.  Large  or  complex  rudder  inputs 
are  indicative  of  complex  nose  notion  which  will  increase  pilot  workload 
during  tracking  maneuvers.  Reference  10  (10:542)  presents  guidelines 
for  allowable  rudder  shaping,  however,  they  are  based  on  historical  data 
for  coordination  in  the  landing  configuration.  Rather  than  evaluate 
coapliance  with  these  guidelines,  a  slightly  different  approach  will  be 
taken  based  on  siailar  principles.  While  specific  categorization  of 
flying  qualities  will  not  be  possible,  a  better  understanding  of  the 
effects  of  the  various  aigenstructuras  is  possible. 

The  rudder  crossfeed,  Ycf,  necessary  to  Maintain  zero  sideslip  in  a 
turn  can  be  defined  as: 

<rP  8  Ycf  Fa  (50) 

For  a  particular  aircraft,  YCf  nay  be  evaluated  as: 
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Vcf  » 


Numerator  of  0  to  Fa  tranafar  function 


(51) 


Nuaerator  of  0  to  6rP  tranafar  function 

Ycf  can  be  calculated  froa  the  the  14th  order  eigenstructures,  and 
the  frequency  response  plotted.  An  understanding  of  the  dynaaica  aids 
in  interpretation.  Sideslip  in  a  turn  can  arise  froa  two  sources, 
adverse/ prover as  yaw  due  to  aileron  deflection,  and  yaw  coupling  due  to 
roll.  Pilots  typically  can  tolerate  sore  adverse  yaw  than  proverae, 
since  the  rudder  deflection  is  in  the  sane  sense  as  the  stick  input  and 
is  sore  natural.  The  eigenstructure  will  only  effect  the 
adverse/ proverae  yaw  indirectly.  The  aileron  to  rudder  feedback  ia 
designed  to  ainiaize  the  adverse  yaw,  and  the  higher  dutch  roll 
frequency  resulting  froa  augaentation  will  help  ainiaize  the  aagnitude 
of  any  resulting  sideslip.  Low  values  of  <0/0>d  will  result  in  sore 
sideslip,  since  aore  of  the  dutch  roll  response,  excited  by 
adverae/proverse  yaw,  ia  in  beta.  The  second  contribution  to  sideslip 
is  froa  roll-yaw  coupling,  which  was  ainiaized  through  the  design  of  the 
roll  subsidence  eigenvector. 

Froa  a  handling  qualities  viewpoint,  large  aagnitudea  of  the 
crossfeed  paraaeter,  IYCfl,  in  the  S  to  10  rad/sec  frequency  range 
indicates  the  pilot  aust  initially  apply  a  large  rudder  input  to  keep 
sideslip  at  zero.  Phase  angles  of  the  crosafeed  paraaeter,  arg<YCf), 
greater  than  180  degrees  in  this  frequency  range  indicates  proverse  yaw, 
which  aay  be  objectionable  with  large  !Ycf*-  If  arg(Ycf>  goes  through 
180  degrees  <in  either  direction)  in  the  frequency  range  one  to  five 
rad/sec,  the  pilot  will  have  to  reverse  pedal  deflection  during  the  roll 
which  aay  also  be  objectionable.  Frequency  ranges  are  approxiaate,  and 
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merely  indicate  the  ranges  in  which  pilots  aoat  often  operate  in  high 
gain  tasks. 

The  techniques  chosen  here  are  not  unique.  Other  techniques  are 
available  which  provide  similar  estimation  of  handling  qualities.  The 
main  advantage  of  working  in  the  frequency  domain  is  the  effects  of 
higher  order  dynamics  can  be  seen.  Zt  must  be  remembered  the  gains  only 
control  the  nodal  response  of  the  four  basic  rigid  body  modes.  Other 
modes,  such  as  the  rudder  and  aileron  dynamics  may  move  into  a  frequency 
range  noticeable  to  the  pilot.  The  roll  CAS  dynamics  will  certainly 
have  a  significant  effect  on  the  aircraft  roll  response.  Also,  only 
very  simple  assumptions  about  pilot  behavior  are  required. 


An  example  of  bank  angle  tracking  using  classical  dynamics  will 
illustrate  the  technique.  For  classical  airplanes,  roll  bandwidth  is 
only  a  function  of  Ara,  and  w<j.  When  there  is  no  roll -yaw  coupling 
<w0  a  wd),  the  bandwidth  will  be  phase  margin  limited  and  will  be 
exactly  equal  to  as  shown  in  Figure  10.  As  is  well  known,  when  wp 
is  greater  than  wd,  closing  the  loop  reduces  system  stability.  This  is 
also  reflected  in  the  lower  bandwidth  shown  in  Figure  11. 

It  has  already  been  shown  how  w0/wd  is  influenced  by  the  dutch  roll 
eigenvector,  thus  it  is  expected  variations  in  (0/0>d  will  influence 
roll  stability.  Of  course  the  damping  and  frequency  of  the  dutch  roll 
will  also  influence  this  loop  closure,  however  the  augmentation  should 
have  increased  these  to  satisfactory  values. 
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Frequency  response  plots  for  the  unaugaented  YA-7D  are  shown  in 
Appendix  C  for  the  unaugaented  YA-7D  and  the  seven  lateral-directional 
eigenstructures  described  in  the  previous  section.  Results  are 
sunaarized  in  Table  VIII. 

At  0.4  Mach  the  unaugaented  YA-7D  has  acceptable  bandwidth  and  tine 
delay  characteristics,  but  has  a  reversal  in  the  gain  and  phase  curves 
near  the  bandwidth  frequency  due  to  a  poorly  daaped  dutch  roll  with 
aoderate  adverse  yaw  (Figure  Cl).  The  pilot  aay  have  difficulty  in  high 
gain  tasks  because  of  the  rapid  decrease  in  phase  aargin  as  gain  ia 
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Figure  11.  Effect  of  Roll-Yaw  Coupling  on  Roll  Dynaaica 


increaaed  in  thia  frequency  range.  Figurea  C2  to  C8  ahow  the  reaponae 
ia  auch  iaproved  with  augaentation.  Thera  ia  no  longer  a  ahelf  in  the 
aagnitude  reaponae,  and  the  effective  tiae  delay  ia  acceptable.  Even 
though  the  roll  aode  waa  daaigned  to  give  a  bandwidth  of  four  radiana 
per  aecond,  the  prefiltera  in  the  atick  and  other  higher  order  effecta 
have  alowed  the  reaponae  conaiderably .  With  the  aileron  to  rudder 
feedback,  adverae  yaw  ia  well  controlled,  although  not  perfectly. 

Little  variation  in  the  reaponae  ia  noted  with  l0/0ld  of  1.5.  With 
10/.  Id  of  three,  arg(0/0>d  now  haa  a  noticeable  effect  on  both  bandwidth 


58 


Table  VIII 


Bandwidth  and  Tiae  Delay  of  Eigenstructurea 

Mach 

<0/0>d 

Bandwidth 

Tiae  Delay 

0.4 

unaugaented 

2.17 

0.066 

0 

1.43 

0.077 

l.S  arg(0) 

1.51 

0.093 

1.3  arg(60) 

1.50 

0.081 

l.S  arg(120) 

1.43 

0.067 

3.0  arg<0) 

2.10 

0.105 

3.0  arg(60) 

2.93 

0.062 

3.0  arg(120) 

1.43 

0.066 

0.6 

unaugaented 

1.83 

0.062 

0 

1.57 

0.074 

l.S  arg(0) 

1.64 

0.081 

1.5  arg<60) 

1.61 

0.083 

1.5  arg<120> 

1.54 

0.074 

3.0  arg<0) 

1.72 

0.091 

3.0  arg<60) 

1.65 

0.096 

3.0  arg<120) 

1.31 

0.077 

and  time 

delay.  The  shape  of 

response  has  a  noticeable 

shelf  at 

arg(0/p)d  of  120  degrees.  It  is  sxpsctsd  that  handling  qualities  will 
be  satisfactory  with  all  eigenstructurea  with  the  possible  exception  of 
(0/ 0 )  (j  equal  to  3  arg(120) . 

As  shown  in  Figure  C9  in  Appendix  C,  the  roll  characteristics  of 


the  YA-7D  at  0.6  Mach  are  already  satisfactory.  There  is  only  a  slight 


amount  of  provaraa  yaw,  so  variations  in  (0/0>d  has  almost  no  offset. 


The  magnitude  responses  shown  in  Figures  CIO  to  C16,  have  the  desired 
K/s  shape,  and  bandwidth  and  time  delay  are  satisfactory.  All  the 
feedback  gains  should  produce  satisfactory  handling  qualities. 

Evaluation  of  Turn  Coordination 

Figures  C17  to  C32  show  frequency  response  plots  for  the  rudder 
crossfeed  parameter  at  0.4  and  0.6  Mach.  The  magnitude  and  phase 
responses  show  only  slight  variations  with  the  different 
eigenstructures . 

At  0.4  Mach,  the  unaugmented  YA-7D  shows  an  initial  rudder 
deflection  needed  to  counter  the  adverse  yaw,  followed  by  a  slight 
decrease,  then  steadily  increasing  rudder  as  the  turn  develops. 

Although  the  response  shows  typical  high  frequency  characteristics,  the 
low  frequency  response  shows  a  poorly  coordinated  turn.  All  augmented 
eigenstructures  show  similar  responses,  with  deflection  needed  to 
counter  the  adverse  yaw,  followed  by  decreasing  deflection  as  the  turn 
develops. 

At  0.6  Mach,  the  unaugmented  YA-7D  needs  only  slight  rudder 
initially,  but  again  has  poor  low  frequency  characteristics.  With 
augmentation  the  low  frequency  characteristics  are  such  improved. 

In  summary,  variations  in  <0/0)d  have  only  slight  effects  on  turn 
coordination.  Whan  adverse/proverse  yaw  is  small,  variations  in  (0/0)d 
also  have  small  effects  on  the  roll  response.  When  adverse  yaw  is 
present,  increasing  1 0/0 1 q  can  have  an  adverse  impact  on  the  shape  of 
the  magnitude  response,  particularly  for  phase  angles  near  120  degrees. 


VIII.  Simulator  Raaulta 

The  siaulator  for  aircraft  Flight  Teat  and  Oevelopaent  (SAFTD)  was 
a  low  fidelity  aiaulator  conaiating  of  a  fixed  generic  fighter  cockpit 
with  a  cathode  ray  tube  diaplay.  Prograaaable  aoftware  enabled 
aiaulation  of  a  variety  of  aircraft.  The  YA-7D  aiaulation  uaed  a  five 
degree  of  freedom  linear  aerodynaaic  model  with  a  digital  eaulation  of 
the  DIGITAC  control  ayatea.  Open  loop  taaka  described  in  Appendix  F 
were  flown  by  project  pilota  and  reaulta  recorded  on  two  eight  channel 
atripchart  recorders.  These  data  were  uaed  to  verify  the  control  laws 
produced  the  desired  response.  Cloaed  loop  taaka  were  not  evaluated  in 
the  aiaulator  due  to  excessive  tiae  delay  in  the  visual  diaplay  and  poor 
fidelity  of  the  aiaulation. 

Initial  SAFTD  reaulta  indicated  the  control  laws  perforaed  as 
predicted.  The  cloaed  loop  eigenstructure  aa  aeasured  by  the  roll  mode 
tiae  constant  <rr),  dutch  roll  frequency  (w<j>  and  damping  ratio  <Td>» 
and  (0/p)d»  agreed  well  with  design  parameters.  Dutch  roll  damping 
tended  to  be  slightly  low,  but  was  still  within  15  percent  of  the 
ainiaua  requirement  and  design  condition  (0.4)  listed  in  HIL-F-8785C 
(S) .  The  stick  force  command  gain  produced  acceptable  roll  rate  per 
pound  gradients  <pss/Fs).  Tables  VIII  and  IX  suaaarize  tha  results  at 
0.4  and  0.6  Mach.  The  lateral-directional  characteristics  of  the 
unaugaented  YA-7D  and  standard  aircraft  with  augmentation  are  included 
for  coaparison. 

Qualitatively,  control  systea  designs  at  both  0.4  and  0.6  Mach  flew 
well  with  no  adverse  characteristics.  The  roll  CAS  command  gain,  Ks, 
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was  evaluated,  and  a  constant  value  of  0.02  rad/sec  per  pound  selected 
for  all  control  laws. 

Representative  tiae  response  plots  for  a  step  lateral  stick  input 
and  release  from  wings  level  sideslip,  generated  in  the  SAFTD,  are 
included  in  Appendix  C. 


Table  IX 


Simulator 

Lateral 

-Directional 

Characteristics  -  0.4 

Mach 

Design  <0/0)d* 
(aag/phase  ang) 

Wd 

Id 

Tr 

Pss'Fa 

(0/0>d 

Unaugaented 

1.8 

0.16 

0.40 

5.6 

3.4  arg(46) 

DFCS  on2 

2.9 

0.70 

0.30 

6.7 

1.9  arg<72) 

0/0 

3.2 

0.35 

0.33 

6.0 

0  arg(0) 

1.5  /  0 

3.2 

0.35 

0.30 

6.4 

1.7  arg(l7) 

1.5/60 

3.2 

0.37 

0.28 

5.2 

1.8  arg(43) 

1.5  /  120 

3.1 

0.38 

0.25 

5.2 

1.4  arg(106) 

3.0  /  0 

3.6 

0.40 

0.28 

5.6 

3.0  arg(0) 

3.0  /  60 

2.9 

0.40 

0.25 

4.7 

3.2  arg(46) 

3.0  /  120 

2.9 

0.37 

0.25 

4.4 

3.6  arg(105) 

^  Note:  Design  conditions  were  rr  =  0.25,  w<j  =  3.0,  and  id  =  0.40 
2  Note:  Includes  effects  of  pref liters  in  CAS. 


Simulator  Lateral -Directional  Characteristic*  -  0 


Design  <#/0>d^ 
(sag/ phase  ang) 

Wd 

Id 

rr 

Paa'F» 

Unaugaented 

1.9 

0.18 

0.40 

6.3 

DFCS  on2 

2.9 

0.70 

0.40 

7.5 

0/0 

3.2 

0.33 

0.30 

10.3 

1.5  /  0 

3.0 

0.40 

0.30 

10.0 

1.5  /  60 

3.1 

0.38 

0.25 

10.0 

1.5  /  120 

3.1 

0.42 

0.30 

10.3 

3.0  /  0 

3.2 

0.38 

0.30 

10.0 

3.0  /  60 

3.3 

0.43 

0.33 

10.0 

3.0  /  120 

3.2 

0.45 

0.25 

10.0 

&  Mach 

<#/8>d 

3.3  arg(38) 
1.0  arg<120) 
0  arg(0) 
1.7  arg<24> 
2.0  arg(57) 

2.2  arg(112) 
3.1  arg(9) 

3.3  arg(51) 
3.1  arg(114) 


1  Note:  Design  conditions  were  rr  *  0.25,  w<j  *  3.0,  and  *  0.40 

2  Note:  Includes  effects  of  prefilters  in  CAS. 
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IX.  Flight  Test  Results 


General 

A  limited  flight  test  progras  baaed  on  the  control  laws  developed 
in  this  report  was  conducted  at  the  Air  Force  Flight  Test  Center 
(AFFTC),  Edwards  Air  Force  Base  (AFB),  California,  between  8  October  and 
17  November  1987.  Seventeen  sorties  (28.3  flight  hours)  were  flown  in 
the  YA-7D  DIGITAC  aircraft.  Seventeen  additional  chase  sorties  <29.7 
flight  hours)  were  flown  in  A-7D  aircraft.  All  tests  were  flown  within 
liaitationa  contained  in  the  Flight  Manual.  USAF  Series  A-7D  Aircraft 
(12)  and  the  Partial  Flight  Manual,  YA-7D  Serial  Humber  67-14583  (15). 

A  more  detailed  description  of  the  test  limitations  and  complete 
description  of  test  procedures  is  given  in  the  HAVE  EIGEN  Test  Plan 


(16). 

Obiectlves 

The  overall  objective  of  the  HAVE  EIGEN  test  program  was  to 
evaluate  the  lateral-directional  handling  qualities  of  the  YA-7D  at 
15,000  feet  pressure  altitude  (He),  0.4  and  0.6  Mach,  with  the  control 
laws  designed  from  eigenstructure  assignment.  Specific  objectives  were 
to: 


1.  Confirm  the  closed  loop  eigenstructure  of  the 
lateral -directional  axes  agreed  with  the  design 
eigenstructure . 

2.  Correlate  pilot  handling  quality  ratings  to  <0/0>c,' 

3.  Evaluate  the  use  of  eigenstructure  assignment  as  a  control 


law  synthesis  technique 
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Aircraft  Description  and  Instrumentation 

The  YA-7D  teat  aircraft,  USAF  serial  number  (S/N)  67-14563,  was  a 
prototype  for  the  A-7D  single  seat,  light  attack  aircraft  and  was 
powered  by  a  single  Allison  TF-41-A-1  engine.  A  coaplete  description  of 


A-7D  is  given  in  the 


(12).  The  test  aircraft  was  modified 


for  the  DIGITAC  test  program  as  described  in  the  Partial  Flight  Manual 
(15).  An  instrumentation  pod  on  the  aircraft's  right  inboard  wing 
station  provided  data  acquisition  and  telemetry  (TH)  capability.  An 
inert  Hk-82  bomb  was  installed  opposite  the  OZGITAC  instrumentation  pod 
for  symmetry.  The  two  outboard  stations  on  each  wing  had  pylons,  but 
carried  no  stores. 

The  test  aircraft  had  a  control  augmentation,  stick  steering  flight 
control  system.  For  the  DIGITAC  program,  the  analog  computer  of  the 
AFCS  was  replaced  with  a  programmable  DFCS.  The  mechanical  flight 
control  system  was  unchanged  by  the  installation  of  the  DFCS.  For  this 
test,  the  longitudinal  control  laws  and  lateral  CAS  prefilters  for  the 
standard  YA-7D  were  used.  The  lateral -directional  control  laws  derived 
from  eigenstructure  assignment  were  added  to  the  DFCS. 

The  weapon  sight  system  of  the  YA-7D  DIGITAC  aircraft  consisted  of 
a  gunsight  computer  electronic  unit,  an  inertial  sensor  unit,  a  control 
panel,  a  signal  generator,  and  a  standard  A-7D  Head-up  Display  (HUD) 
(AN/AVQ-7(V) ) .  All  controls  on  the  HUD  control  panel  functioned  as 


described  in  the 


(12)  with  the  exception  of  the  TEST, 


SCALES,  and  BARO/RDR  ALT  switches,  which  were  deactivated. 

The  HUD  was  equipped  with  a  camera  that  recorded  HUD  video  onto 


16am  color  film.  The  camera  was  armed  by  a  switch  on  the  the  pilot's 
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right  rear  consols.  The  fils  cssera  was  activated  by  depressing  the 
trigger  switch  to  the  first  detent,  or  by  activating  the  camera  run 
switch  on  the  left  canopy  bow. 

The  test  aircraft  was  also  equipped  with  a  Yaw,  Angle-of-attack, 
Pitot-Static  (YAPS)  noaa  boom,  Ampex  AR-700  magnetic  tape  recorder,  TH 
transmitter,  Millikan  DBM-2  Head-Up  Display  camera,  and  sensitive 
cockpit  instruments.  A  complete  description  of  the  data  acquisition 
system  was  in  the  Partial  Flight  Manual  (15).  Specific  parameter 
measurement  device  descriptions  and  resolutions  were  in  the 
Instrumentation  Handbook  (17).  Additional  information  on  test 
instrumentation  is  included  in  Appendix  E. 

The  weight  of  the  aircraft  varied  from  approximately  31,000  pounds 
at  takeoff  to  24,000  pounds  at  landing  with  2,000  pounds  of  fuel 
remaining.  The  center  of  gravity  <CG>  moved  aft  from  27.1  percent  mean 
aerodynamic  chord  (MAC)  to  27.6  percent  MAC  as  fuel  was  burned.  The 
center  of  gravity  also  dropped  from  2  inches  above  the  waterline  at 
engine  start,  to  the  waterline  at  the  end  of  the  mission.  The  weight 
changes  and  center  of  gravity  movement,  as  well  as  longitudinal  control 
system  characteristics,  had  negligible  effects  on  the  lateral- 
directional  flying  qualities  of  the  aircraft. 

Test  Methods  and  Conditions 

Testing  was  conducted  in  three  phases.  Ground  testing  of  the 
aircraft  was  conducted  at  twice  the  normal  loop  gain  to  check  for  limit 
cycle  and  structural  resonance  oscillations.  Simultaneously,  the  open 


loop  simulator  testing  described  in  the  previous  section  was  performed 


to  evaluate  how  closely  the  control  laws  achieved  the  desired  response. 
And  finally,  open  and  closed  loop  airborne  testing  was  flown. 

Open  and  closed  loop  testing  was  performed  with  each  control  law. 
Stripchart  data  were  used  to  estlsate  the  roll  aode  tlse  constant,  dutch 
roll  frequency  and  damping,  and  (0/9)4.  Further  information  on  open 
loop  testing  la  given  in  Appendix  F. 

Three  closed  loop  tests  were  performed;  two  HUD  tracking  tasks  and 
an  air-to-air  tracing  task  with  another  A-7D  aa  a  cooperative  target. 

The  project  pilots  rated  each  control  law  in  each  taak  at  least  twice. 
The  test  engineer  assigned  the  points,  and  to  prevent  biaaing,  the  , 
pilots  did  not  know  which  phi/beta  ratio  was  being  tested.  The  specific 
criteria  and  description  of  the  tasks  is  given  in  more  detail  in 
Appendix  F. 

Test  and  Evaluation 

Ground  Testing.  Initial  ground  testing  of  the  control  systems 
revealed  small  amplitude,  high  frequency  oscillations  in  the  ailerons 
with  several  of  the  control  laws.  Smaller  oscillations  were  also 
present  in  the  rudder.  Investigation  revealed  high  gains  in  the  lateral 
accelerometer  to  aileron  channel  were  amplifying  structural  vibrations 
producing  a  structural  resonance  in  the  flight  control  system  centered 
at  10.5  Hertz  (Hz).  Frequency  analysis  showed  actual  modes  of  the 
resonance  to  have  frequency  content  at  10.5,  13  and  15  Hz.  The  first 
lateral  bending  mode  of  the  YA-7D  was  at  10.5  Hz  and  the  natural 
frequency  of  the  lateral  accelerometers  was  at  13  Hz.  It  was  felt  these 
were  the  causes  of  the  structural  resonance.  Unfortunately,  the  rather 


broad  frequency  band  prevented  uaing  a  notch  filter  in  the  accelerometer 


channel . 


The  10  rad/aec  first  order  noise  filter,  described  in  Section  III, 


was  not  initially  included  as  part  of  the  design.  When  the  noise  filter 


was  added  to  the  ground  testing,  it  eliainated  the  feedback  at  noraal 


gain  levels.  Standsrd  AFFTC  procedure  was  to  ground  test  at  twice 


noraal  gain  to  provide  at  least  six  decibels  <Db>  of  gain  aargin  froa 


the  resonance.  At  these  increased  gain  levels  a  saall  amount  of 


structural  resonance  was  still  present.  Finally,  a  non-linear  digital 


filter  already  used  in  the  standard  DIGITAC  yaw  channel  was  added  to  the 


first  order  lag  filter.  The  non-linear  filter  had  very  little  phase  or 


gain  loss  below  five  Hz  and  therefore  had  little  effect  on  the  dynaaic 


aodea  of  interest.  Above  5  Hz,  the  gain  dropped  off  rapidly,  providing 


11  Ob  of  gain  margin  above  10  Hz.  This  was  sufficient  to  eliminate  all 


structural  resonances.  The  gains  of  the  control  lews  were  recomputed 


with  the  10  rad/sec  lag  filter  added  in  the  acceleroaeter  channel,  and 


the  configuration  reflown  in  the  siaulator  with  both  lag  filter  and 


non-linear  digital  filter  in  place.  Simulator  results  showed  the 


non-linear  filter  had  no  noticeable  effect  on  the  lateral -directional 


modes  of  interest.  Open  loop  data  with  the  newly  computed  gains  and 


noise  filters  were  within  measurement  accuracy  of  previous  simulator 


results.  No  further  structural  resonance  problems  were  encountered  on 


the  ground  or  in  flight. 


Flight  Test  Results  -  0.4  Mach.  The  control  laws  derived  from 


eigenstructure  assignment  produced  unflyable  characteristics  at  0.4 


Hach.  Although  some  of  the  control  laws  were  better  than  others,  all 


% 


I 


s 


suffered  fro*  vary  poor  dutch  roll  characteristic*,  in  soae  cases 
divergent.  This  was  particularly  puzzling  because  initial  flight  teat 
results  at  O.S  Mach  agreed  with  siaulator  results. 

The  problea  was  finally  traced  to  inaccurate  assuaptions  in  the 
YA-7D  linear  aodel  used  in  the  derivation  of  the  control  lawa.  In  the 
linear  aodel,  it  was  asauaed  the  lateral  accaleroaeter  was  4.78  feet 
forward  of  the  center  of  gravity.  Corrections  to  the  lateral 
accaleroaeter  output  were  Bade  in  the  aodel  to  account  for  the  yaw 
acceleration  which  would  be  sensed  by  the  accaleroaeters.  However,  not 
aodeled  in  the  linear  aodel  or  SAFTD  was  the  vertical  position  of  the 
acceleroaeters,  approximately  17  inches  below  the  center  of  gravity. 
Thus,  in  addition  to  sensing  yaw  acceleration,  the  acceleroaeters  alao 
sensed  roll  acceleration.  Although  the  aoaent  ara  waa  relatively  saall, 
roll  acceleration  can  be  quite  high,  aore  than  100  degrees/sec2.  The 
SAFTD  was  used  to  evaluate  affects  on  the  dutch  roll  of  having  the 
acceleroaetera  17  inches  below  the  CG.  Results  are  shown  in  Table  XI. 
While  the  dutch  roll  period  and  phi /beta  ratio  were  not  changed 
significantly,  daaping  was  greatly  reduced.  The  effects  at  0.4  Hach 
were  such  greater  than  at  0.6  Mach.  This  was  to  be  expected  since  the 
lateral  acceleration  gains  were  aore  than  five  tiaes  greater.  Mora 
iaportantly,  the  control  laws  at  0.4  Mach  were  now  found  to  be  very 


sensitive  to  changes  in  airspeed.  A  20  knot  inert 


enough  to 


drive  the  dutch  roll  divergent  in  soae  cases.  The  dutch  roll  response 
was  also  no  longer  a  second  order  function,  but  becaae  non-linear. 
These  effects  corresponded  closely  to  flight  test  results  at  0.4  Mach. 


Table  XI 


Simulator 

Result*  Showing  Effect* 
Poaition  17  Inches 

of  Lateral 
Below  CG 

Acceleroaatar 

Hach 

Design 
sag  / 

<0/e>d 

phase 

wd 

Id 

<#/0>d 

0.4 

0 

/ 

0 

3.0 

0.32 

0  arg(O) 

1.3 

/ 

0 

2.5 

0.28 

2.2  arg(14) 

1.5 

/ 

60 

2.7 

0.25 

1.8  arg(60) 

1.5 

/ 

120 

2.9 

0.10 

1.5  arg(114> 

3.0 

/ 

0 

2.3 

0.28 

3.3  arg(O) 

3.0 

/ 

60 

2.4 

0.22 

2.6  arg(53) 

3.0 

/ 

120 

2.9 

0.01 

2.3  arg(98) 

0.6 

0 

/ 

0 

3.2 

0.35 

0  arg(O) 

1.5 

/ 

0 

2.8 

0.27 

1.8  arg<23) 

1.5 

/ 

60 

3.0 

0.29 

1.9  arg(57) 

1.3 

/ 

120 

3.1 

0.27 

1.5  arg<86) 

3.0 

/ 

0 

2.7 

0.32 

3.9  arg(22> 

3.0 

/ 

60 

2.9 

0.29 

3.4  arg(55) 

3.0 

/ 

120 

3.2 

0.26 

2.9  arg<88) 

Those  control  law*  found  dlvargent  in  flight  war*  also  found  divargant 
in  tha  aiaulator. 

The  liaitad  acopa  of  tha  flight  teat  investigation  did  not  parait 
correcting  the  gain*  for  tha  actual  lateral  acceleroaatar  position  and 
retaating.  However,  tha  close  correlation  between  aiaulationa  with  tha 
correct  accalaroaater  position  and  flight  teat  raaulta  at  0.4  Hach 
indicat*  the  problem  was  correctly  identified. 


The  rather  dramatic  effects  of  the  lateral  accelerometer  vertical 
position  error  in  the  model  can  be  analyzed  by  examination  of  the 
equations  for  the  YA-7D  lateral  acceleration  (Ny).  The  lateral 
acceleration  at  the  CG  modeled  in  Section  III  for  the  YA-7D  at  0.4  Mach 
and  15,000  feet  is: 

My  *  0.8Sr  -  48.00  ♦  0.55p  ♦  2.476aa  ♦  13.37<5r  (52) 

Ny  is  primarily  a  function  of  sideslip  and  rudder  deflection.  At  the 
accelerometer,  4.8  feet  forward  of  the  CG,  a  component  of  yaw 
acceleration  is  sensed.  The  accelerometer  output  becomes: 

Ny  =  -0.47r  -  41.00  ♦  0.098p  -  2.55daa  ♦  2.84<Sr  (53) 

The  accelerometer  output  is  now  primarily  a  function  of  sideslip,  and 
the  signal  may  be  used  in  place  of  sideslip  with  good  results.  This  was 
the  relationship  used  in  the  linear  model.  At  0.S  Mach  the  relationship 
is  similar. 

If  the  accelerometer  is  now  moved  4.8  feet  forward  of  the  CG  and  17 
inches  down,  it  will  now  sense  a  component  of  roll  acceleration  also. 

The  accelerometer  output  becomes: 

Ny  =  -1.75r  -  24.50  ♦  2.25p  ♦  14.73<iaa  -  0.46dr  (54) 

The  output  is  now  primarily  a  function  of  sideslip  and  aileron 
deflection.  If  this  is  not  modeled,  the  feedback  gains  will  not  produce 
the  desired  results  as  was  evident  in  both  simulator  and  flight  test 
results.  Although  the  gains  were  small  enough  so  significant  errors 


were  not  introduced  at  0.6  Mach,  inatabilitiea  were  introduced  at  0.4 


Mach. 

Flight  Teat  Results  -  0.6  Mach.  Although  the  problea  with  the 
lateral  accelerometer  position  was  also  present  at  0.6  Mach,  it  had  much 
less  effect  because  the  control  systea  gains  were  auch  saaller.  The 
open  and  closed  loop  testing  described  in  Appendix  F  were  accoaplished. 

Open  Loop  Test  Results.  Open  loop  tests  were  performed  to 
verify  the  flight  control  laws  produced  the  desired  aircraft  response. 
The  closed  loop  lateral -directional  eigenstructure  (as  measured  by  the 
roll  aode  time  constant,  dutch  roll  frequency  and  damping,  and  (0/0) d  > 
were  measured  as  described  in  Appendix  F.  Results  are  shown  in  Table 
XII.  Flight  test  results  generally  agreed  within  15  percent  of  design 
conditions  and  simulator  results.  This  was  an  excellent  correlation 
considering  wind  tunnel  aerodynamic  data  used  to  derive  the  control  laws 
was  in  error  as  shown  in  Table  XIII.  Flight  test  results  indicated  the 
aircraft  had  a  much  larger  N0  than  estimated  by  the  aerodynamic  data. 

It  was  also  noted  by  the  project  test  pilots  that  the  aircraft  had 
noticeable  adverse  yaw  not  predicted  by  the  aerodynamic  data. 
Nevertheless,  the  control  laws  derived  in  the  previous  sections  were 
successful  in  attaining  the  desired  open  loop  response  in  both  ground 
simulations  and  flight  testing. 

An  unusual  characteristic  of  the  control  laws  was  discovered  during 
open  loop  testing.  The  control  laws  with  l0/0ld  of  0  and  (0/0>d  of  1.5 
arg(0)  required  unnatural  aileron  control  during  sideslips.  Although 
unusual  and  in  violation  of  paragraph  3. 3. 6. 3  of  MIL-F-6785C  (9),  this 
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Table  XII 


Flight  Teat  Lateral 

-Directional 

Characteristics  - 

-  0.6  Mach 

Design  <0/0)d* 
(nag/phase  ang) 

Wd 

Id 

Tr2 

Pss'Fs 

<#/0>d 

(nag/ phase  ang) 

Unaugaented 

2.6 

0.16 

0.50 

N/A 

2.6  /  58 

DFCS  on 

2.9 

0.45 

0.36 

8.3 

2.7  /  63 

0/0 

3.4 

0.40 

0.44 

9.7 

0/0 

1.5  /  0 

3.4 

0.40 

0.42 

8.3 

1.5  /  14 

1.5  /  60 

3.4 

0.32 

0.36 

8.0 

1.3  /  60 

1.5  /  120 

3.6 

0.35 

0.35 

7.9 

1.2  /  113 

3.0  /  0 

3.2 

0.36 

0.44 

7.9 

2.8  /  0 

3.0  /  60 

3.2 

0.35 

0.44 

8.9 

2.8  /  52 

3.0  /  120 

3.4 

0.32 

0.44 

8.7 

2.6  /  100 

1  Note*.  Design  conditions  were  rr  *  0.25.  w<j  3 

2  Note:  Includes  effects  of  prefilters  in  CAS. 

3.0.  and 

Id  3  0.40. 

Tabla  XIII 

Comparison  of 

Math  Model 

with  Flight  Teat  Results 

at  0.6  Mach 

Tr 

Wd  Id 

<0/6>d 

Math  Model 

0.43 

2.1  '  0.15 

3.4  arg(40) 

Flight  Test 

0.50 

2.8  0.16 

2.6  arg(58) 

did  not  produce  any  objectionable  characteristics  during  the  closed  loop 

teaks . 


Based  on  open  loop  data,  all  flight  control  laws  would  be  expected 
to  produce  satisfactory  flying  qualities.  The  dutch  roll  frequency  was 
higher  than  expected  and  the  damping  ratio  slightly  lower,  but  were 
still  near  ainiaua  requirements  of  HIL-F-8785C  (9).  The  effective  roll 
mode  time  constants,  which  included  effects  of  lag  prefilters  in  the 
CAS,  and  stick  force  command  gradient,  were  satiaf actory . 

Closed  Loop  Test  Results.  The  three  project  pilots  evaluated 
each  control  law  using  three  closed  loop  tasks  a  minimum  of  two  timea. 
The  exact  procedures  are  described  in  detail  in  Appendix  F.  Evaluations 
of  the  unaugmented  aircraft  and  standard  DIGITAC  DFCS  were  also 
performed  for  comparison.  The  raw  data  are  given  in  Tables  El  through 


Figures  El  through  E3  present  data  comparing  pilot  repeatability 
and  pilot-to-pilot  repeatability.  In  general,  each  pilot  was  consistent 
in  his  ratings  within  plus  or  minus  one  Cooper-Harper  rating.  Also,  the 
inter  pilot  ratings  generally  agreed  within  the  sane  tolerance. 

The  following  narrative  will  discuss  each  closed  loop  task  and 
pilot  ratings  in  detail,  as  well  as  give  an  overall  summary. 

All  pilots  found  HUD  Task  1  the  most  demanding  and  difficult  to 
perform.  As  the  pilot  increased  the  gain  and  excited  the  dutch  roll, 
capturing  the  target  became  extremely  difficult.  All  pilots  also  agreed 
they  had  to  reduce  their  gain  just  to  accomplish  the  task.  Host 
complaints  centered  on  mechanization  of  the  taak.  The  target  was  not 
inertially  stabilized  and  moved  toward  the  pipper  as  the  aircraft  was 


XT. 


banked  without  regard  to  aircraft  noae  position.  The  illogical  target 
•oveaent  caused  the  task  to  be  ineffective,  except  for  saall  movements 
of  less  than  five  milliradians  (mils)  around  the  pipper.  In  principle, 
the  task  was  perhaps  the  best  for  evaluating  lateral  and  directional 
characteristics,  but  all  pilots  ranked  it  lowest  of  the  three  closed 
loop  tasks  for  its  ability  to  discriminate  flying  qualities.  Due  to 
software  development  delays,  the  task  was  not  available  until  late  in 
the  flight  test  period,  and  only  limited  evaluations  were  performed. 
Insufficient  time  remained  to  correct  the  software  problems  and  refly 
the  testa. 

Despite  technical  problems  with  the  task,  pilots  were  able  to 
identify  differences  in  the  control  system  designs.  Figure  12  on  the 
following  page  shows  pilot  rating  as  a  function  of  phi/beta  magnitude 
and  phase  angle.  Pilots  were  very  consistent  in  ranking  zero  deg.-ee 
phase  angle  as  the  best,  with  120  degrees  as  the  worst.  Results  show  no 
consistent  variation  of  pilot  rating  with  magnitude,  although  none  of 
the  pilots  liked  a  phi/beta  magnitude  of  zero. 

Many  pilot  comments  centered  on  residual  yaw  oscillations. 

Although  adverse  yaw  was  not  a  complaint,  there  was  a  sufficient  amount 
to  excite  the  dutch  roll.  The  pilots  felt  the  damping  ratio  was  too 
low,  and  that  increased  damping  would  have  improved  the  Cooper-Harper 
ratings. 

HUD  Task  2  was  not  expected  to  show  much  variation  in  the  pilot 
ratings  since  it  primarily  was  a  roll  task  and  all  control  laws  produced 
similar  roll  mode  time  constants  and  roll  response  to  stick  force 
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Figure  12.  HUD  Task  1  -  Cooper-Harper  Ratings 
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gradients.  It  was  expected  all  control  lawa  would  produce  satisfactory 


flying  qualities.  Results  are  shown  in  Figure  13  on  the  following  page. 
The  pilots  were  unanimous  in  objecting  to  the  roll  response.  Aircraft 
response  was  abrupt,  causing  poor  predictability.  As  pilot  gain 
increased,  cockpit  roll  and  yaw  accelerations  became  very  high  which 
made  tracking  difficult.  Again,  project  pilots  felt  increased  dutch 
roll  damping  would  have  helped.  This  conclusion  was  surprising  because 
the  dutch  roll  damping  ratio  was  near  the  HIL-F-8785C  requirement,  yet 
all  pilots  felt  the  damping  ratio  was  too  low. 

None  of  the  control  laws  produced  satisfactory  flying  qualities  in 
HUD  Task  2,  including  the  standard  DIGITAC  DFCS,  even  though  the  stick 
gradients  and  roll  node  time  constant  were  good.  Pilots  A  and  C  even 
preferred  the  unaugmented  aircraft  most  because  of  its  good 
predictability,  although  stick  forces  were  high.  Two  possible 
explanations  were  explored. 

Excessive  time  delay  was  first  suspected  as  causing  the 
degradations  in  pilot  ratings.  All  control  laws  had  an  effective  time 
delay  of  approximately  110  milliseconds  (msec).  Approximately  60  msecs 
was  due  to  hydraulic  actuators,  the  digital  FCS,  and  higher  order 
effects.  Host  of  the  delay  was  due  to  the  three  rad/sec  linear 
prefilter  in  the  CAS.  The  time  constant  for  the  CAS  prefilter  was 
reprogrammed  to  be  cockpit  selectable  to  0.300,  0.125,  and  0.025 
seconds.  While  reducing  the  prefilter  time  constant  reduced  the 
effective  time  delay,  it  also  reduced  the  effective  roll  node  time 
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constant.  This  increased  the  roll  response,  but  the  effective  roll  sode 
tine  constant  was  still  greater  than  0.25  seconds.  HUD  Task  2  was  then 
reflown  to  evaluate  the  different  prefilters.  Reducing  the  tiae 
constants  of  the  prefilter  reduced  the  tiae  delay  below  60  a sec,  but 
degraded  pilot  ratings  also.  With  the  0.025  second  prefilter  the 
aircraft  becaae  very  prone  to  lateral  pilot  induced  oscillations  leading 
to  Cooper-Harper  ratings  as  high  as  eight.  Cockpit  accelerations  were 
worse  and  roll  response  even  aore  abrupt. 

Investigation  then  focused  on  stick  dynaaica.  The  YA-7D  stick  was 
a  force  coaaand  stick  with  little  lateral  aoveaent  because  of  the  saall 
cockpit.  Stripchart  data  showed  significant  variation  in  stick  force 
even  when  the  pilot  was  atteapting  to  hold  the  stick  steady  during 
rolling  aaneuvera.  It  is  believed  the  coabination  of  a  sensitive  stick 
with  high  cockpit  accelerations  produced  the  poor  pilot  ratings.  For 
exaaple,  if  the  pilot  put  in  a  five  pound  force,  the  aircraft  would 
produce  lateral  accelerations  in  the  cockpit  on  the  order  of  three 
ft/sec?.  Aasuaing  the  stick  and  pilot's  ara  weighed  10  pounds,  this 
produced  a  force  variation  of  one  pound,  or  20  percent.  When  pilots 
specifically  looked  for  variations  in  stick  force,  all  coaaented  on  the 
large  effect  of  the  cockpit  accelerations.  Pilot  B  reflew  several  of 
the  control  laws  with  both  hands  on  the  stick  to  ainiaize  the  effects  of 
stick-ara  dynaaica.  Even  though  he  was  using  the  aost  sensitive  stick 
(with  a  CAS  prefilter  tiae  constant  of  0.025  seconds)  he  noted  a  aarked 
iaproveaent.  The  Cooper-Harper  rating  for  HUD  Task  2  dropped  froa  an 
eight  to  a  three.  Thus,  it  is  felt  the  stick  dynaaica,  rather  than  a 
problea  in  the  control  laws  produced  the  low  pilot  ratings  in  the  roil 


response.  In  fact,  the  pilots  felt  stick  dynamics  were  the  most 
important  factor  in  the  DIGITAC's  lateral  handling  qualities. 


Results  in  Figure  13  show  little  variation  of  pilot  rating  with 
phi/beta  magnitude  or  phase  angle.  With  a  phi/beta  magnitude  ratio  of 
three,  slight  roll  oscillations  were  excited  by  adverse  yaw,  but  they 
were  well  damped  and  did  not  significantly  degrade  pilot  ratings. 

Pilot  opinions  of  the  air-to-air  tracking  task  differed.  Pilot  A 
disliked  it  and  felt  he  was  not  able  to  increase  his  gain  sufficiently 
because  of  pendulum  effects  inherent  in  the  sight.  However,  pilot  C 
liked  it  the  best  of  all  closed  loop  tasks  and  felt  it  most  illustrated 
the  differences  in  the  control  systems.  This  is  reflected  in  the 
results  shown  in  Figure  14  on  the  following  page.  Pilot  C  was  able  to 
use  the  task  to  break  out  differences  in  the  control  laws  very 
consistently.  All  pilots  were  consistent  in  disliking  the  phi/beta 
phase  angle  of  120  degrees,  which  agreed  with  test  results  from  HUD  Task 
1.  Again,  there  was  no  strong  breakout  of  pilot  rating  with  phi/beta 
magnitude  ratio.  Dutch  roll  damping  ratio  was  sufficient  in  this  task. 
Many  of  the  pilot's  complaints  centered  on  the  abrupt  roll  response, 
which  has  already  been  discussed. 

Figure  15  on  page  82  shows  Cooper-Harper  ratings  for  the  three 
closed  loop  tasks  using  the  average  of  all  three  pilots'  ratings. 

Results  are  similar  to  those  previously  discussed.  The  phi /beta 
magnitude  ratio  did  not  have  a  noticible  effect  on  the  flying  qualities, 
but  the  pilot  ratings  consistently  degraded  for  a  phi /beta  phase  angle 
of  120  degrees.  HUD  Task  1  showed  a  phase  angle  of  zero  degrees  to  be 
the  best,  while  the  other  tasks  did  not  show  a  breakout  between  zero  and 
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Figura  14.  Air-to-Air  Tracking  -  Coopar-Harpar  Ratings 
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60  degrees.  The  pilots  felt  with  s  phsse  angle  of  120  degrees  the  dutch 
roll  tended  to  get  out  of  phase  with  the  pilot's  inputs. 

The  drasatic  effects  of  lateral  stick  dynamics  influenced  pilot 
ratings  in  all  tasks.  It  was  felt  that  breakout  of  pilot  rating  with 
phi/beta  aagnitude  and  phase  angle  say  have  actually  been  greater  than 
reported,  but  poor  stick  dynaaics  overshadowed  the  aoderate  differences 
in  aircraft  response  created  by  the  control  systea  designs.  For 
exaaple,  the  pilots  felt  the  aircraft  response  with  phi/beta  of  zero  was 
the  aost  abrupt,  and  overly  sensitive.  Zt  was  difficult  to  identify  how 
auch  of  this  was  due  to  the  natural  response  of  the  aircraft  and  how 
auch  due  to  stick  dynaaics. 

The  effects  of  lsteral  stick  dynaaics  also  illustrated  that 
achieving  good  flying  qualities  is  auch  aore  than  coapliance  with  open 
loop  requireaents.  Guidance  given  in  HIL-F-8785C  is  insufficient  to 
insure  satisfactory  handling  qualities.  All  coaponents  between  the 
pilot's  input  and  the  aircraft  response  aust  be  considered  in  detail. 

Evaluation  of  Eloenatructure  Assignaent 

The  results  of  the  flight  tests  verified  eigenstructure  assignaent 
as  a  valid  technique  to  derive  flight  control  laws.  Despite  gross 
uncertainties  in  the  aerodynaaic  data  and  errors  introduced  by  a  digital 
flight  control  systea,  the  control  laws  produced  the  desired  open  loop 
responses.  The  ease  with  which  new  control  laws  could  be  derived  when 
changes  had  to  be  incorporated  proved  to  be  the  aost  positive  feature  of 
eigenstructure  assignaent.  Using  the  interactive  coaputer  software  in 
Appendix  A,  derivation  of  14  new  control  laws  due  to  the  addition  of  the 
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noise  filters  took  less  than  one  hour.  Changes  in  sass  properties, 
aerodynamic  data,  ate.,  could  be  incorporated  as  aasily. 

Several  practical  probleas  arose  which  illustrated  serious 
deficiencies  with  eigenstructure  assignment.  Lateral  acceleration 
feedback  was  very  noisy  and  sensitiva  to  center  of  gravity  position. 
Although  the  structural  resonance  and  position  problems  with  the  lateral 
acceleration  signal  were  not  unique  to  eigenstructure  assignment,  they 
did  show  inaccurate  modeling  coupled  with  high  gains  can  lead  to 
unexpected  results. 

The  major  disadvantage  of  eigenstructure  assignment  was  the  lack  of 
intuitive  interpretation  of  the  control  laws,  making  specific  solutions 
to  problems  very  difficult.  Many  hours  of  analysis  and  simulation  were 
required  to  identify  the  vertical  displacement  of  the  lateral 
accelerometer  as  the  cause  for  the  poor  performance  at  0.4  Mach. 
Furthermore,  it  waa  recognized  early  that  dutch  roll  damping  was 
insufficient,  but  increased  damping  could  not  be  accomplished  with  an 
increase  of  a  single  gain,  as  might  be  done  in  a  conventional  control 
systam.  The  interaction  of  all  feedback  loops  must  be  considered  to 
produce  the  desired  response.  Results  from  adjustment  of  a  single  gain 
were  unpredictable  making  fine  tuning  of  the  control  laws  based  on  pilot 


comments  impossible 


This  research  provided  such  practical  experience  with 
eigenstructure  aseignaent.  The  aost  attractive  aspect  of  the  aethod  is 
the  ease  with  which  gains  can  be  derived.  Once  a  desired  eigenatructure 
is  selected,  and  a  state  space  aodel  developed,  gain  calculation  is 
quick  and  easy.  Developaent  of  even  coaplicated  state  space  aodels  is 
relatively  easy,  since  a  building  block  approach  can  be  used  to  expand 
the  aodel.  Changes  to  the  aodel  or  desired  sigenstructure  can  be 
incorporated  quickly,  and  new  gains  calculated. 

Eigenstructure  aseignaent  was  successful  in  achieving  the  desired 
lateral -directional  eigenstructures.  Despite  a  poor  aerodynaaic  aodel, 
and  errors  in  aodeling  the  lateral  acceleroaeter  output,  the  desired 
eigenvalues  and  eigenvectors  were  achieved.  The  gains  produced  the 
desired  (0/0>d  over  a  wide  range  of  aegnitudes  and  phase  angles. 

Both  aathaatical  analysis  and  flight  test  data  showed  the  handling 
qualities  of  the  YA-7D  to  be  relatively  insensitive  to  changes  in 
(0/0) d<  This  gives  the  designer  leeway  in  selecting  eigenstructures 
which  can  yield  different  gains  and  yaw-roll  coupling.  However,  the 
coaplexity  of  the  feedback  laws  pose  serious  practical  difficulties. 

The  control  laws  lack  intuitive  interpretation,  aaking  isolation  of 
handling  qualities  probleas  difficult.  Failure  analysis,  which  was 
never  addressed  in  this  Halted  study,  would  also  be  coaplicated.  The 
YA-7D  lateral -directional  flying  qualities  showed  no  iaproveaent  over 
those  of  the  standard  DFCS,  which  used  conventional  control  laws.  Thesa 


problems,  although  certainly  not  insurmountable,  would  aeea  to  aaka 
eigenetructure  Impractical  for  conventional  aircraft. 

Eigenatructure  assignment  may  offer  significant  advantages  for 
unconventional  aircraft.  With  an  additional  control,  such  as  direct 
side  force,  derivation  of  control  laws  with  conventional  methods  no 
longer  offers  a  clear  advantage.  Analysis  with  conventional  techniques 
would  be  very  complex,  while  derivation  of  control  laws  with 
eigenatructure  assignment  would  be  hardly  more  complicated  than  for  the 
two  control  case. 


Most  of  the  analysis  in  this  report  is  valid  and  directly 
applicable  with  the  addition  of  direct  side  force.  Not  only  does  this 
additional  control  allow  one  more  of  the  eigenvector  components  to  be 
exactly  placed,  but  it  allows  modification  of  components  which  could  not 
be  effected  by  the  two  conventional  controls.  The  clearest  example  of 
this  is  the  dutch  roll  eigenvector.  It  was  shown  in  Section  V  that  the 


dutch  roll  eigenvector  was  almost  totally  determined  by  the  ratio  of 
only  two  components,  <0/9)4.  However,  in  that  analysis  it  was  assumed 
Yv  and  consequently  <r/0)<i  could  not  be  greatly  changed  by  the  closed 
loop  controls.  With  direct  sideforce  control  this  assumption  is  no 
longer  true.  The  designer  can  now  control  both  <0/0)0  and  <r/0)d,  which 
allows  the  handling  qualities  to  be  modified  in  a  manner  not  possible 
with  conventional  aircraft.  Sideslip  angle  can  be  controlled  directly, 
and  turn  coordination  can  be  improved.  It  also  allows  control  of  the 
lateral  acceleration,  which  nay  improve  handling  qualities  by  reducing 
pilot  station  acceleration.  In  fact,  the  aircraft  can  now  be  made  to 
roll  about  any  desired  axis.  It  is  not  clear  how  to  select  <r/0)q  to 


achieve  the  best  coebination  of  sideslip  angle  and  lateral  acceleration 
control.  And  direct  sideforce  control  makes  selection  of  (0/0)4  even 
more  confusing.  This  ares  deserves  future  study  to  determine  if 
handling  qualities  can  be  usefully  improved  with  direct  side  force 
control,  and  the  relationship  of  <r/0>4  to  such  improvements. 

The  final  conclusions  from  this  study  are  certainly  not  new  to 
experienced  aircraft  designers  and  not  unique  to  eigenstructure 
assignment,  but  are  important  to  emphasize  nonetheless.  Satisfactory 
handling  qualities  can  only  be  achieved  by  considering  the  aircraft  as  a 
whole.  Everything  that  affects  the  pilot  must  be  considered.  In  the 
case  of  the  YA-7D,  cockpit  acceleration  coupled  with  the  stick  design  to 
produce  unusual  lateral  dynamics.  The  study  also  pointed  out  the 
importance  of  flight  testing.  Regardless  of  the  care  used  in  the  math 
model  and  control  law  design,  satlsfsctory  handling  qualities  can  only 
be  guaranteed  through  carefully  planned  and  conducted  flight  tests. 

Eigenstructure  assignment  is  a  valid  technique  for  derivation  of 
control  laws.  Although  of  questionable  value  for  conventional  * 
airplanes,  it  promises  more  usefulness  in  unconventional  aircraft. 


APPENDIX  A.  FORTRAN  Program  for  Eigenatructure 
Assignment  Using  Output  Feedback 

This  FORTRAN  ispleaentation  of  eigenatructure  assignment  follows 
the  description  in  Section  II.  International  Mathematics  and  Statistics 
Library  (IMSL)  subroutines  were  used  to  calculate  (Xil  -  A)'1  and  the 
generalized  inverse,  L*i  as  well  as  matrix  products  and  additions.  The 
only  difficulty  is  in  incorporating  complex  eigenvalues  and  eigenvectors 
into  the  algorithm.  This  can  be  done  by  converting  certain  complex 
operations  to  real  ones  and  using  available  library  routines  to  perform 
complex  operations  as  described  below. 

Accuracy  is  usually  best  if  <Xj,I  -  A)-1  is  calculated  directly 
using  complex  arithmetic.  It  is  then  fairly  easy  to  multiply  by  B  to 
produce  the  complex  matrix  Li  «  (hi  -  A) "^B  .  Up  to  this  point  the 

operations  have  been  the  sane  for  complex  and  real  eigenvalues,  although 
the  routines  used  will  of  course  differ.  For  real  eigenvalues  the 
generalized  inverse,  L*  is  found  next  using  a  singular  value 
decomposition.  However,  for  complex  matrices  library  subroutines  do  not 
exist  to  find  the  generalized  inverse.  Even  if  available  it  would  still 
not  be  desirable  because  it  would  then  not  be  possible  to  specify  only 
the  real  portion  or  only  the  imaginary  portion  of  a  given  eigenvector 
component.  To  get  around  this  problem  a  real  matrix,  Li'  is  defined 


Li' 


a 


Re<Li>  I  -Im(Li> 
Im(Li)  I  Re(Li> 


<55> 
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where  Re(Li)  end  Im<Li.)  are  the  real  and  laaglnary  portions  of  Li, 
respectively.  Li  la  an  n  X  a  complex  matrix  ao  Li'  is  a  2n  X  2m  real 
matrix.  The  corresponding  complex  eigenvector,  vj,,  becomes 

yi'  *  t  vir  V2R  ...  VnR  yij  V2I  . . .  Vni  1  <56) 

where  ynR  and  yni  are  the  ntti  real  and  imaginary  components  of  the 
eigenvector,  yi .  Thus,  yj.'  is  of  length  2n.  The  rows  in  Li' 
corresponding  to  the  rows  in  yi'  which  are  to  specified  are  then  used  to 
find  the  generalized  inverse.  The  achievable  eigenvector  and  wi  are 
found  exactly  as  for  real  eigenvalues,  except  of  course  they  are  both  of 
length  2n.  They  can  then  be  divided  into  real  and  imaginary  portions 
again  so  each  complex  eigenvalue  will  take  up  two  rows  of  wi  and  V 
instead  of  one.  Therefore  the  calculation  of  F  is  not  affected. 

The  FORTRAN  program  listed  below  uses  eigenstructure  assignment  to 
calculate  F  assuming  there  are  two  controls  and  four  output  quantities 
which  are  linear  combinations  of  the  state. 


O  O  O  O  O  Ci  o 


C  PROGRAM  NAME:  EIGEN4 

C  AUTHOR:  MICHAEL  COSTIGAN 

C  DATE:  16  JUL  66 

C  PURPOSE:  TO  CALCULATE  A  FEEDBACK  MATRIX  F 

C  THAT  WILL  ASSIGN  CLOSED  LOOP  EIGENVALUES 

C  AND  EIGENVECTORS  USING  OUTPUT  FEEDBACK. 

C  INPUTS: 

C  TITLE  -  30  CHARACTER  TITLE  PRINTED  IN  OUTPUT 

C  N  -  ORDER  OF  THE  SYSTEM 

C  A  -  AN  N  X  N  STATE  MATRIX 

C  B  -  AN  N  X  2  CONTROL  MATRIX 

C  C  -  AN  N  X  4  OUTPUT  MATRIX 

C  THE  ABOVE  MATRICES  ARE  READ  FROM  FILE  STATE 
C  LABEL  -  30  CHARACTER  LABEL  PRINTED  IN  OUTPUT 

C  LAMDR  -  DESIRED  COMPLEX  DUTCH  ROLL  E-VALUE 

C  LAMSP  -  DESIRED  SPIRAL  MODE  E-VALUE 

C  LAMRS  -  DESIRED  ROLL  MODE  E-VALUE 

C  VDDRR  -  REAL  PORTION  OF  DESIRED  DUTCH  ROLL  E-VECTOR 

C  VDDRI  -  IMAG  PORTION  OF  DESIRED  DUTCH  ROLL  E-VECTOR 

C  VDSP  -  DESIRED  COMPONENTS  OF  SPIRAL  E-VECTOR 

C  VDRS  -  DESIRED  COMPONENTS  OF  ROLL  E-VECTOR 

C  IDRR  -  N  X  1  VECTOR  IDENTIFYING  ROWS  OF  REAL  DUTCH 

C  ROLL  E-VECTOR  TO  BE  SPECIFIED 

C  IDRI  -  N  x  1  VECTOR  IDENTIFYING  ROWS  OF  IMAG  DUTCH 

C  ROLL  E-VECTOR  TO  BE  SPECIFIED 

C  I DSP  -  N  X  1  VECTOR  IDENTIFYING  ROWS  OF  SPIRAL 

C  E-VECTOR  TO  BE  SPECIFIED 

C  IDRS  -  N  X  1  VECTOR  IDENTIFYING  ROWS  OF  ROLL 

C  E-VECTOR  TO  BE  SPECIFIED 

C  THE  ABOVE  VARIABLES  IDENTIFY  THE  DESIRED  EIGENVALUES 

AND  EIGENVECTORS  AND  ARE  READ  FROM  FILE  EIGENV 


REAL  A<20,20) ,B(20,2) ,C(5,20) , LAMSP, LAMRS, VDDRR C20) , 
•VDDRI (20) , VDSP (20) ,VDRS(20) ,VDRR(20) ,VDRI(20) ,VSP(20) , 
•VRS(20) , V(20,5) ,AAUG(20,20) ,BAUG(20,20) ,WK(350) ,F(2,5) 
•,W(2,5) 

INTEGER  N,ISUM(5) , IDRR (20) , IDRI (20) , I DSP (20) , IDRS (20) 
COMPLEX  LAMDR , EVEC ( 20 , 20 ) , EVAL ( 20 ) , ZN , EVECA (20,20), 
•EVALA(20) ,ZNA 
CHARACTER  TITLE»30,LABEL*30 

READ  INPUT  DATA  FROM  FILES 

OPEN (4, FILE*' STATE' ,STATUS»'OLD' ) 

REWIND (4) 

READ (4,*) TITLE 


o  o  ana  nan  ana  ana  ana  ana 
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READ(4,*)N 

READ(4, •)  < ( A < I , J)  ,  J*1,N)  ,  1*1  ,N) , ( <B< I , J) , J  =  1 ,2  ) , 
*1*1, N) , < (C(I, J)  ,  J*1,N)  ,  1*1,4) 

CLOSE (4) 

0P£N<4,FILE*'EIGENV' ,STATUS*'OLD' ) 

REWIND (4) 

READ < 4,*) LABEL 

READ ( 4 , • ) LAHDR , LAMSP , LAHRS , (IDRR(I) ,2*1,N> , 

•<IDRI(I) , 1*1 ,N) , (IDSP(I) , 1*1 ,N> , (IDRS(I) , I*1,N) 

DO  10  1*1, N 

ISUM  < 1 ) *ISUM < 1 ) +IDRR ( I ) 

ISUH<2) *ISUN(2) +IDRI ( I ) 

ISUM<3)*ISUM(3)*IDSP< I) 

10  ISUM(4)*ISUM<4) *IDRS( I ) 

READ<4, •) (VDDRR(I) ,1*1, ISUM (1) ) , (VDDRI (I) ,1*1 , ISUM (2) ) 
*, (VDSP(I) ,I*1,ISUN(3) ) , <VDRS(I) ,1*1, ISUM (4) ) 

CLOSE (4) 

CALCULATE  SPIRAL  MODE  E-VECTOR 

CALL  REVEC (A,B,N,W,IDSP,VDSP, LAMSP, VSP, ISUM (3) ,3) 
CALCULATE  ROLL  MODE  E-VECTOR 

CALL  REVEC (A,B,N,W,IDRS,VDRS,LAMRS,VRS, ISUM (4) ,4) 

CALCULATE  DUTCH  ROLL  MODE  E-VECTOR 

CALL  DRVEC <  A , B , N , W , IDRR , IDRI , VDDRR , VDDRI , LAMDR , VDRR , 

• VDRI , ISUM < 1 ) , ISUM < 2 ) ) 

FORM  MATRIX  OF  EIGENVECTORS 

DO  20  1*1, N 
V<I,1)*VDRR<I) 

V(I,2)*VDRI(I) 

V(I,3)*VSP<I) 

20  V(I,4)*VRS(I> 

CALCULATE  FEEDBACX  MATRIX,  F 

CALL  FCALC(C,W,N,V,F> 

CALCULATE  AUGMENTED  A  MATRIX 

CALL  VMULFF(F,C,2,4,N,2,5, AAUG,20, IER) 

CALL  VMULFF ( B , AAUG , N , 2 , N , 20 , 20 , B AUG , 20 , I ER ) 

DO  30  1*1, N 
DO  30  J*1,N 

30  AAUG <I,J)*A(I,J) +BAUG < I , J ) 


CALCULATE  E-VALUES  AND  E-VECTORS  OF  AUGMENTED  SYSTEM 


OOO  O  O  o  o  o  o 


CALL  EIGRF <AAUG,N,20,2,EVAL,£VEC,20,WK, IER) 

IF(WK(1)  .GT.  1. ) PRINT* , 'WARNING  -  EIGENVECTORS' 

*,'  OF  AUGMENTED  A  MATRIX  MAY  NOT  BE  ACCURATE' 

CALCULATE  E-VALUES  AND  E-VECTORS  OF  UNAUGMENTED  SYSTEM 

CALL  EIGRF<A,N,20,2,EVALA,EVECA,20,WK,IER) 

IF(WK(1>  .GT.  1 . ) PRINT* , ' WARNING  -  EIGENVECTORS' 

•'OF  UNAUGMENTED  SYSTEM  MAY  NOT  BE  ACCURATE' 

NORMALIZE  E-VECTORS 

DO  40  J*1,N 
K*1 
M*1 

35  ZN*EVEC<K,J) 

IF <CABS(ZN)  .LT.  l.E-4  .AND.  K.LE.N)  THEN 
K*K*1 
GO  TO  35 

ELSEIFUC  .GT.  N)  THEN 

PRINT* , 'WARNING  -  THERE  IS  A  ZER  E-VECTOR' 

CALL  EXIT 
ENDIF 

45  ZNA*EVECA<M,J) 

IF (CABS(ZNA)  .LT.  l.E-4  .AND.  M.LE.N)  THEN 
M*M*1 
GO  TO  45 

ELSEIF <M  .GT.  N)THEN 

PRINT*, 'WARNING  -  THERE  IS  A  ZERO  E-VECTOR' 

CALL  EXIT 
ENDIF 

DO  40  1*1, N 

EVEC  < I , J ) *EVEC ( I , J ) /ZN 
40  EVECA(I,J)*EVECA(I,J)/ZNA 

PRINT  CALCULATED  VALUES 

WRITE (6, 50) TITLE 
WRITE  <6,50) LABEL 

WRITE(6,55) <  <F(I, J)  ,  J*l,4)  ,  1*1,2) 

50  FORMAT </17X, A) 

55  FORMAT < //20X , ' FEEDBACX  GAIN  MATRIX  F  ' 
*//2U0X,4(F8.4,2X)/)//) 

WRITE (6, GO) 

DO  65  1*1, N 

65  WRITE (6, 70) (AAUG(I , J) , J*1 ,N) 

60  FORMAT <//20X, 'AUGMENTED  A  MATRIX'//) 

70  F0RMAT<2<6( 1X,F10.3)/) //) 

WRITE<6,75) 

%  75  FORMAT <//20X,' DUTCH  ROLL  MODE  -  REAL' 

*  /2X, 'EIGENVALUE', 


92 


*3X, 'DESIRED  EIGENVECTOR' ,5X, 'ACHIEVED  EIGENVECTOR'//) 
K*1 

IFdDRR(l)  .EQ.DTHEN 

WRITE(6,85)REAL(LAMDR) ,VDDRR(K) ,VDRR<1) 

K*K*1 

ELSE 

WRITE<6,90)REAL(LAHDR),VDRR(1) 

ENDIF 

DO  80  1*2, N 

IF<IDRR(I). EQ.DTHEN 

WRITE(6,95)VDDRR<K),VDRR<I) 

K*K»1 

ELSE 

WRITE(6,100)VDRR<I) 

ENDIF 

80  CONTINUE 

85  FORMAT<4X,F7.3,10X,F10.3,11X,F10.3) 

90  F0RHAT(4X,F7.3,15X,'X',15X,F10.3> 

95  F0RMAT(21X,F10.3,11X,F10.3) 

100  F0RMAT(26X, 'X' ,15X,F10.3) 

WRITE<6,110) 

110  F0RMAT(///2OXf 'DUTCH  ROLL  MODE  -  I MAG' 

*//2X, 'EIGENVALUE', 

*3X, 'DESIRED  EIGENVECTOR '5Xf 'ACHIEVED  EIGENVECTOR'//) 
K«1 

IF(IDRI (1)  .EQ.  1)THEN 

WRITE <6, 85 ) AIMAG ( LANDR) , VDDRI <  K ) , VDRI ( 1 ) 

K*K*1 

ELSE 

WRITE (6 , 90) AIMAG  <LAMDR) , VDRI < 1 ) 

ENDIF 

DO  120  1*2, N 

IF(IDRKI). EQ.DTHEN 

WRITE (6 , 95) VDDRI < X ) , VDRI < I ) 

K*K*1 

ELSE 

WRITE ( 6 , 100 ) VDRI ( I ) 

ENDIF 

120  CONTINUE 

WRITE(6,130) 

130  FORMAT (///24X, 'SPIRAL  M0DE'//2X, 'EIGENVALUE' ,3X, 
•'DESIRED  EIGEN VECTOR ',5X,' ACHIEVED  EIGENVECTOR'//) 

K-l 

IF(IDSP(1)  .EQ.  1 ) THEN 
WRITE (6,85) LAMSP , VDSP  <  K ) , VSP  <1 ) 

X*K*1 

ELSE 

WRITE (6, 90) LAMSP, VSP (1) 

ENDIF 

DO  140  1*2, N 

IF  (IDSP(I).  EQ.DTHEN 

WRITE <6, 95) VDSP (K) ,VSP( I) 


o  o  o  n  o  non 


WRITE(6,100)VSP(I) 

ENDIF 

140  CONTINUE 

WRITE <6, 150) 

ISO  FORMAT ( /  / / 24X , ' ROLL  MODE' //2X, 'EIGENVALUE' ,3X, 
•'DESIRED  EIGENVECTOR' ,5X,' ACHIEVED  EIGENVECTOR'//) 

K*1 

IF < IDRS ( 1 )  .EO.  1)THEN 

WRITE <6, 8S  >  LAMRS , VDRS  <  K ) , VRS  < 1 ) 

K*K*1 

ELSE 

WRITE (6,90) LAMRS , VRS  < 1 ) 

ENDIF 

DO  160  1*2, N 

IF ( IDRS ( I ) . EQ . 1 ) THEN 

WRITE (6,95) VDRS <K) ,VRS(I) 

K«K*1 

ELSE 

VRIT£(6,100)VRS<I) 

ENDIF 

160  CONTINUE 

WRITE (6, 170) 

DO  160  K-l.N 

WRITE (6 . 190 > EVALA < K ) , EVECA < 1 ,K ) 

DO  180  M-2.N 

180  WRITE (6,200) EVECA(H,K) 

WRITE <6, 210) 

DO  220  K-1,N 

WRITE(6,190)EVAL(K) ,EVEC(1,K) 

DO  220  H»2,N 

220  WRITE(6,200)EVEC(M,X) 

170  FORMAT (//1SX, 'OPEN  LOOP  SPECTRAL  DECOMPOSITION'// 

•10X , ' EIGENVALUES' , 15X , ' EIGENVECTORS' / ) 

210  FORMAT (//1SX, 'CLOSED  LOOP  SPECTRAL  DECOMPOSITION'// 
•10X, 'EIGENVALUES' , 15X, ' EIGENVECTORS' / ) 

190  F0RHAT(//2X,F10.3,3X,F10.3,'i',5X,F10.3,3X,F10.3,'i') 
200  F0RMAT(31X,F10.3,3X,F10.3,'i') 

STOP 

END 


SUBROUTINE  REVEC 

• (A,B,N,W, IDRE, VDRE,LAMRE, VRE, INUM, IRUN) 

THIS  SUBROUTINE  CALCULATES  THE  LEAST  SQUARES 
PROJECTION  OF  THE  DESIRED  EIGENVECTOR  COMPONENTS 
INTO  THE  ALLOWABLE  SUBSPACE.  AT  LEAST  TWO 
COMPONENTS  MUST  BE  SPECIFIED. 


£ 


ooo  o  o  o  oooortoo  non 


REAL  WK<550) , 

•A(20,20) ,B<20,2) , VDRE(20) , LAMRE, VR£( 20) , LAMDAI (20,20) , 
•LAMINV(20,20) ,L<20,2> ,LTILDA<20,2) ,L2TILDA<20,2) , 

•S(20) ,LTILIV<2,20) ,LPR0J(20,20) ,VPROJ<20) ,Z<20) ,W<2,5> 
INTEGER  N,IDGT,IDRE(20) , IER, INUN , IR0W< 20) , MROW, IRUN 
IDGT-0 

CALCULATE  (LANDA*I  -  A)  INVERSE 

DO  10  1*1, N 
DO  10  J*1,N 

IF(I.EQ.J)  THEN 

LAMDAI (I, I) * LAMRE- A (I, I) 

ELSE 

LAMDAI ( 1, J) *-A< I, J) 

ENDIF 

10  CONTINUE 

CALL  L1NV2F (LAMDAI, N, 20, LAMINV,IDGT,WK, IER) 

IF ( IER . EO . 129. AND . IRUN . EQ . 3) PRINT* , ' WARNING  -  THE  ' 

•, 'SPIRAL  LAMDA*I-A  MATRIX  IS  SINGULAR' 

IF ( IER. EQ. 129. AND. IRUN. EQ. 4)  PRINT*, 'WARNING  -  THE  ' 

*, 'ROLL  MODE  LAMDA*I-A  MATRIX  IS  SINGULAR' 

CALCULATE  LAMDA*1-A  INVERSE  X  B 

CALL  VMULFF(LAMINV, B,N,N, 2,20,20, L, 20, IER) 

CALCULATE  ROWS  OF  LAMDA*1-A  INVERSE  X  B  TO  BE  USED 
IN  CALCULATING  LEAST  SQUARES  PROJECTION 

MROW-1 
DO  20  1*1, N 

IF(IDREU)  .EQ.  1)  THEN 
IROW(MROW)*I 
MR0W*MR0W+1 
ENDIF 

20  CONTINUE 

DO  30  I*1,INUM 
DO  30  J*l,2 

LTILDA(I,J)*L(IROW(I) ,J) 

30  L2TILDA ( I , J ) *LTILDA ( I , J ) 

CALCULATE  GENERALIZED  INVERSE 

CALL  LGINF (LTILDA,20, INUM,2,0,LTILIV,2,S,WK, IER) 

IF ( IER. EQ. 33) PRINT*, 'WARNING  -  GENERALIZED  INVERSE  IS  ' 
», 'ILL  CONDITIONED' 

IF<IER.EQ.29)PRINT*, 'WARNING  -  GENERALIZED  INVERSE' 

*,'  NOT  OBTAINED  -  COMPUTATION  TERMINATED' 

FIND  PROJECTION  OF  DESIRED  VECTOR 


oooo  ooo  oooooo  ooy 


CALL  VMULFF (L2TILDA , LTILI V , INUH , 2 , I NUN , 20 , 2 , LPROJ , 

*20, IER) 

CALL  VMULFF ( LPROJ , VDRE , INUM , IMUM , 1 , 20 , 20 , VPROJ , 20 , 1ER ) 

CALCULATE  LINEAR  COMBINATION  OF  SUBSPACE  AND 
REMAINING  COMPONENTS  OF  VECTOR 

CALL  VMULFF (LTILI V , VPROJ ,2 , INUM , 1 , 2 , 20 , Z , 20 , IER ) 

CALL  VHULFF(L,Z,N,2,1,20,20,VRE,20,IER> 

DO  40  1*1, INUH 
40  VRE(IROW(I))*VPROJ(I) 

DO  SO  1*1,2 
50  W(I,IRUN)*Z(I) 

RETURN 

END 


SUBROUTINE  DRVEC ( A , B , N , V , IDRR , I DRI , VDDRR , VDDRI , LAMDR , 
• VDRR , VDRI , IRNUM , I I NUM ) 

THIS  SUBROUTINE  COMPUTES  THE  LEAST  SQUARES  PROJECTION 
OF  THE  COMPLEX  DUTCH  ROLL  EIGENVECTOR.  AT  LEAST  TWO 
COMPONENTS  OF  EACH  OF  THE  REAL  AND  IMAG  PARTS  MUST  BE 
SPECIFIED. 

REAL  A(20,20) ,B(20,2) ,VDRR<20> , VDRI (20) ,VDDRR(20) , 
•VDDRK20)  ,WK(150)  ,LTILDA (40,4)  ,L2TILDA (40,4)  , 
•VDPRIM(40) ,S(8,0) ,LPRIN(40,40) , 

"LTILIV(4,40> , LPROJ ( 40, 40 ) , VPROJ <40 ) ,Z(40) , 
•VDPRIM(40),W(2,S) 

INTEGER  N, 

• IDRR < 20) , IDRI (20) , INUM , IRNUH , I INUM , MROW , I ROW ( 20 ) 
COMPLEX  LAMDR, LAMDAI (20,20) ,LAMINV(20,20) ,WA(550) , 
•L(20,20) 

CALCULATE  LAHDA*I-A  INVERSE 

DO  10  1*1, N 
DO  10  J-l,l 

IF(I.EQ.J)  THEN 
LAMDAI (I, I) -LAMDR- A (I, I) 

LAMINV(I,I)*(1.0,0. ) 

ELSE 

LAHDAI(I,J)*-A(I,J) 

LANINVd,  J)-(0.  ,0. ) 

ENDIF 

10  CONTINUE 

CALL  LEQ2C( LAMDAI ,N,20,LAMINV,N,20,0,WA, WK, IER) 
IF(IER.EQ.129)PRINT*, 'WARNING  -  LAMDA*I-A  FOR  THE' 

*,'  DUTCH  ROLL  MODE  IS  SINGULAR' 

IF(IER.EQ.130)PRINT*, 'WARNING  -  LAMDA*I-A  FOR  THE' 


o  o  ft  oooo  o  n  o  ooo 


.•nor*  xv  xv*v  xv  xv  vv  yv  vwv  v_  wxv  ■r_v-.'-,,'1.'vv'v.v'.^r.’v.’vyw^’V’.’v^'^._. 


*,'  DUTCH  ROLL  MODE  IS  ILL  CONDITIONED' 

CALCULATE  LAMDA*I-A  INVERSE  X  B 

DO  20  1*1, N 
DO  20  J*l,2 
Ld,J>«<0.,0.) 

DO  20  M«1,N 

20  Ld,J)*Ld,J>*LAMINVd,M>*B(H,J> 

CALCULATE  ROWS  TO  BE  USED  IN  FINDING  LEAST  SQUARES  FIT 

MR0W*1 
DO  30  1*1, N 

IF<IDRR(I) .  EQ.  1)  THEN 
IROW(MROW)*I 
BROV-MROW-1 
ENDIF 

30  CONTINUE 
DO  40  1*1, N 

IF(IDRId)  .EQ.  1)  THEN 
IROW<HROW)*I*N 
HROW*MROWd 
ENDIF 

40  CONTINUE 

INUM-IRNUMdINUH 

FORM  REAL  MATRIX  REPRESENTATION  TO  FIND  GENERALIZED 
INVERSE 

DO  50  1*1, N 
DO  50  J*l,2 

LPRIMd,  J)=REAL(Ld, J) ) 

LPRIM(I, J+2)=-AIMAG(L(I, J) ) 
LPRIHd*N,J)*AIHAG<Ld,J>> 

50  LPRIMd«-N,J+2>*LPRIHd,J> 

DO  55  I*1,1NUM 
DO  55  J-1,4 

LTILDA(I,J)*LPRIM(IROW(I) ,J> 

55  L2TILDA (I , J ) *LTILDA (I , J ) 

DO  60  I*1,IRNUM 
60  VDPRIM (I ) * VDDRR <1 ) 

DO  70  1*1,1INUH 
70  VDPRIMd*IRNUM)«VDDRId) 

CALCULATE  GENERALIZED  INVERSE  OF  AUGMENTED  MATRIX 

CALL  LGINF <LTILDA,40,INUH,4,0,LTILIV,4,S,WK,IER) 
IF(IER.EQ.33)PRINT», 'WARNING  -  GENERALIZED  INVERSE  OF' 
*, 'DUTCH  ROLL  AUGMENTED  MATRIX  IS  ILL  CONDITIONED' 

IF (IER.EQ. 29) PRINT*, 'WARNING  -  GENERALIZED  INVERSE  OF' 
«, 'DUTCH  ROLL  AUGMENTED  MATRIX  NOT  OBTAINED  - 


O  O  O  OOOOOOO  OOOO  OOO 


COMPUTATION  TERMINATED' 

CALCULATE  PROJECTION  OF  DESIRED  VECTOR 

CALL  VHULFF (L2TILDA , LTILIV , INUM , 4 , INUM , 40 , 4 , LPROJ , 
•40, IER) 

CALL  VMULFF 

*< LPROJ , VDPRIM , INUM , INUM , 1 , 40 , 40 , VPROJ , 40 , IER ) 

FIND  LINEAR  COMBINATION  OF  THE  SUBSPACE  AND  CALCULATE 
THE  UNSPECIFIED  E-VECTOR  COMPONENTS 

CALL  VMULFF (LTILIV, VPROJ ,4, INUM, lf 4, 40 ,2, 40, IER) 

CALL  VMULFFCLPRIH, Z,2«N,4,1, 40, 40, VDPRIM, 40, IER) 

DO  75  1*1, INUM 
75  VDPRIM ( IRON < I ) ) *VPROJ  < I ) 

DO  80  1*1, N 

VDRR<D<VDPRIM<D 
80  VDRI<D*VDPRIH<N*I) 

DO  90  1*1,2 
V<I,1)*Z<I) 

90  V<I,2)*Z<I+2) 

RETURN 

END 


SUBROUTINE  FCALC<C,W,N,V,F) 

THIS  SUBROUTINE  CALCULATES  THE  FEEDBACK  MATRIX,  F 
TO  PLACE  THE  EIGENVALUES  AND  EIGENVECTORS 

REAL  C<5,20) ,W<2,5) ,V<20,5) ,F<2,5> ,WK<350) , 
*CV<5,20),CVINV<20,20) 

INTEGER  N,IDGT 

CALCULATE  FEEDBACK  MATRIX  F 

CALL  VMULFF <C,V,4,N,4,5,20,CV,5,IER> 

IDGT*0 

CALL  LIMV2F <CV,4,5,CVINV,IDGT,WK,IER) 

IF<IER.EQ.129  .OR.  IER.EQ.13DPRINT*, 'WARNING  -  ', 
•'THE  CV  MATRIX  IS  EITHER  SINGULAR  OR  VERY' 

•'  ILL-CONDITIONED' 

CALL  VMULFF <tf ,CVINV,2,4,4,2,5,F,2, IER) 

RETURN 


Listed  below  are  aesple  outputs  of  the  eigenatructure  assignment 
program.  These  examples  are  from  the  12th  order  models  of  the  YA-7D 


described  in  the  text.  The  desired  eigenstructures  are  also  derived  in 
the  text  with,  <0/0>d  of  zero  in  these  cases.  Also  included  as  part  of 
the  output  are  spectral  decompositions  of  the  open  loop  and  closed  loop 
systems.  It  is  easily  verified  that  the  desired  eigenvalues  and 
eigenvectors  were  exactly  achieved. 


1 


YA-7D  at  0.4  Mach 
phi /beta  *  0 


FEEDBACK  GAIN  MATRIX  F 

0.4596  0.0401  0.1066  0.0149 
1.1719  0.0796  -0.0635  -0.0063 


I 

A 

V 

I 


DUTCH  ROLL  MODE  -  REAL 

EIGENVALUE  DESIRED  EIGENVECTOR  ACHIEVED  EIGENVECTOR 


-1.200  X  3.851 

1.000  1.000 

X  -0.000 

0.  0.000 

X  -0.610 

X  0.299 

X  -0.344 

X  4.953 

X  0.622 

X  9.928 

X  4.570 


99 


u 


a 


Of 


I  Mi  • 


aw mwi 


-57.349 


DUTCH  ROLL  NODE  - 


EIGENVALUE  DESIRED  EIGENVECTOR 


ACHIEVED  EIGENVECTOR 


2.750 


X 

1.000 

X 

0. 

X 

X 

X 

X 

X 

X 

X 

X 


-1.805 

1.000 

0.000 

0.000 

-1.667 

-4.052 

-1.651 

1.036 

-3.882 

6.370 

-0.456 

-36.299 


8 

' 


SPIRAL  NODE 


EIGENVALUE  DESIRED  EIGENVECTOR 


ACHIEVED  EIGENVECTOR 


-0.025 


X 

0. 

X 

1.000 

X 

X 

X 

X 

X 

X 

X 

X 


0.076 

-0.000 

-0.025 

1.000 

0.007 

-0.011 

0.007 

-0.000 

-0.011 

0.000 

-0.002 

-0.086 


ROLL  NODE 


EIGENVALUE  DESIRED  EIGENVECTOR 


ACHIEVED  EIGENVECTOR 


-4.000 


0. 

0. 

1.000 


-0.001 

0.003 

1.000 


X 

X 

X 

X 

X 

X 

X 

X 

X 


-0.250 

0.224 

0.153 

0.179 

-0.717 

0.135 

-0.539 

0.163 

-0.236 


OPEN  LOOP  SPECTRAL  DECOMPOSITION 


EIGENVALUES 


-1.000  0.  i 


10.000  0.  i 


-1.311  0.  i 


EIGENVECTORS 


0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

1.000  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

0.  0. 

1.000  0. 

1.000  0. 

0.266  0. 

-11.566  0. 

6.624  0. 

-0.000  0. 

-0.000  0. 

0.  0. 


i 

i 

i 

i 

1 

1 

1 

1 

1 

i 

1 

i 


1 

i 

i 

i 

i 

i 

i 

1 

i 

i 

i 

1 

1 

i 

1 

i 

1 

i 

i 
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-0.037 


0.  i 


1.000 

0.140 

-0.494 

13.273 

-0.000 

-0.000 

0. 

0. 

0. 

0. 

-0.036 

-6.275 


0.  1 
0.  i 
0.  i 
0.  i 
0.  i 
0.  i 
0.  i 
0.  1 
0.  i 
0.  i 
0.  i 
0.  i 


-0.354  1.7211  1.000  0.  1 

0.157  0.5801 

-2.853  -1.3871 

-0.446  1.7501 

0.000  -0.0001 

0.000  -0.0001 

0.  0.  1 

0.  0.  1 

0.  0.  1 

0.  0.  1 

0.606  0.1891 

-11.541  -22.7871 


-0.354  -1.7211  1.000  0.  1 

0.157  -0.5801 

-2.853  1.3871 

-0.446  -1.7501 

0.000  0.0001 

0.000  O.OOOi 

0.  0.  i 

0.  0.  1 

0.  0.  1 

0.  0.  1 

0.606  -0.1891 

-11.541  22.7871 


33.330 


0.  1 


1.000 

0.018 

-0.785 


0.  1 
0.  i 
0.  1 
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-20.000 


0.  i 


-65.950 


67.3321 


-65.950 


-67.3321 


0.024 

0.  i 

-0.000 

0.  1 

12.923 

0.  1 

0. 

0.  1 

0. 

0.  1 

0. 

0.  1 

0. 

0.  1 

0.924 

0.  1 

-15.419 

0.  i 

1.000 

0.  1 

0.054 

0.  i 

-11.887 

0.  1 

0.594 

0.  1 

-17.947 

0.  1 

-0.000 

0.  1 

0. 

0.  1 

0. 

0.'  1 

0. 

0.  1 

0. 

0.  1 

-0.491 

0.  i 

-41.611 

0.  1 

1.000 

0.  1 

-0.005 

0.0081 

-0.753 

-0.0121 

0.006 

0.0051 

-0.000 

O.OOOi 

25.733 

-26.4891 

0. 

0.  i 

0. 

0.  1 

28.328 

77.9091 

-7114.024 

-3230.7461 

0.907 

0.0051 

-12.124 

-0.8901 

1.000 

0.  1 

-0.005 

-0.0081 

-0.753 

0.0121 

0.006 

-0.0051 

-0.000 

-O.OOOi 

25.733 

26.4891 

0. 

0.  1 

0. 

0.  1 

28.328 

-77.9091 

-7114.024 

3230.7461 

0.907 

-0.005i 

-12.124 

0.8901 

-65.950  67.3321 


-65.950  -67.3321 


1.000 

-0.005 

-0.753 

0.006 

0.000 

25.733 

0.001 

-0.173 

28.328 

-7114.066 

0.907 

-12.124 


1.000 

-0.005 

-0.753 

0.006 

0.000 

25.733 

0.001 

-0.173 

28.328 

-7114.066 

0.907 

-12.124 


0.  1 
0.0081 
-0.0121 
0.0051 
-O.OOOi 
-26.4891 
0.0021 
-0.0781 
77.9101 
-3230.7581 
0.0051 
-0.890i 


0.  1 
-0.008i 
0.0121 
-0.005i 
O.OOOi 
26.4891 
-0.002i 
0.0781 
-77.9101 
3230.7581 
-0.0051 
0.8901 


CLOSED  LOOP  SPECTRAL  DECOMPOSITION 


EIGENVALUES 


EIGENVECTORS 


-67.371 


67 . 963i 


1.000 

-0.005 

-0.900 

0.005 

-2.160 

25.437 

6.522 

-6.964 

28.897 

-7339.977 

0.889 

-12.291 


1.000 

-0.005 

-0.900 

0.005 


104 


0.  i 
0.0081 
-0.235i 
0.0081 
-0.4131 
-26.9121 
-6.3621 
871.8801 
79.3551 
-3382. 228i 
-0.0221 
-1.1641 


0.  1 
-0.0081 
0.2351 
-0.008i 


.  f.  c 


.  «*.  >•.  •  .  A  ■ 


\v-w\v\v.v; 


-67.371 


-67.9631 


65.050 


66.6121 


65.050 


-66.6121 


33.747 


0. 


-2.160 

0.4131 

25.437 

26.9121 

6.522 

6.3621 

-6.964 

-871.8801 

28.897 

-79.3551 

-7339.977 

3382.2281 

0.889 

0.0221 

-12.291 

1.1641 

1.000 

0.  1 

-0.005 

0.0071 

-0.595 

0.1711 

0.006 

0.003i 

1.938 

0.0971 

26.148 

-26.1561 

-4.689 

6.2371 

-110.424 

-718.0681 

27.390 

77.1501 

-6920.886 

-3194. 049i 

0.926 

0.0281 

-11.942 

-0.6601 

1.000 

0.  1 

-0.005 

-0.0071 

-0.595 

-0.1711 

0.006 

-0.0031 

1.938 

-0.0971 

26.148 

26.1561 

-4.689 

-6.2371 

-110.424 

718.0681 

27.390 

-77.1501 

-6920.886 

3194.0491 

0.926 

-0.0281 

-11.942 

0.6601 

1.000 

0.  1 

0.017 

0.  1 

-0.580 

0.  1 

0.017 

0.  1 

0.574 

0.  1 

13.324 

0.  1 

-0.395 

0.  1 

13.324 

0.  1 

-0.167 

0.  1 

5.627 

0.  i 

0.949 

0.  1 

-15.043 

0.  1 
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% 


I*" 


.ij.  >jwt i  u 


ryyv>vv 'jvwvv  .v,'.v.v.v.\".v.\\v  vtv v vv v ^ v.- v tt v. -c, v.* 


i\ 


-14.597  0.  i 


-9.248  0.  i 


-1.200  2.7501 


-1.200  -2.7501 


1.000 

0. 

1 

0.019 

0. 

i 

28.103 

0. 

i 

-1.925 

0. 

i 

33.247 

0. 

1 

20.385 

0. 

i 

8.981 

0. 

1 

-131.105 

0. 

i 

11.457 

0. 

1 

-167.253 

0. 

1 

4.769 

0. 

1 

52.053 

0. 

1 

1.000 

0. 

1 

0.091 

0. 

1 

5.993 

0. 

1 

-0.648 

0. 

i 

4.666 

0. 

1 

5.794 

0. 

1 

2.508 

0. 

1 

-23.197 

0. 

1 

4.187 

0. 

i 

-38.722 

0. 

1 

1.938 

0. 

i 

16.990 

0. 

i 

1.000 

0.  1 

0.113 

0.3131 

-0.000 

-O.OOOi 

-0.000 

0.0001 

0.037 

-0.4161 

0.468 

-0.8331 

0.092 

-0.3861 

0.951 

0.7151 

0.520 

-0.7641 

1.476 

2.3461 

1.019 

0.3591 

-8.588 

-13.4521 

1.000 

0.  1 

0.113 

-0.3131 

-0.000 

0.0001 

-0.000 

-O.OOOi 

0.037 

0.416i 

0.468 

0.8331 

0.092 

0.386i 

0.951 

-0.715i 

0.520 

0.764 

1.476 

-2.346. 
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vi 


a 
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» 


► 


a 
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YA-7D  at  0.6  Mach 


phi/b«ta  *  0 


FEEDBACK  GAIN  MATRIX  F 


0.3091 

0.5725 


0.0155 

0.0160 


0.0465 

-0.0033 


0.0057 

-0.0032 


DUTCH  ROLL  MODE 


EIGENVALUE  DESIRED  EIGENVECTOR 


ACHIEVED  EIGENVECTOR 


-1.200 


X 

1.000 

X 

0. 

X 

X 

X 

X 

X 

X 

X 

X 


3.619 

1.000 

-0.000 

-0.000 

-0.582 

0.429 

-0.399 

3.491 

0.536 

3.261 

4.560 

-112.282 


DUTCH  ROLL  MODE  - 


EIGENVALUE  DESIRED  EIGENVECTOR 


ACHIEVED  EIGENVECTOR 


2.750 


X 

1.000 

X 

0. 

X 

X 

X 

X 

X 

X 

X 

X 


-1.779 

1.000 

-0.000 

0.000 

-1.080 

-1.510 

-1.095 

0.218 

-1.421 

3.185 

-0.445 

-61.930 
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SPIRAL  NODE 

EIGENVALUE 

DESIRED  EIGENVECTOR 

ACHIEVED  EIGENVECTOR 

-0.025 

X 

0.051 

0. 

-0.000 

X 

-0.025 

1.000 

1.000 

X 

0.003 

X 

-0.005 

X 

0.003 

X 

-0.000 

X 

-0.005 

X 

0.000 

X 

-0.001 

X 

-0.077 

ROLL  NODE 

EIGENVALUE 

DESIRED  EIGENVECTOR 

ACHIEVED  EIGENVECTOR 

-4.000 

0. 

-0.001 

>> 

0. 

0.003 

1.000 

1.000 

X 

-0.250 

X 

0.063 

X 

0.010 

X 

0.050 

X 

-0.201 

X 

0.009 

X 

-0.036 

X 

0.049 

X 

-0.935 

OPEN  LOOP  SPECTRAL 

DECOMPOSITION 

EIGENVALUES 

EIGENVECTORS 

-1.000 

0.  1 

0.  0. 

i 

0.  0. 

i 

0.  0. 

i 

0.  0. 

i 

0.  0. 

1 

0.  0. 

i 

❖ 

0.  0. 

1 
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v\-  /  ■>  >>  ■■> 

.  «■.  ■f.  -r. 

R 

R 


-0.310 


-2.0391 


20.000 


0.  1 


33.330 


-0.915 

1.4881 

0.000 

0.0001 

-0.000 

-0.0001 

0.000 

0.0001 

-0.000 

-0.0001 

0.000 

0.0001 

0.000 

0.0001 

0.802 

0.3211 

-15.386 

-40.0251 

1.000 

0.  1 

0.076 

-0.507i 

-2.752 

2.3281 

-0.915 

-1.488i 

0.000 

-O.OOOi 

-0.000 

0.0001 

0.000 

-O.OOOi 

-0.000 

O.OOOi 

0.000 

-0.0001 

0.000 

-O.OOOi 

0.802 

-0.3211 

-15.386 

40.0251 

1.000 

0.  1 

-0.132 

0.  1 

682.003 

0.  1 

-34.095 

0.  1 

454.378 

0.  1 

-0.000 

0.  i 

0.000 

0.  1 

-0.000 

0.  1 

0.000 

0.  1 

0.000 

0.  1 

27.738 

0.  1 

1946.593 

0.  1 

1.000 

0.  1 

0.022 

0.  1 

-1.021 

0.  i 

0.031 

0.  i 

0.000 

0.  1 

6.724 

0.  i 

0.000 

0.  i 

-0.000 

0.  i 

0.000 

0.  1 

0.000 

0.  1 

0.991 

0.  i 

-14.186 

0.  1 

65.950 


67.3321 


65.950 


-67.3321 


65.950 


67.3321 


65.950 


-67.3321 


-0.113 

0.0681 

-0.001 

-O.OOli 

1.000 

0.  i 

-0.007 

-0.0081 

2.334 

-2.0431 

-0.346 

2.7931 

1.516 

12.5521 

-945.130 

-725.7591 

-5.304 

-3.4321 

560.662 

-130.7771 

-0.077 

0.0881 

2.664 

-0.6921 

-0.113 

-0.0881 

-0.001 

O.OOli 

1.000 

0.  i 

-0.007 

0.0081 

2.334 

2.043i 

-0.346 

-2.7931 

1.516 

-12.5521 

-945.130 

725.7591 

-5.304 

3.4321 

560.662 

130.7771 

-0.077 

-0.0681 

2.664 

0.6921 

1.000 

0.  1 

-0.001 

0.0091 

-0.463 

-2.9991 

-0.019 

0.0261 

-6.362 

-8.3601 

13.531 

-13.5991 

42.760 

-2.2121 

-2671.058 

3024. 966i 

14.229 

40.643i 

-3674.957 

-1722.3301 

0.989 

-0.105i 

-10.419 

-5.7481 

1.000 

0.  1 

-0.001 

-0.0091 

-0.483 

2.9991 

-0.019 

-0.0261 

-6.362 

8.3601 

13.531 

13.5991 

42.760 

2.2121 

-2671.056 

-3024.966i 

14.229 

-40.6431 

-3674.957 

1722.3301 
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.*■  .n  , 


0.969 

-10.419 


0.1051 

5.7481 


4 


CLOSED  LOOP  SPECTRAL  DECOMPOSITION 


EIGENVALUES 


-66.236 


67.2511 


-66.236 


-67.2511 


-66.564 
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This  appendix  contains  frequency  reaponae  plota  of  0/Fa  and  Ycf  for 
the  unaugaented  and  augmented  YA-7D  at  0.4  and  0.6  Mach.  The  responses 
were  generated  froa  the  14th  order  eigenatructurea  described  in  the 
text.  The  unaugaented  YA-7D  includes  only  the  mechanical  flight  control 
systea,  while  the  augmented  eigenatructurea  contain  effects  froa  both 
the  mechanical  ayatea  and  roll  CAS.  The  only  difference  in  the 
augmented  eigenatructurea  ia  the  design  <0/0)<j- 
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Figure  Cl.  0/Fa  Frequency  Reaponae  for  Unaugaented  YA-7D  -  0.4M 
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Figure  C3.  0/Fa  Frequency  Response  for  (0/0>d  =1.5  arg<0)  -  0. 
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Figure  C14.  0/Fa  Frequency  Reaponae  for  (0/0>d  *  3.0  arg(0)  -  0. 
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Figure  C15.  0/Fa  Frequency  Reaponae  for  <0/0>d  *  3.0  arg<60)  -  0 
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APPENDIX  D.  Siaulator  Tlae  Reaponae  Plot* 
for  Selected  Elqenatructurea 

Tlae  reaponae  plota  were  generated  ualng  the  Siaulator  for  Aircraft 
Flight  Teat  and  Developaent  (SAFTD).  Both  free  and  forced  reaponaea  at 
the  aoat  aevere  flight  condition,  0.4  Mach,  are  included  for  the 
unaugaented  YA-7D,  (0/0)<i  of  zero,  and  l#/0ld  of  1.5,  and  3.0,  at  phaae 
anglea  of  0,  60,  and  120  degreea.  At  0.6  Mach,  only  the  free  and  forced 
reaponae  of  the  unaugaented  YA-7D  and  <0/0)d  of  zero  are  ahown.  The 
variablea  are  in  body  axea,  and  include  effecta  of  aenaor  dynaaica  and 
the  digital  flight  control  ayatea. 

Inputa  were  applied  by  teat  taaa  pilota.  The  free  reaponae  waa  a 
releaae  froa  a  winga  level  aidealip  of  approxiaately  three  degreea.  The 
forced  reaponae  waa  a  atep  lateral  atick  input. 
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>.*.  Figure  Dl.  Unaugaantad  YA-7D  Fraa  Raaponsa  -  0.4N 
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Figura  D2.  Unaugaantad  YA-7D  For cad  Raaponaa  -  0.4N 
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Figura  D2  (CONTINUED).  Unaugaantad  YA-7D  Forcad  Raaponaa  -  0.4N 
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Figure  D3  (CONTINUED).  Fr««  Response  (4/0)4  *  0  -  0.4N 
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Figure  D3.  Free  Response  (0/0)<j  “1.5  arg(O) 
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Figura  D5  (CONTINUED).  Fraa  Response  (#/0)d  *  1.5  arg(O)  -  0.4H 
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Figure  D6  (CONTINUED).  Forced  Response  (0/0)d  «  1.5  arg(O) 
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Figure  D7  (CONTINUED).  Free  Response  <0/p>d  *  1*5  arg<60) 
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Figure  D8.  Forced  Response  (0/0) d  *  1.5  arg(SO) 
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Figure  D9  (CONTINUED).  Free  Response  (4/0)4  »  1.5  arg(120) 
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Figure  D12.  Forced  Reeponse  (#/0)<i  *  3.0  arg(O)  -  0.4N 
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Figure  D14  (CONTINUED).  Forced  Response  (0/0)d  =  3.0  arg(60)  -  0.4N 
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Figure  D15.  Free  Response  <*/p>d  *  3.0  arg(120)  -  0.4N 
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Figure  D15  (CONTINUED).  Free  Response  (0/0>d  *  3.0  arg(120)  -  0 
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Figure  D16.  Forced  Response  (#/0)<j  3  3.0  arg(120)  -  0.4N 
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Figure  D18  (CONTINUED). 


Unaugeented  YA-7D  Forced  Response  -  0.6H 
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Figure  19.  Free  Response  <0/0)4  =  0  -  0.6N 
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APPENDIX  E.  Flight  Test  Data 


Figuras  El,  E2,  and  E3  show  intar  and  intra  pilot  rapaatability 
for  HUD  Tracking  task  1,  HUD  tracking  task  2,  and  tha  air-to-air 
tracking  task,  raspactivaly.  Tha  intar  pilot  coaparisons  ara  a  plot  of 
tha  pilot's  highast  rating  varaua  lowast  rating  for  aach  control  law. 

Tha  dashad  lina  raprasants  a  diffaranca  of  one  Cooper-Harper  rating. 

Tha  intra  pilot  coaparisons  ara  a  plot  of  one  pilot's  rating  versus  a 
second  pilot's  rating  for  each  control  law.  The  method  used  to 
dateraine  a  pilot's  rating  for  aach  control  law  is  discussed  in  Appendix 
F.  The  dashad  lina  raprasants  a  difference  of  plus  or  sinus  one 
Cooper-Harper  rating. 
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Table  El 


HUD  Tracking  Task 

1  Cooper-Harper  Ratings 

CONFIGURATION 

PILOT  A 

PILOT  B 

PILOT  C 

Unaugsented 

YA-7D 

5,5 

N/A 

5,5 
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6,6,6 

4, 5,5, 4 

Phi /Beta 

Nag 

Phi/Beta 

Phase 
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7,6 

7,7,6 
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1.5 
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5,6 
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Figura  E2.  Intar  and  intra  pilot  Coopar-Harpar  Rating 
Comparisons  for  HUD  Tracking  Task  2 
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Table  E2 


HUD  Tracking  Task  2  Cooper-Harper  Ratings 


CONFIGURATION 


Unaugsented  YA-7D 


Standard  YA-7D 


Phi /Beta 
Mag 


Phi/Beta 

Phase 


PILOT  A 


4. 5,4. 5 


5,4.5 


4, 5, 4. 5 


5,4.5 


5,5,4 
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4,4,4 


3, 5, 4, 4 


5,5,5 

2.3.3 

3.3.3 


4,5,5 


PILOT  C 


4, 4, 4, 4 


4, 4, 3, 4, 5. 5, 5 

3. 4. 3. 4. 4. 4. 3 

5. 4. 4. 5. 5.3. 4 
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4. 4. 4. 5. 3. 5. 3 

3. 5. 4. 5. 5. 3. 4 
4, 5, 4, 4, 4, 4, 3 


PILOT  RATING  PLT  A  LOWEST  PILOT  RATING 

Figura  E3.  Intar  and  intra  pilot  Coopar-Harpar  Rating 
Coapariaona  for  air-to-air  tracking  task 
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Table  E3 


Air-to-air 

Tracking  Task  Cooper-Harper  Ratings 

CONFIGURATION 

PILOT  A 

PILOT  B 

PILOT  C 

Unaugeented 

YA-7D 

4. 5. 5. 5 

7 

N/A 

Standard  YA- 

-7D 

4,4, 4,4. 5 

4, 4, 4, 4 

4, 4, 4, 4 

Phi/Beta 

Nag 

Phi/Beta 

Phase 

0 

0 

5, 4. 5, 4. 5, 4 

4,4,4 

5, 6, 6, 5 

1.5 

0 

4,5,5 

4,4,4 

6,5,5, 4 

1.5 

60 

3, 3, 5, 5, 4 

3*3, 4, 4 

5, 5, 4, 4 

1.5 

120 

4. 5, 4. 5, 4 

5,4.5 

7, 6, 6, 5 

3.0 
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3, 4, 5, 4. 5 

3,3,4 

4, 5, 4, 3 
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5,3,4 

4. 4, 5, 3 
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120 

4. 5, 5, 4, 5 
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7, 6, 6, 6 

jwwwray.’yvya  *vyiryvTKJfv*x**J'»^K^x*XA  v~  y’jcwvwrj'auvvirjrgv  w  vry  v»ir*  v  v  v  -  *  v 


APPENDIX  F.  Flight  Tact  Data  Acquisition 
and  Analysis 

Instrumentation  and  Data  Raduction 

The  Ampex  AR-700  magnetic  tape  recorder  onboard  the  aircraft 
recorded  all  126  DIGITAC  parameters  for  later  playback  and  review.  The 
available  parameters  were  listed  in  the  Partial  Flight  Manual  <15).  The 
primary  parameters  required  were  lateral  stick  force  and  position,  roll 
and  yaw  rate,  lateral  acceleration,  and  roll  and  sideslip  angles.  The 
data  were  played  back  after  a  mission  on  two  eight  channel  stripchart 
recorders  in  the  USAF  Test  Pilot  School  Telemetry  station.  The  TN 
station  was  also  used  for  realtime  reception  and  display  instrumentation 
data.  The  specific  parameter  /  channel  assignments  were  selected  by  the 
test  engineer  for  each  mission  based  on  data  requirements. 

Open  loop  testing  was  conducted  to  confirm  the  pole  locations  of 
the  lateral-directional  axes.  TN  data  were  hand  reduced  from  stripchart 
output  to  determine  roll  node  time  constant,  dutch  roll  frequency  and 
damping,  and  phi/beta  magnitude  and  phase  angle.  Figures  FI  and  F2  on 
the  following  pages  show  how  these  system  parameters  were  determined 
from  time  history  data. 

Closed  loop  testing  of  each  control  law  configuration  was  conducted 
to  assign  Cooper-Harper  ratings,  as  described  in  The  Use  of  Pilot  Rating 
in  the  Evaluation  of  Aircraft  Handling  Qualtlea  <13).  The  first  test 
point  of  each  mission  evaluated  standard  YA-7D  lateral-directional 
control  laws.  This  procedure  standardized  miasion-to-mission  and 
pilot-to-pilot  performance.  Camera  film  of  the  closed  loop  Head-up 
Display  was  reviewed  to  standardize  pilot  aggressiveness  and  ensure 
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Figura  F2.  Determination  of  Open  Loop  Dutch  Roll  Paraaetera 


pilot  performance  supported  the  Cooper-Harper  rating  given  in  flight 


Qualitative  data  were  collected  using  the  pilot  coaaent  shown  card  in 
Appendix  G ,  and  ratings  influenced  by  external  factors  eliminated.  The 
raw  data  for  all  Cooper-Harper  ratings  assigned  during  testing  are  in 
Tables  El,  E2,  and  E3  for  HUD  tracking  task  1,  HUD  tracking  task  2,  and 
the  air-to-air  tracking  task,  respectively. 

For  each  pilot/task/control  law  combination,  a  single  Cooper-Harper 
rating  was  determined  by  the  test  engineer.  In  most  cases  this  rating 
was  the  mode;  however,  in  a  few  cases  the  number  of  ratings  or 
distribution  warranted  a  rating  based  on  pilot  comments  and  engineer 
judgment. 

Cooper-Harper  ratings  were  compared  between  pilots  to  ensure 
consistent  evaluations.  Figures  El,  E2,  and  E3  show  inter  and  intra 
pilot  Cooper-Harper  rating  comparisons  for  HUD  Tracking  task  1,  HUD 
tracking  task  2,  and  the  air-to-air  tracking  task,  respectively. 
Cooper-Harper  ratings  were  correlated  by  task  for  each  pilot,  as  well  as 
a  three  rating  pilot  average,  with  the  different  control  law 
configurations.  Although  averaging  Cooper-Harper  ratings  does  not  have 
a  strict  mathematical  basis,  it  was  used  here  to  show  trends  and  make 
presentation  of  data  more  meaningful.  The  averaged  ratings  only  show 
relative  pilot  preferences,  and  no  longer  correspond  to  the  criteria  on 
the  Cooper-Harper  scale.  Each  control  law  represented  various 
combinations  of  design  magnitude  ratio  and  phase  angle  between  roll 
angle  and  sideslip  angle.  The  deaign  phi/beta  ratios  were  used  rather 
than  actual  since  the  actual  were  generally  within  IS  percent  of  design 
and  thia  allowed  more  meaningful  presentation  of  data. 
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Teat  Methods  and  Conditions 

The  phi/beta  teat  satrix  was  completed  at  15,000  feet  He  and  0.6 
Hach  with  Automatic  Maneuvering  Flaps  (AMF)  off.  Mach  0.6  corresponded 
to  305  Knots  Calibrated  Airspeed  (KCAS)  and  9.72  units  (4.63  degrees) 
angle  of  attack  at  27,000  pounds  gross  weight.  Test  points  were 
accomplished  within  ♦  5  Knots  Indicated  Airspeed  (K1AS)  and  ♦  1000  feet 
He  of  the  specified  conditions.  The  test  matrix  was  not  completed  at 
0.4  Mach  due  to  the  technical  problems  the  software. 

The  Simulator  for  Aircraft  Flight  Test  and  Development  verified  the 
flight  control  gains  produced  the  desired  aircraft  response.  This  also 
evaluated  the  effects  of  nonlinearities,  such  as  the  Digital  Multimode 
Flight  Control  System  and  servo  rate  limits  not  modeled  in  the 
derivation  of  the  control  laws.  Lsteral  stick  step  input  deflections 
and  release  from  wings  level  sideslips  were  performed  for  each  control 
law  configuration.  In  addition,  each  pilot  practiced  the  open  loop 
inputs  and  became  familiar  with  the  HUD  tracking  tasks  prior  to  flight. 

Open  loop  testa  were  conducted  under  the  direction  of  the  test 
engineer  using  realtime  TM.  For  each  control  law  conf iguraton,  two 
inputs  were  made.  The  first  was  a  lateral  stick  deflection  roll  from  45 
degrees  of  bank  through  level  flight  to  45  degrees  of  bank  in  the 
opposite  direction.  The  test  engineer  monitored  the  stick  input  and 
roll  rate  output  to  insure  a  near  square  step  input  and  an  output 
magnitude  which  allowed  accurate  roll  mode  time  constant  determination. 
The  second  input  was  a  controls  free  release  from  a  wings  level,  steady 
sideslip  of  approximately  three  degrees.  A  precise  lateral  trim  shot 
was  performed  beforehand  to  minimize  trim  effects  in  the  bank  angle 
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response.  The  test  engineer  monitored  the  release  to  ensure  bank  angle 
and  sideslip  angle  response  magnitudes  allowed  accurate  determination  of 
dutch  roll  frequency  and  damping  ratio,  and  phi/beta  magnitude  and  phase 
angle.  Both  open  loop  tests  were  repeated  as  necessary  to  obtain 
quality  inputs. 

The  test  engineer  assigned  specific  closed  loop  test  points  for 
each  mission  based  on  data  requirements.  The  first  teatpoint  of  each 
mission  evaluated  standard  YA-7D  lateral -directional  control  laws.  To 


prevent  biased  data,  the  pilots  did  not  know  which  phi/beta  ratio  they 
were  testing.  Two  HUD  tracking  tasks  and  air-to-air  tracking  task, 
described  in  the  following  section,  were  performed  using  stick  inputs 
only.  After  each  closed  loop  task,  the  pilot  comment  card  in  Appendix  G 
was  completed  and  Cooper-Harper  ratings  assigned  using  the  scale  in 
Figure  G2.  An  A-7D  safety  chase  aircraft  was  used  as  a  target  for  the 
air-to-air  tracking  task.  Each  test  team  pilot  evaluated  each 
configuration  at  least  twice. 


HUD  tracking  Task  1.  The  closed  loop  tracking  task  consisted  of  an 
artificial  target  driven  to  a  lateral  offset  in  a  step  function  by  a 
random  number  generator.  A  centered,  fixed  reticle  consisting  of  five 
and  ten  milliradian  (mil)  diameter  rings  was  also  generated  (Figure  F3) 
The  closed  loop  HUD  tracking  task  was  to  acquire  the  ail  center  of  the 
target  within  the  five  ail  reticle  as  expeditiously  as  possible. 

The  target  moved  after  the  pilot  tracked  the  target  for  one  second.  The 
task  ran  for  90  seconds. 


Desired  Performance:  centar  of  within  tha  fiva  ail  ring 
in  laas  than  3  seconds,  SOX  of  tha  tine,  and  raaainad  within 
tha  ring. 

Adequate  Parforaanca:  cantar  of  within  tha  fiva  ail  ring 
in  laaa  than  6  aaconda,  90x  of  tha  tine,  and  raaainad  within 
tha  ring. 


Figura  F3.  HUD  Tracking  Taak  1  Syabology 

HUD  Tracking  Taak  2.  Thia  cloaad  loop  tracking  taak  conaiatad  of  a 
roll  coaaand  rafaranca  lina  with  tic  aarka  at  ♦  5  dagraa  bank  anglaa 

froa  laval  (Figura  F4) .  Tha  roll  coaaand  lina  waa  drivan  by  a  randoa 
nuabar  to  bank  anglaa  within  ♦  45  dagraaa  of  wing  laval  auch  that  the 
pilot  generated  rata  waa  approxiaataly  40  degraaa/aecond.  Tha  pilot 
attempted  to  aatch  tha  wing  rafaranca  lina  (which  waa  parallel  to  the 
teat  aircraft'a  winga)  with  tha  roll  coaaand  reference  line  over  a  taak 
period  of  90  aaconda. 

Air-to-air  tracking  taak.  Tha  target  aircraft  aatabliahed  laval 
flight  at  15,000  feat  He  and  tha  daaign  Hach  nuabar  for  tha  control 
ayataa  being  evaluated  and  called,  "Target  Ready."  Whan  tha  TN  engineer 
called,  "TH  ready,"  tha  teat  aircraft  called,  "Teat  ready,  aanauvar 
one."  Tha  target  aircraft  continued  atraight  and  laval  for  aaneuvar  one 


Deal rad  Performance:  wings  aligned  within  ♦  5  degrees  of  the 
artificial  horizon  90k  of  the  time. 


Adequate  Performance:  wings  aligned  within  ♦  10  degrees  of 
the  artificial  horizon  90X  of  the  time. 


Figure  F4.  HUD  Tracking  Task  2  Symbology 


and  maneuver  two.  When  the  test  aircrsft  called,  “Test  ready,  maneuver 
three,'*  the  target  aircraft  began  a  series  of  roll  reversals  at  a  random 
interval  between  one  and  three  seconds  at  a  roll  rate  less  then  20 
degreea/aecond.  The  target  did  not  bank  greater  than  40  degrees,  nor 
did  it  exceed  1.5  g'a.  Maneuvering  continued  until  the  teat  aircraft 
called,  "Maneuver  complete,"  or  any  participant  called,  “Knock  it  off." 

Maneuvering  was  conducted  within  the  following  data  band: 


Altitude:  15,000  ♦  1000  feet  He 

Airspeed:  Design  Mach  *-20,  -5  KZAS 

Load  factor:  1.0  -  1.5  g 
Power :  As  required 


A  five  and  ten  ail  ring  reference  mark  was  drawn  on  the  canopy  bow 
of  the  teat  aircraft  approximately  two  inches  aft  of  the  canopy  bow 
along  the  centerline. 


Three  aaneuvers  were  perforaed  in  a  buildup  approach.  The  first 
two  aaneuvers  did  not  receive  Cooper-Harper  ratings,  but  were  designed 
to  expose  poor  handling  qualities  before  aoving  to  a  high  gain  task 
close  to  the  target  aircraft.  The  test  aircraft  positioned  co-speed, 

200  feet  behind,  20  feet  below  the  jetwash,  and  offset  200  feet 
laterally  behind  the  target  aircraft.  When  the  target  aircraft  was 
ready,  the  test  aircraft  called,  “Maneuver  one,"  and  aoved  aggressively 
to  the  six  o'clock  position,  200  feet  aft  of  the  target  while 
aaintaining  longitudinal  and  vertical  spacing.  After  coapleting  the 
aaneuver,  the  test  aircraft  called,  “Maneuver  one  coaplete.“  The  test 
aircraft  then  adjusted  position  to  one  ship  length  the  target,  20  feet 
below  the  jetwash,  and  offset  50  feet  laterally.  On  the  call,  “Maneuver 
two,“  the  test  aircraft  aoved  aggressively  to  position  the  canopy 
reference  aark  on  the  intersection  of  the  target  tailpipe  and  vertical 
stabilizer  while  aaintaining  longitudinal  and  vertical  spacing.  After 
attaining  position,  the  test  aircraft  called,  “Maneuver  three,”  and 
naneuvered  aggressively  to  keep  the  reference  aark  on  the  target 
tailpipe/vertical  stabilizer  intersection.  After  sufficient  tracking 
tiae,  the  test  aircraft  called,  "Maneuver  coaplete.”  The  test  aircraft 
was  liaited  to  leas  than  2  g's  and  was  cautious  to  avoid  the  target's 
jetwaah.  Rating  criteria  during  aaneuver  three  was: 

Desired  perforaance:  within  5  ails  of  intersection  of  target 
tailpipe  and  vertical  stabilizer,  90s  of  the  tiae. 

Adequate  Perforaance:  within  10  nils  of  intersection  of 
target  tailpipe  and  vertical  stabilizer,  90X  of  the  tiae. 


APPENDIX  G.  Miscellaneous  Information 

Sortie  Summary 

A  auaaary  of  the  flight  times  and  personnel  are  given  in  Table  Gl. 


Table  Gl 


Sortie  Suaaary 


Date 

Pilot 

Duration 

Chaae 

Chase  Duration 

8  Oct 

Coatigan 

1.3 

Pitotti 

1.4 

8  Oct 

Gregory 

2.0 

Pitotti 

2.1 

9  Oct 

Gregory 

1.3 

Pitotti 

1.4 

13  Oct 

Mitchell 

1.7 

Pitotti 

1.8 

13  Oct 

Gregory 

1.8 

Waggoner 

1.9 

14  Oct 

Mitchell 

1.7 

Pitotti 

1.8 

14  Oct 

Mitchell 

1.7 

Pitotti 

1.8 

15  Oct 

Mitchell 

1.5 

Waggoner 

1.6 

15  Oct 

Gregory 

1.8 

Pitotti 

1.9 

20  Oct 

Costigan 

1.5 

Pitotti 

1.5 

21  Oct 

Coatigan 

1.5 

Waggoner 

1.6 

21  Oct 

Coatigan 

1.3 

Waggoner 

1.4 

27  Oct 

Gregory 

1.8 

Pitotti 

1.8 

27  Oct 

Mitchell 

1.6 

Waggoner 

1.7 

12  Nov 

Coatigan 

2.0 

Waggoner 

2.1 

13  Nov 

Gregory 

2.2 

McCorry 

2.2 

17  Nov 

Mitchell 

1.6 

Waggoner 

1.7 

Pilot  Experience 

Three  pilots  participated  in  thia  flying  qualities  investigation. 
All  three  were  students  at  the  USAF  Test  Pilot  School  during  the  teat 
program.  Their  backgrounds  were: 


Pilot  A  -  2,000  total  flight  hours  with  1,300  hours  in  the  T-38  and  600 
hours  in  operational  F-4E's.  Operational  experience  was 
priaerily  air-to-air. 

Pilot  B  -  1,200  total  flying  hours  with  1,100  hours  in  the  F-lll. 
Operational  experience  was  priaarily  air-to-ground. 


Pilot  C  -  3,200  total  flying  hours  with  1,400  hours  in  F-4's  and  1,100 

in  F/A-18"s.  Operational  experience  included  both  air-to-air 
and  air-to-ground. 

Pilot  Cards 

The  test  teas  developed  cards  to  be  carried  on  each  mission  to 
standardize  pilot  comments  and  to  facilitate  correlating  data.  The 


pilot  card  is  shown  in  Figure  Gl.  Figure  G2  shows  the  Cooper -Harper 
rating  scale  used  by  the  pilots  to  evaluate  the  aircraft  in  each  task 


« i.t  k>»: 


PILOT  COMMENT 

CARD 

DATE 

MISSION 

PILOT 

CHASE 

Multinode  Control  Panel  Settings: 

i 

i 

i 

i 

i 

i 

i 

<n 

GU 

1 

1 

1 

1 

1 

1 

1 

N 

O, 

P4 

Fuel  Weight _ 

— 

Configuration  (phi/beta) __ 

(filled  by 

engineer) 

TASK  HUD  TASK  1 

HUD  TASK  2 

AIR-AIR  TRACKING 

Cooper-Harper  Rating  1 

2  3  4  5  6 

7  8  9  10 

Did  external  factors  affect  the  teat? 

Yes  No 

Turbulence  Cain 

Slight 

Moderate 

Other 

- - — _ 

Did  the  longitudinal  characteristics  affect  the  teat? 

Yea  No 

Response  characteristics 

Roll  Response 

Sluggish 

Satisfactory 

Abrupt 

Roll  Predictability 

Inprscise 

Satisfactory 

Precise 

Roll  Oscillations 

Negligible 

Noticeable 

Objectionable 

Adverse/Proverse  Yaw 

Negligible 

Noticeable 

Objectionable 

Yaw  Predictability 

Inprscise 

Satisfactory 

Precise 

Pilot  Consents: 

Trackinq  Task  Score 

Engineer  Consents: 

Figure  Gl.  Pilot  Consent  Card 


ADEQUACY  FOR  SElECTEO  TASK  OR 
REQUIRED  OPERATION* 


AIRCRAFT  CHARACTERISTICS  •  *  SEATED  TASK  OR  REQUIRED  OPERATION* 


Is  it 

sonsfoctory  wnhouf' 
v.  improvement // 


/Is  odequal*\ 
/  performance  \ 
QeianoOle  with  a  WeraOKC 
\  OKU  aorklood?/ 


it  cantranat)K/ 


EmcatlSM 

Pilot  compensation  not  o  loctor 

tor 

Highly  dcsiroble 

desired  performance 

Good 

Pilot  compensation  not  a  loctor 

lot 

Negligible  deficiencies 

desired  performance 

Far  -  Som«  mildly 

Minimal  pilot  compensation  requred 

*0f 

unpleosant  deficiencies 

desired  performance 

Pilot  daemons 


Dafmition  of  required  oparotion  involves  designotion  of  flight  phase  and/or  subphosot  with 
accompanying  conditions. 
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